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significance,  but  illustrate  the  types  of  analysis  to  be  performed  when 
these  techniques  are  applied  to  more  data. 

Measured  scattering  coefficient  based  on  instantaneous  radar  antenna 
footprint  was  cross-correlated  with  wave  height  and  slope  for  measurements 
at  the  Noordwijk  tower  on  21  September  1979-.  The  scattering  coefficient 
leads  the  wave  height  by  about  104°  and  leads  the  slope  by  about  77°.  A 
new  modulation  index  concept,  y(t),  the  ratio  of  the  measured  scattering 
coefficient  (based  on  instantaneous  radar  antenna  footprint)  to  the  theore¬ 
tical  scattering  coefficient  (based  on  uniform  capillary  modulation  and 
instantaneous  resonance  with  the  K-spectrum)  has  been  presented.  The  RMS 
value  of  p(t)  is  0.876  with  standard  deviation  0.349  for  HH-polarization  and 
0.673  with  standard  deviation  of  0.044  for  VV-polarization.  The  cross¬ 
correlation  of  p(t)  and  wave  height  shows  that  the  peak  of  the  modulation 
leads  the  wave  height  by  about  153*.  Cross-correlation  of  y ( t )  and  slope 
shows  that  the  peaks  of  the  modulation  are  located  at  about  the  minima 
of  the  dominant  slope,  165®  away  from  both  sides  of  the  maximum.  An 
algorithm  for  transformation  from  temporal  to  spatial  domain  has  been 
developed,  and  the  modulation  transfer  function  based  on  Plant's  definition 
[1979]  has  been  calculated. 
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1.0  ABSTRACT 


Measurements  of  the  radar  cross-section  of  the  sea  were  made  from 
the  Noordwijk  platform  off  the  Dutch  coast  during  September-November 
1979  as  part  of  Project  MARSEN  using  an  FM  microwave  spectrometer  over 
the  8-18  GHz  frequency  range.  A  set  of  techniques  for  studying  the 
modulation  of  the  capillary  waves  on  the  ocean  surface  was  developed 
for  use  in  analyzing  the  entire  observation  set  and  is  presented  here 
as  applied  to  a  single  set  of  observations.  Thus,  the  results  given 
here  are  only  of  limited  oceanographic  significance,  but  illustrate  the 
types  of  analysis  to  be  performed  when  these  techniques  are  applied  to 
more  data. 

Measured  scattering  coefficient  based  on  instantaneous  radar 
antenna  footprint  was  cross-correlated  with  wave  height  and  slope  for 
measurements  at  the  Noordwijk  tower  on  21  September  1979.  The 
scattering  coefficient  leads  the  wave  height  by  about  104°  and  leads 
the  slope  by  about  77°.  A  new  modulation  index  concept,  |i(t),  the 
ratio  of  the  measured  scattering  coefficient  (based  on  instantaneous 
radar  antenna  footprint)  to  the  theoretical  scattering  coefficient 
(based  on  uniform  capillary  modulation  and  instantaneous  resonance  with 
the  K-spectrum)  has  been  presented.  The  SMS  value  of  p(t)  is  0.876 
with  standard  deviation  0.349  for  HH-polarization  and  0.673  with 
standard  deviaton  of  0.044  for  W-polarization.  The  cross-correlation 
of  p(t)  and  wave  height  shows  that  the  peak  of  the  modulation  leads  the 
wave  height  by  about  133°.  Cross-correlation  of  p(t)  and  slope  shows 
that  the  peaks  of  the  modulation  are  looated  at  about  the  minima  of  the 
dominant  alope,  165°  away  from  both  sides  of  the  maximum.  An  algorithm 
for  transformation  from  temporal  to  spatial  domain  has  been  developed. 


and  the  Modulation  transfer  function  based  on  Plant's  definition  [1979] 
has  been  calculated. 
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2.0  INTRODUCTION 


Radar  backacatter  from  the  ocean  can  be  need  to  aonitor  the  wind 
veloeity  with  a  scatteroaeter  that  obaervea  at  two  os  wore  directiona 
relative  to  the  wind  vector,  to  determine  wave  heights  and  ocean 
elevation  (and  from  it  cnrrenta)  with  an  altimeter,  and  to  aonitor 
features  of  the  wave  spectrum,  bottom  topography,  internal  waves,  ship 
wakes,  and  oil  spills  with  an  imaging  radar.  Since  the  radar 
backscatter  represents  a  confusing  "clutter''  for  radars  observing 
targets  such  as  ships  and  periscopes,  studies  of  the  radar  backscatter 
have  been  conducted  since  the  early  days  of  radar  before  the  direct 
applications  listed  above  were  developed. 

The  relatively  recent  advent  of  spaceborne  scatterometers, 
altimeters,  and  synthetic-aperture  imaging  radars  has  led  to  greatly 
increased  interest  in  the  nature  of  radar  backscatter  from  the  sea. 
Furthermore,  the  use  of  these  instruments  calls  for  different  kinds  of 
understanding  (consequently  theory  and  measurements)  than  those 
required  for  earlier  clutter  measurements.  To  aid  in  improving  this 
understanding,  several  major  international  measurement  programs  have 
been  conducted  in  recent  years.  Project  MAR SEN,  of  whioh  this  research 
is  a  part,  is  one  of  the  more  recent  programs.  This  program  was 
conducted  in  the  fall  of  1979  in  two  instrumented  areas  in  the  North 
Sea:  off  the  German  coast  near  the  Helgoland  Bight  and  off  the  Dutch 
coast  near  Noordwijk.  This  is  a  report  of  part  of  one  of  the  studies 
conducted  from  the  Noordwijk  platform  of  the  Dutch  Department  of  Water 
Control . 

Measurements  of  radar  backscatter  from  instrumented  towers  in  the 
sea  are  based  on  the  need  for  observations  under  more  carefully 


controlled  conditions  than  axe  possible  vitb  measurements  from  aircraft 
and  spacecraft,  and  to  study  vitb  a  greater  variety  of  conditions  than 
can  be  encountered  with  measurement*  from  aircraft  vbose  flight  times 
are  limited.  A  goal  of  measurement  programs  such  as  MAE SEN  is  to 
perform  many  measurements  of  different  kinds  from  towers  and  aircraft 
so  that  the  synergism  of  data  comparisons  may  lead  to  better 
understanding  than  vould  be  obtained  vitb  isolated  measurement  programs 
by  small  groups  of  investigators. 

Radar  measurements  were  made  at  the  Noordvijk  platform  by  three 
groups:  The  University  of  Eansas,  a  Dutch  group  of  experimenters  from 

the  Physics  Research  Laboratory  of  National  Defence  and  the  Technical 
University  of  Delft,  and  a  French  group  from  the  Institut  Francaise  du 
Petrole.  A  JPL  group  made  measurements  with  a  capillary  wave  gauge  and 
near-surface  vind  probes,  and  a  British  group  made  measurements  of 
breaking  waves  with  special  buoys  deployed  from  the  platform.  In 
addition,  overflights  were  made  by  a  NASA  aircraft  carrying  a 
scatterometer  and  by  several  imaging  radars  from  the  U.S.,  Germany, 
France,  and  the  Netherlands.  The  object  of  the  aircraft  measurements 
over  the  Noordvijk  platform  was  to  allow  comparison  with  the  platform 
measurements. 

Major  goals  of  the  University  of  Eansas  measurements,  most  of 
which  were  shared  by  the  other  radar  investigators  on  the  platform 
were: 

1.  To  establish  the  variation  of  the  capillary-wave  amplitude  over 
the  underlying  larger  waves,  and  to  establish  the  portion  of  the 
underlying  waves  contributing  the  strongest  radar  signals.  This 
is  important  because  the  nature  of  a  synthetic-aperture  image  of 
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tli*  larger  waves  depends  on  the  local  velocities  associated  with 
the  strongest  scattering  centers,  and  these  velocities  are 
different  for  different  parts  of  the  larger  waves.  The  nature  of 
the  University  of  Kansas  radar  allowed  it  to  be  used  as  a  ''radar 
wave  gauge"  to  measure  the  wave  height  at  the  center  of  the  radar 
beam,  so  the  type  of  analysis  conducted  to  meet  this  objective  was 
unique — no  other  radar  measurements  have  been  made  where  the  wave 
height  was  measured  exactly  at  the  radar  observation  spot. 

To  establish  the  variations  of  the  capillary-wave  component  of  the 
ocean-wave  spectrum  for  different  wind  speeds  and  directions 
relative  to  that  of  the  wind.  This  is  important  because  the  radar 
backscatter  from  the  sea  at  angles  of  20°  or  more  from  vertical  is 
known  to  be  governed  by  resonance  with  these  capillary  and 
near-oapillary  waves.  In  faot,  the  determination  of  the  spectrum 
from  the  radar  measurements  is  based  on  this  known  resonance 
phenomenon.  Tith  the  University  of  Kansas  radar's  capability  of 
measuring  over  wavelengths  from  about  1.7  cm  to  3.3  cm  and  that  of 
the  French  radar  from  3.3  cm  to  25  cm,  this  experiment  provided  an 
unparallelled  opportunity  to  observe  the  capillary  spectrum. 

To  determine  empirically  the  variation  of  radar  scattering 
coefficient  from  the  sea  over  a  vide  range  of  angles  of  incidence, 
angles  relative  to  the  wind  direction,  frequencies,  and 
polarisations.  This  is  importsnt  in  scstterometry,  radar  design, 
and  clutter  rejection  studies.  Although  many  previous 
measurements  of  this  nature  have  been  performed  from  aircraft  and 
stationary  platforms,  none  has  allowed  measurements  with  such  a 
variety  of  radar  parameters  at  one  time.  Furthermore,  although 


aircraft  and  spacecraft  aeasureaents  have  allowed  developaent  of 
sncceaaful  wind-vector  aeasureaent  ays teas  using  scattoroaeters, 
one  of  the  greatest  unknowns  is  the  exact  relationship  between 
scattering  coefficient  and  wind  speed — the  previous  aeasureaents 
have  led  to  considerable  seatter  in  the  values  of  the  coefficients 
in  this  relationsip. 

4.  To  establish  the  variation  of  the  quantities  determined  under 
objectives  1-3  when  an  oil  spill  is  present  on  the  surface.  Oil 
is  known  to  daap  the  shorter  coaponents  of  the  ocean-wave 
spectrum,  and  this  fact  is  used  as  the  basis  for  oil-spill 
aonitoring  radar  systeaa  naed  by  several  countries,  although 
quantitative  aeasureaents  of  the  effect  are  scarce. 

The  study  reported  here  is  related  to  object  (1)  above,  the 
determination  of  the  modulation  of  the  capillary  waves  associated  with 
the  underlying  larger  waves.  This  report  is  priaarily  concerned  with 
developaent  of  the  aethodology  for  a  later  study  of  the  aore  complete 
data  set.  However,  the  preliminary  results  reported  for  a  single  set 
of  data  have  value  in  their  own  right. 

One  of  the  approaches  coaaonly  used  for  describing  this  modulation 
is  to  calculate  a  ' 'modulation  transfer  function' '  (KTF)  in  the 
frequency  domain.  The  KTF  amplitude  shows  the  underlying-wave 
frequencies  at  which  the  capillary  aaplitude  correlates  best  with  the 
underlying  waves,  and  the  phase  indicates  the  location  of  the  relevant 
capillary  ooaponents  on  the  individual  frequency  coaponents  of  the 
underlying  waves. 

In  this  report  an  alternate  tiae-doaain  approach  is  presented  that 
is  based  on  the  crosa-oorrelation  function  between  the  radar  signal  and 


til*  underlying  mu.  The  eross-correlaton  function  tends  to  show  the 
principal  effects  in  a  nore  eoapaet  and  understandable  fora  than  the 
MTF.  So  that  this  relation  nay  be  understood  better,  the 
cross-correlation  function  is  calculated  for  the  capillary  swdulation 
alone  by  using  as  an  input  the  ratio  of  the  observed  signal  to  the  one 
that  would  be  observed  if  the  capillary  vaves  were  of  unifora  aaplitude 
(changes  in  the  slope  of  the  underlying  wave  would  cause  all  variations 
in  the  received  signal  if  this  were  true).  This  " instantaneous 
aodulation' *  can  only  be  calculated  in  the  tine  doaain. 

During  the  course  of  the  investigation  of  the  data  the  presence  of 
strong-signal  ''events''  was  observed.  These  are  presuaably  the  saae 
as  the  ''sea  spikes'*  that  have  been  reported  by  nuaerous  previous 
investigators.  To  deteraine  the  capillary-wave  effects,  one  aust  first 
reaove  these  anomalously  large  signals  froa  the  time  series  of  the 
signal.  Methods  for  doing  this  have  been  devised,  and  the  result  of 
the  reaoval  is  a  significant  iaproveaent  in  the  repeatability  of  the 
nodulation-study  results. 

This  study  was  intended  priaarily  to  develop  aethodology,  so  only 
one  data  set  was  used — this  was  for  a  wind  speed  of  8.8  a/ sec  with  the 
radar  looking  in  the  upwind  direction.  The  doainant  wave  period  was 
5.8  seoonds,  and  the  aaziaua  of  the  aeasured  scattering  coefficient  was 
found  to  occur  1.8  seoonds,  or  104°,  before  the  crest  of  the  onooaing 
wave — it  is  therefore  on  the  front  face  of  the  wave  and  nearer  the 
trough  than  the  crest.  The  loostion  of  the  aaziaua  instantaneous 
aodulation  (ratio  of  observed  signal  to  that  predicted  if  the  capillary 
waves  were  uniformly  distributed)  is  153°  ahead  of  the  crest,  a 
location  alaost  in  the  trough.  Furthermore,  a  secondary  aaziaua  was 
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found  at  62°  behind  the  crest,  or  veil  np  on  the  back  face  of  the  wave. 

Correlations  were  also  perfomed  between  the  scattered  signal  and 
the  instantaneous  slope  of  the  waves.  The  slope  spectrum  has  a  higher 
peak  frequency  than  the  wave-height  spectrum,  so  these  correlations 
lead  to  different  separations  between  maximum  of  the  slope  and  aaxiaua 
radar  signal  than  would  be  observed  if  a  single  sine  wave  were 
observed.  That  is,  one  should  be  able  to  obtain  the  location  of  the 
peak  signal  on  the  slope  of  a  single  sine-wave  component  from  the 
height  record  alone,  but  with  a  coaplex  wave  pattern  this  is  no  longer 
possible.  Interpretation  of  the  relation  between  peak  radar  scatter 
and  slope  remains  to  be  worked  out. 

Most  of  the  analyses  were  perfomed  on  tine  records,  since  the 
wave  heights  were  measured  at  the  fixed  center  of  the  radar  footprint 
as  time  histories.  However,  one  oan  determine  the  spatial  picture  over 
a  limited  distance  from  a  time  record  by  suitably  summing  a  series  in 
distance  using  the  Fourier  coefficients  of  the  time  series.  The  method 
for  doing  this  was  also  developed  here  and  will  be  used  in  future 
studies  of  the  measurements. 
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Hi*  phenomenon  of  sen  scatter  and  its  theoretical  basis  has  been 
of  interest  since  the  advent  of  the  practical  radar  systaa.  It  has 
been  generally  agreed  [Vright  (1968),  Chan  and  Fang  (1973).  Va  and  Fang 
(1972),  Long  (1974)1  that  the  scattering  propertiea  of  the  sea  sarfaee 
can  be  explained  fairly  well  in  tens  of  a  two-scale  slightly  roagh 
sarfaee  aodel,  where  the  small-scale  capillary  waves  are  assnmed  to 
satisfy  the  small  pertarbatioa  assumption,  and  the  large-scale  wavea 
are  assamed  to  satisfy  the  Kirchhoff  approximation.  This  theory, 
verified  to  a  large  extent  by  experiment,  indicatea  that  baekseatter  at 
microwave  frequencies  is  from  those  wavelength  components  of  the 
ocean-wave  spectrum  that  are  resonant  to  the  radar  probing  wavelength 
(Figare  3.1).  In  general,  the  acattering  cross-section  per  unit  area, 
<S°,  depends  on  the  local  angle  of  incidence  (Figare  3.2),  and  also  on 
the  small  capillary  waves  which  are  tilted  by  the  large-scale  wave 
slope.  This  angle  depends  on  the  position  of  the  eapillary  waves  with 
respect  to  the  large-scale  waves.  The  amplitades  of  the  small 
capillary  waves  are  modified  by  a  variety  of  hydrodynamic  interactions 
with  the  large-scale  waves.  The  local  scattering  eross-section  depends 
on  this  aaplitade.  In  fact,  in  first-order  scattering  theory,  the 
cross-section  and  power  speetram  aaplitade  are  proportional (Sect ion 
4.3.1).  Therefore,  it  is  the  modalation  of  scattering  by  the 
large-scale  waves  whioh  makes  them  detectable  by  microwave  radar.  If 
the  small  capillary  wave  amplitade  is  modalated,  so  is  the  stress,  k 
modalated  stress,  if  of  the  proper  phase,  resalts  in  larger 
wave-growth.  Table  3.1  shows  the  list  of  symbols  for  the  following 
analysis. 
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Wind  Spetd  (  m/sec  ) 

Figure  .VI:  Comparison  between  theoretical  and  measured  app  for  vertical 
and  horizontal  polarizations  at  40°  incidence  angle  for  down¬ 
wind  observation  with  a  13-9  GHz  radar.  This  illustrates 
close  correspondence  between  experiment  and  theory  based  on 
resonance  with  capillary  waves. 


Figure  3.2:  Dependence  of  scattering  coefficient  on  local  angle  of 
incidence,  9'  and  small  capillaries  which  are  tilted  by 
the  large-scale  wave  slope. 
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TABLE  3. 1 
LIST  OF  SYMBOLS 


u*  Ai r  friction  ve loci ty 

0  Angle  of  incidence  (angle  between  vertical  and  radar-wave 

di rection) 

kx  Component  of  wave  number  in  x  direction 

cg,cg°  Group  speed  of  wind-generated  wave  with  and  without  plunger¬ 
generated  wave 

u  Horizontal  component  of  orbital  velocity  of  plunger-generated 

wave 

e(t)  Linearly-detected,  rectified,  received  signal 

0'  Local  angle  of  incidence  (angle  between  normal  to  local 

surface  and  radar-wave  direction) 

M  Measured  fractional  modulation 

r  Measured  peak-peak  modulation  of  R(0) 

kQ  Microwave  number 

P  Modulation  index  of  received  signal 

uQ  Modulation  of  u 

u(t)  Periodic  constituent  of  e(t) 

f,(k)  Perturbation  of  surface  displacement  spectrum 

Phase  angle 

Phase  speed  of  plunger-generated  wave 
0  Phase  speed  of  wind-generated  wave  in  absence  of  plunger¬ 

generated  wave 

Plunger-generated  wave  number 
Plunger-generated  wave  slope 
Radian  frequency  of  plunger-generated  wave 
Radian  frequency  of  wind-generated  wave 
(t)  Random  (wind-generated)  constituent  of  e(t) 

Received  power 

0  Received  power  in  absence  of  plunger 

Straining  constant 
Straining  function 
Straining  function 
Time  lag 

Wind-generated  wave  growth  rate 
r  Wind-generated  wave  relaxation  rate 

Wind-wave  interaction  functional 
,Fq  Wind-wave  surface  displacement  spectrum 


Keller  and  Wright  [1975]  tnggaatad  a  two-seale  seattaring  theory 
called  'the  relaxation  tine  nodal,'  and  state  tnat  tne  snail  capillary 
waves,  if  perturbed  fron  equilibrium,  will  relax  bach  at  an  exponential 
rate  called  'the  relaxation  rate.'  The  relaxation  rate  is  a  function  of 
the  snail  capillary-wave  spa ct ran. 

They  found  that  these  waves  grew  exponentially  over  several  orders 
of  swgnitude  in  spectral  intensity.  The  energy  influx  fron  the  wind  to 
a  wave  of  given  wave  nunber  nuat  be  balanced  by  transfer  to  other  wave 
nnnbers  for  dissipation  or  growth  .of  the  wave.  The  transport  equation 
gives  this  balance.  In  the  case  of  growth  along  a  single  cartesian 
coordinate  x  in  the  direction  of  the  wind,  it  nay  be  written: 


2£+  Ca  i_EL  -  J!F+ 

i-t 


(3.1) 


where  F(k,x,t)  is  the  surface  displacenent  s pec t ran,  eg  is  the  group 
velocity,/^  is  the  exponential  growth  rate,  H(F,k.u*)  is  a  function 
which  accounts  for  nonlinear  transfer  to  other  waves  and  non-linear 
dissipation,  all  evaluated  at  wave  nunber  k.  The  quantities  E~  and 
F  nay  be  obtained  directly  in  wave-tank  experinents  fron  the 
intensities  of  the  first-order  Bragg  line  in  the  Doppler  spectra 
(Figure  3.3a).  ft  is  the  initial  tewporal  growth  rate,  provided 
nonlinear  interactions  are  negligible  in  the  initial  stages  of  growth 
[Keller  and  Wright,  1975].  The  quantities  and 

the  naxinun  value  of  F  observed 


(C9C*£) -/fifty'  u 

rv.  8  '  b  ¥  ttm 

under  the  conditions  of  the  neasureaent 

soae  cases  reported  by  Plant  and  Wright  [1977J.  The  differenced)  C '  )-C 

J  i 

- 


under  the  conditions  of  the  neasureaent,  are  shown  in  Figure  3.3b  for 

VE' 

can  be  aatehed,  for  the  cases  shown,  by  values  of 
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^nl/F«*jt  calculated  fro*  the  second-order  gravity-capillary  wave— vave 
interaction.  Then  the  nonlinear  interactions  bring  tbe  vave  to  a 
steady  state  or  equilibria*,  tbe  vertical  flax froa  tbe  vind  is 
large  coapared  to  horizontal  flax  *>(£)•  Tbe  difference  is 
dissipated  either  by  tbe  breaking  of  tbe  vave  itself  or  by  tbe  transfer 
to  other  vave  nnabera  for  subsequent  dissipation,  but  in  either  case, 
tbe  equilibria*  is  spatially  localized.  It  is  tbe  perturbation  of  this 
equilibria*  by  tbe  fluid  notions  associated  vitb  tbe  large  vaves  vhicb 
aodulates  tbe  saall-vave  aaplitade.  and  in  consequence,  tbe  scattering 
cross-section. 


wmo  wave  ooem.cn  spectra 


Figure  3.3a 

Evolution  of  Doppler  spectra  of  backscattered  microwaves 
with  increasing  wind  speed  at  fixed  fetch.  The  wavelength 
was  6.9  cm  (f  »  4.35  GHz),  the  depression  angle  was  30°, 
and  the  Bragg  wavelength  was  4.05  cm. 


'  +■  |j||  '  i 
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Figure  3.3b 

Terms  in  the  equilibrium  transport  equation  for  gravity-capillary 
waves.  The  curve  labeled  Sni/Fm  is  the  difference  of  the  other 
two  curves,  i.e.,  Snj/Fm  =  (Cg (3F/3x)- 2F)/Fm. 

The  Bragg  wavelength  is  denoted  by  Xg. 


A  phenomenological  nodel  of  the  perturbation  of  this  equilibrium 
by  the  horizontal  component  of  orbital  speed  of  a  large  nave  [Keller 
and  Vright.  1975]  can  be  constructed  after  including  the  etfect  of  the 
horizontal  current  u(x,t)  in  the  transport  equation 

+c.  to (3.2) 

U>-  CJ»+  &  (3.3) 
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where: 


Cq  -  bjO 


(3.4) 


The  quantity  is  the  angular  freqency  of  the  saall  wave  under 
consideration.  The  tern  on  the  right  containing  (k),  the  straining 
constant,  is  called  the  radiation  stress  and 


tit)  =  i-  <r9‘/c 


(3.5) 


o  o 

where  and  £  are  the  group  and  phase  speed,  respectively,  of  the 

saall  wave  in  the  absence  of  the  current  u. 

Let  F0(k,x)  be  the  surface  displaceaent  spectrua  in  the  absence  of 
the  current  u(x,t)  which  we  now  assuae  to  have  the  fora 


(3.6) 


This  is  the  current  which  perturbs  the  equilibriua  spectrua.  In 
consequence  we  expand  the  resultant  surface  displaceaent  spectrua  in  a 
perturbation  series  in  (uq/c),  where  c  *J^*/k  is  the  phase  speed  of  the 


large  wave. 


r=r*+LUo/c)fi(t)c  „ . 


(3.7) 


The  key  assuaption  is  that  tne  increaent  in  the  nonlinear  transfer 
can  be  written: 


s  ft*  C  s%r)  (  «•/£)  A 


(3.8) 


v..r-v,=  iwatiStetia  » 
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This  is  called  the  relaxation  time  appproximatioa.  It  is  now 
straightforward  to  solve  (3.2)  aad  (3.7)  to  first  order  ia^ig/c^ 

L  ,dM<L 

/C  '*  (jD-//fr)z 

T  „„ 

"  c- 

where  the  relsxstioa  rate 

fir-  (  ^/S?)  ~  fi  <3.io) 

Ar> 

has  beea  iatrodaced  [fright,  1977]. 

The  first  order  ia  (uq/c)  teras  ia  the  expaasioa  adequately 
describes  the  respoase  of  the  wind-generated  waves  for  wiads  of  less 
thaa  7  to  8  mj sec  aad  ug/c  <  0.1. 

Modalatioa  of  backscattering  cross-sectioas  have  beea  aee«o»£a  at 
9.35  GHz  aad  depressioa  aagle  of  45  degrees  (Bragg  waveleagth  *  2.3  ca) 
by  Cellar  aad  fright  [1975,  1976].  The  liaearly  detected  sigaal, 

proportioaal  to  the  scattered  electroaagaetic  field  was  ratified  aad 
filtered  to  reaove  frequency  a  aad  phase  modalatioa  aad 
cross-correlated. 

The  aeasared  auto-covariance  function,  R(^5  exhibited  a  spike  at 
«  0,  due  to  the  longer  short-gravity  waves  of  the  wind-generated  wave 
system,  a  mean  value  proportional  to  tae  square  of  the  mean  scattered 
field  aad  a  periodic  portion  due  to  modulating  plunger-geaerated  wave: 
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ect)  =  /Ct)  +  ffatt) 


A  1.0).  / \t ) 


(3.11) 


(3.12) 


(3.13) 


it  cut)  =  Cos  (.-ci-e) 


(3.14) 


K.(x)=j.  J  '^(t)Ati-r)oU i-£\j)lc.i(nio  <3.1 

'•  j.  V  ' 


The  ratio  of  the  peak-to-peak  amplitude  of  the  periodic  portion  of 
R(Z*),r,  to  11(0),  which  latter  quantity  is  proportional  to  the  aain 
scattered  power  is  called  the  fractional  modulation,  ■  (Figure  3.4). 
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(3.16) 
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.here /  i  is  friction. 1  ,04011110,. 


(3.17) 
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Figure  3.^ 

Autocorrelation  functions  of  the  microwave  return 
(left)  and  wind-wave  spectra  (right)  for  various  wave 
amplitudes  at  u*  *  16.5  cm  sec" I 


If  we  represent  p,  power  return  as 


P=  P.  (i+  /m 


(3.18) 


wliere 


/W=  (B  +  0%-  T*+  ZOt)(4*/c) 


end 

✓=  r*n-'(<°sry3) 


(3.20) 


(3.21) 


Positive  end  negative  signs  refer  to  npwind  end  downwind-looking 
antenna  orientations,  respectively. 

Tie  calculated  theoretical  curves  and  fractional  modulation,  are 
nonetheless  an  excellent  fit  to  the  measured  value  for  winds  of  less 
than  8  ml  sec  and  ug/c  <  0.1.  For  these  winds,  larger  wave  aaplitndes, 
the  fractional  modulations  saturated,  i.e.,  they  ceased  to  increase 
with  increasing  modulation  wave  amplitude,  (Figures  3.5a,b). 


Figure  3.5a 

Fractional  modulation  vs.  air  friction  velocity  for 
0.575  Hz  wave,  UQ/C  «  0.091;  solid  circles,  upwind;  open 

circles,  downwind 


Figure  3.5b 

Comparison  of  theoretical  and  measured  fractional  modulations 
looking  upwind  vs  orbital  velocity  of  0.575  Hz  waves.  Solid 
data  points  and  solid  line  are  for  u*  ■  16.5  cm  sec*1. 

Open  data  points  and  dashed  line  are  for  u-  *  30  cm  sec"'. 
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Keller  tad  Wright  [1976]  stated  that  at  tha  higher  wind  spaed 
there  is  definitely  another  phenomenon  present  unaccounted  for  by  the 
relaxation  theory.  They  believe  that,  at  present,  the  taeory  is 
phenomenological  and  the  relaxation  rate  is  an  adjustable  parameter. 
The  reasons  for  choosing  certain  values  for  this  and  other  incompletely 
known  parameters  become  clearer  upon  consideration  of  relaxation  of  a 
pair  of  waves  which  simulate  short  gravity  waves  coupled  to  capillary 
waves  by  the  second-order  resonant  interaction  [Valenzuela  and  Laing. 
1972.  Plant  and  Wright.  1977.  Plant  1979]. 

In  a  consequent  study  of  sea  surface  backscatter  phenomena.  Plant, 
Keller  and  Wright  [1978]  conducted  an  experiment  to  measure  the 
modulation  of  small  capillary  wind-generated  waves  by  those  with 
typical  ocean  wave  periods.  They  stated  that  tne  horizontal  component 
of  the  orbital  velocity  of  shoaling  waves  much  exceeds  the  vertical 
component  exception  in  the  region  near  the  wave  crests.  Hence  the 
relative  importance  of  smdulation  due  to  tilting  is  much  reduced  over 
most  of  the  large-scale  wave,  so  the  shoaling  waves  beyond  the  surf 
zone  provide  easily  understood  modulating  conditions.  The  power 
spectra  of  e(t),  rectified  scattered  field  strength  e(t),  and  the 
line-of-  sight  speed  v(t),  here  show  contributions  both  from  the  crests 
and  low-amplitude  qua si-periodic  portions  of  the  waves.  The  narrow 
peaks  in  the  spectra  were  used  to  obtain  quantitative  measurements  of 
the  modulation.  They  also  were  used  to  calculate  'fractional 
modulation/^'  (Figures  3.6a,b  and  c). 

In  the  most  recent  work  by  Plant  et  al.  [1980]  the  modulation 
transfer  function  is  calculated. 

Let  p(t),  the  square  of  the  homodyne-detected  [King,  1978] 
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■iemiTf  voltage,  which  in  torn  is  proportional  to  the  backscattered 
electromagnetic  field,  ba  oallad  tba  iastaataaeoasly  received  povar. 
Lot  ▼( t)  ba  tke  iastaataaaoas  line-of-sight  spaad.  If  6pT  is  tba 
oross-speetrna  of  iastaataaeoasly  received  power  sad  iastaataaaoas 
liae-of-sight  spaad  sad  is  tba  satospactraa  of  tba  latter,  define 


/yy)  ('-G. )  =  (■&•) 


(3.23) 


Tba  pbasa  spaad  CCA)  sad  aasa  received  power  P  normalize  tba 
aodalatioa  transfer  faactioa,  s(A)  asking  it  diaansionless.  GpV(A) 
is.  of  coarse,  coaqplex.  To  fix  ideas,  sappose  tae  plaae  of  incidence 
(tba  plane  containing  tba  aicrowava  propagation  vector)  and  tba  aoraal 
to  tba  and is tar bad  ocean  snrface.  also  contain  the  direction  of  a 
aonocbroaatic,  aodnlating  sarfaca  wave  of  freqaency.^^  and  wave  nnaber 
E.  propagating  in  water  of  depth  D,  sad  the  depression  angle  is  0  , 
than: 


VU -UVo  CnCG-t  +  fs) 


(3.24) 


(f,)m  tcLirh  (kO)  Hvr (ff) 


(3.25) 


vVu,  =  **4  \k»)) 


(3.26) 


and 
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1 


(3.27) 


P(,  +  »"(«'/<; )  Cot  (nti-f)') 

fhtr,  9S  u  tut  phase  U|l*  by  which  the  lia.-of-.ight  .peed  leed.  the 
horizontal  coaponent  of  orbital  .peed,  a-ao  (At)  (Figaro  3.7). 


Figure  3.7 

Modulation  spectra  of  the  demodulated  return  of  the  9.375  GHz 
radar.  Wind  speed  =»  7.5  .m/s. 

(a)  Line-of-sight  component  of  orbital  speed; 

(b)  phase  of  the  modulation  transfer  function; 

(c)  coherence; 

(d)  modulus  of  the  modulation  transfer  function  divided  by  the 
ocean  wave  phase  speed,  m/c. 

The  transfer  faaotioaa  tad  their  phase  shown  ia  Figaros  3.8  to 
3.11  are  referred  to  the  liao-of-sight  Telocity  as  defined  by  (3.24). 
The  transfer  faaotioa  which  refers  to  the  horizontal  coaponent  of 
orbital  speed  (a  aore  asefal  quantity  for  interpretive  purposes)  is 
obtained  by  aaltiplyiag  a  by  Vq/uq  given  by  (3.26).  Although  the 
transfer  faaotioa,  before  averaging  over  the  wind  speed  bands, 
exhibited  scatter  as  large  as  a  factor  of  2,  the  finally  averaged 


■odulcs  (Figure  3.12a)  aud  phasa  (Figure  3.12b)  show  saooth.  stable 
dependence  on  nodulating  frequency  and  wind  speed. 


Figure  3.8 

Modulus  of  the  1.5  GHz  (Bragg 
wavelength  »  13  cm)  transfer  functions. 


Figure  3-9 

Phase  of  the  1.5  GHz  (Bragg  wavelength  «  13  cm)  transfer 
functions.  Phases  are  measured  from  wave  crests  and  are 
positive  leading  the  crest.  Oashed  lines  are  theoretical 
predictions  assuming  modulation  caused  solely  by 
straining.  [Plant,  Keller  and  Wright,  1980] 


Although  the  total  averaged  MTF  looks  very  saooth  and  corresponds 
reasonably  well  with  the  theory,  we  believe  the  JCTF  concept  is 
especially  difficult  for  radar  t&eorists  who  are  not  involved  in  pure 
theoretical  hydrodynamics  and  fluid  aechanics  equations.  Ve  think  that 
'aodulation  index'.  £  £ )  ,  (Section  4.4.2) ,  would,  if  it  could 
be  treated  well,  be  an  easily  interpreted  concept  that  could  be  used  to 
aake  a  siaple  stateaent  about  phase  relationships  between  tne  peak  of 
capillary  and  large-scale  waves. 
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Figure  3-12 

Wind  speed  dependence  of  9.375  GHz  modulation 
transfer  functions,  (a)  NOSC  tower;  (b)  Nags 
Head,  North  Carolina  [Plant  et  al.,  1978], 
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A  distinction  must  bo  made  between  the  idealized  classical  wave  motions 
of  hydrodynamics  and  the  mnch  more  complex  (and  far  less  well 
understood)  wave  motions  that  occur  in  nature.  There  is  s  distinction 
between  deep  water  and  shallow  water  waves  depending  on  the  ratio 
between  wavelength  and  water  depth. 

4.1.2  Definition  of  Periodicity 

The  solution  of  the  classical  equations  for  ocean  waves  takes  the 
form  of  a  sum  of  periodic  functions  (Section  4.1.3).  A  function  f(t) 
is  periodic  of  period  T  if  and  only  if  f(t)  *  f(t+T)  for  all  t. 
Records  of  the  rise  and  fall  of  the  water  surface  as  a  function  of  time 
at  a  fixed  point  on  the  ocean  surface  do  not  satisfy  the  above 
definition.  Water  waves  in  nature  are  not  periodic.  Hence,  the  sum  is 
in  principle  infinite  and  representable  by  a  continuous  spectrum. 

For  all  types  of  truly  periodic  traveling-wave  motion 

^  -CT  (4.1) 


where  is  the  distance  between  the  crests  of  the  waves,  C  is  the 
speed  with  which  the  crests  move  and  T  is  the  time  required  for  the 
water  surface  to  go  through  one  complete  oscillation  at  a  fixed  point. 
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The  period  end  length  of  a  wave  are  rather  cumbersome  notationa 
and  it  is  possible  to  replace  then  by  the  concepts  of  wave  number  and 
angular  frequency. 


k-  ZJL 

y 

=  ?zJL 
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4.1.3  General  Gravity.  Larse-Scale  Wave  Theory 

The  simplest  wave  motion  is  that  of  a  simple  harmonic  long-crested 
traveling  wave  in  the  +x  direction  whose  equation  would  be  given  by 

hL%  t)  s  a  CosC  22 *)  u-4) 

7- 


It  is  possible  to  write 


h  L/x)t)s  a  Cot  (4.4a) 

In  the  following  derivation  it  is  assumed  that  the  earth  is  flat, 
that  the  water  is  of  constant  depth  h,  that  the  Coriolis  force  is 
negligible,  that  the  density  of  the  air  can  be  neglected  compared  with 
the  density  of  the  water,  that  the  water  is  of  constant  density,  that 
the  body  of  water  is  infinite  in  extent  and  completely  covered  by 
waves,  and  that  the  viscosity  and  surtace  tension  can  be  neglected. 
Moreover,  it  is  assumed  that  there  is  no  variation  in  wave  propagation 
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in  the  y  direction.  The  equations  are  unsolvable  in  closed  font,  as 
they  are  nonlinear.  A  perturbation  technique  must  be  employed  to 
linearize  these  equations  and  to  obtain  a  useful  revalt. 

As suae  that  the  unknown  functions  in  hydrodynamic  equations  can  be 
expanded  in  powers  of  a  small  dimensionless  parameter  £.  ,  so  that  the 
solution  is  a  power  series  in  £  that  will  converge  if  enough  terms  are 
taken.  These  expansions  are  substituted  into  the  equations,  and  all 
terms  to  the  same  power  in  £  are  collected.  It  is  possible,  in 
general,  to  obtain  terms  to  zero  power  in  £  ,  to  the  first  power 
in  £  ,  to  the  second  power  in  ,  and  so  on.  Depending  upon  tne  power 
of  £  used  in  the  expansion,  successively  higher  and  more  nearly 
correct  approximations  to  the  solution  of  the  complete  equation  are 
obtained  [Neumann  and  Pierson,  1966] . 

The  steady,  time- independent  wave  equation  can  be  shown  to  be 

3aC,ttkr)-e££&i!®fa»*Ui*)mO  (4.5) 

where 

C*  =,  JL  &wh  Ckh)  (4.6) 

k 

where  g  is  the  gravitational  acceleration  (9.8  m/sec). 

Equation  (4.6)  relates  the  phase  speed  C  to  the  wsve  number  K  and 
the  depth  of  the  water  h.  For  a  given  depth  there  is  only  one  possible 
value  for  the  wave  speed  given  the  wavelength  (or  the  wsve  number). 
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Equation  (4.6)  can  ba  simplified  for  two  situations  with  reference 
to  the  values  of  h  and  E.  The  hyperbolic  tangent  of  a  function  varies 
between  zero  and  1  for  positive  values  of  the  argnaent.  The  hyperbolic 
tangent  approaches  the  value  1  very  closely  for  values  of  the  argnaent 
that  are  very  large  Thus,  if  the  product  of  E  and  h  is  such  that  the 
hyperbolic  tangent  is  essentially  1,  (4.6)  reduces  to  tne  auch  siapler 
fora  given  by 


c 


Z 

o 


(4.7) 


The  waves  that  satisfy  this  siaplified  fora  of  the  equation  are  called 
short  waves,  or  'deep  water'  waves.  The  second  approziaation  occurs  if 
the  product  of  E  and  h  is  very  saall.  Under  tnese  circuastances  the 
hyperbolic  tangent  of  Eh  is  approximately  equal  to  Eh  and  if  this 
substitution  is  aade,  (4.6)  will  be 


C*.3h 


Waves  that  satisfy  this  equation,  relating  the  phase  speed  to  tne 
depth  since  the  length  of  the  wave  no  longer  determines  the  phase 
speed,  are  called  long  waves,  or  'shallow-water'  waves. 

Table  4.1  shows  the  various  relations  for  deep-water  waves  between 
the  different  paraaeters. 
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Table  4.1 

Relationship  between  wavelength,  phase  speed,  period,  angular  frequency 
wave  nnaber,  and  group  velocity  for  deep  water  linear  periodic 
gravity  waves.  Read  down  to  express  the  quantity  in  the  top  row  in 
terns  of  the  quantity  in  the  left-hand  colnan 
[Nenaann  and  Pierson  (1966] 


4.1.4  General  Solution  for  Larte-Scale  Gravity  Waves 
The  linear  solution  that  has  been  obtained  can  be  used  as  a  basis 
for  second-order  solutions.  One  could  proceed  to  tnird  and  even  higher 
orders.  In  general  one  can  show 


h Lot )  =  Z  Hi  <T»s(£Vir  +  9c)  «... 

L 

It  is  possible  to  transform  the  solution  to  the  form  in  which  the 
presence  of  a  current  to  produce  the  steady  notion  is  elininated.  The 
general  solution  for  the  aoving  wave  profile  would  then  be  given  by 


=•  STa*-  Cos(/Ci  (/x+C;t)+%)  (4.io) 

It  is  worthwhile  to  nention  that  Cj  in  general  is  not  equal  to 

.  Cj  and  CQ  have  nonlinear  relations,  but  it  is  very 
difficult  to  find  the  exact  nonlinear  relationships.  In  the  following 
analysis  we  assume  a  linear  relationship  between  C$  and  C0,  where: 
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Ci  ■  c0 


(4.11) 


4.2  Capillary  Imi 

4.2.1  Capillary  Phenomena 

If  the  surface  taaaioa  of  the  water  ia  important,  a a  it  ia  for 
very  abort  waves,  additional  restoring  force  dee  to  corwatare  of  the 
free  surface  eaters  the  Bernoulli  equation. 

Ia  linearized  theory,  the  equation  for  the  phase  speed  is: 

C.  =  (  J-+  Tk.)'/l 


where  ^  is  the  sea  water  density.  Perturbation  techniques  were 
applied  to  the  problem  by  Pierson  and  Fife  [1961],  Ia  general,  waves 
have  a  phase  speed  that  increases  with  increasing  wave  amplitude. 
However,  in  linear  theory  the  speed  and  amplitude  are  independent. 

When  capillary  waves  interact  with  large-scale  waves  the 
capillaries  appear  at  times  to  be  concentrated  on  tne  forward  face  of 
the  large-scale  wave  just  before  the  sharp  crest.  In  this  case, 
perturbation  techniques  are  rather  difficult  to  apply,  because  the 
effect  of  surface  tension  is  not  uni fora  over  the  whole  disturbed 
surface.  Longue t-Higg ins  [1963}  introduced  surface  tension  as  a 
perturbation  effect  on  large-scale  waves  in  the  region  of  the  sharp 
crest  and  has  shown  the  reason  for  the  occurrence  of  this  phenomenon. 

4.2.2  Capillary  Generation  bv  Wind 

One  of  the  universal  features  of  wind-generated  wave  systems  is 
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the  Modulation  of  velocity  and  amplitude  of  small  capillary  waves  by 
large-scale  gravity  waves.  These  aodnlations  are  of  practical 
importance.  For  example,  the  iaaging  properties  of  aicrowave  radars 
need  for  reaote-senaing  of  ocean  waves  depend  on  the  modulation  of 
surface  seatterer  density,  aaplitude,  and  velocity.  Furthermore,  the 
short  gravity  and  capillary  waves  can  support  the  entire  stress  of  the 
wind  on  the  water  [Larson  and  Wright,  1975].  If  the  large-scale  waves 
aodulate  the  saall  capillary  wave  aaplitudes,  they  also  modulate  the 
stress.  The  component  of  this  modulation  stress  which  is  in  phase  with 
the  horizontal  component  of  orbital  speed  of  the  large-scale  wave  will 
work  on  the  large-scale  wave  and  cause  it  to  grow.  The  mechanism  of 
energy  transfer  between  short  and  long  waves,  discussed  by  Hasselman 
[1971]  also  depends  on  the  phase  and  amplitude  of  this  modulation. 

4.2.3  $ou T<?ay  oj  Module^ iQns 

At  least  three  sources  of  modulation  have  been  noted  so  far: 
modulation  of  the  orientation  of  small  capillary  waves  (called 
tilting),  modulation  of  the  amplitude  of  the  small  capillary  waves 
exemplified  by  the  straining  of  the  capillary  waves  by  the  horizontal 
component  of  the  orbital  velocity  of  the  large-scale  waves 
[Longue t-Higgins  and  Stewart  (1964),  Phillips  (1966)],  and  modulation 
due  to  wave-induced  airflow.  The  last  type  of  modulation,  wave-induced 
airflow,  is  a  principal  result  in  the  recent  work  of  Wright,  Plant  and 
Seller  [1980].  The  dependence  of  modulation  on  modulating  ocean  wave 
frequency  and,  partially,  upon  wind  speed  are  qualitatively  consistent 
with  wave-induced  airflow  as  a  modulator,  but  positive  identification 


of  the  modulation  sources  remains  to  be  made. 


4.3 


4.3.1  History  tad  Basis  _of  Equations 

The  current  scstter  theory  is  •  refinement  of  the  two-scale  model 
developed  in  the  1960's  [Bass  and  Fuks  (1968),  Fuks  (1969),  Valenzuela 
(1967,  1968)].  The  success  of  such  a  theory  in 

explaining  dependence  on  wind  speed,  wind  direction,  angle  of 

incidence  and  frequency  depends  almost  exclusively  upon  the  adequacy  of 
knowledge  of  the  sea  spectrum. 

The  derivation  of  the  average  backscattering  coefficient  for  the 
sea  surface  involves  two  major  steps:  (1)  compute  the  backscattering 
coefficient  for  the  capillary  waves  using  standard  perturbation 
techniques  [Valenzuela  (1967),  fright  (1966)]  and  (2)  account  for  the 
tilting  effect  of  the  large-scale  waves.  The  expression  for  the 
polarized  scattering  coefficients  from  step  (1)  is  well  known 
[Valenzuela  (1967),  Chan  and  Fung  (1973)]  and  is 

f)z  \*ff  \ZWLe,  f)  (4.13) 

where  for  horizontal  polarization  p“h  and  for  vertical  polarization 
p*v,  k  is  the  electromagnetic  wave  number,  Q  is  the  incidence  angle, 
is  the  aspect  angle  relative  to  the  upwind  direction,  and  is 

the  roughness  spectrum  of  the  sea  surface.  The  coefficients,  etpp,  for 
different  polarization  states  are  defined  as 

^hb  -  Cc>s(*)  (4.14) 


35 


(4.15) 


r 


fv  c»t(«)  ♦  (it-  /)  Tv' SCje/U* 


where  P  is  tbs  incidence  angle,  Qf  is  the  aspect  anglii  relative  to  tne 
upwind  direction,  k'  is  the  wave  amber  in  sea  water.  Rfc.By  are  Fresnel 
reflection  coefficients  for  horizontal  and  vertical  polarizations, 
respectively,  and  Ty-l+Ry.  The  depolarization  scattering  coefficient 
is  also  known  [Valenzuela,  1967]. 


C0 

nmdn 


(4.16) 


-fl* 


where ; 

2 

D-  ^  (fc-srH-nf)  ^  (4.i7) 

Y\  $*»0)Cc>S(f)  (4.18) 

+ 

—  *7)  ir  ^  S**r&)  S**(p)  (4.19) 


To  account  for  the  tilting  effect  of  large-scale  waves,  it  is 
necessary  to  identify  polarization  changes  due  to  tilting  and  then  to 
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iTirt|t  the  resulting  expression  for  the  sesttering  coefficient  of  the 

cspillary  neves  over  the  slope  distribution  of  the  large-scale  waves. 

The  concept  involved  in  considering  polarization  changes  can  be 

described  in  teras  of  two  sets  of  coordinates:  (1)  the  reference 

coordinates  in  which  the  scatter  problea  is  posed,  and  (2)  the  local 

coordinates  in  which  scattering  due  to  capillary  waves  is  coaputed. 
A  A 

Let  v,  h,  be  the  unit  vertical  and  horizontal  polarization  vectors  in 

A  A 

the  reference  fraae  and  v',  h'  be  the  corresponding  unit  vectors  in  the 
local  fraae.  Let  the  incident  field  be  ^  0.  Then  components  of  the 
local  incident  fields  are: 


f  ~  A ,  _ 

sh‘ 


(4.20) 


and 


<■  ‘  *  •  — 

ZV'  -  v.  S. 


(4.21) 


Upon  scattering  by  the  capillary  waves,  each  of  the  above  incident 
field  coaponents  is  depolarized  so  that  tne  scattered  fields  are 


*r  -  *  K+  f ,  £ 

Y  *y  »  ry  V 


(4.22) 


and 
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(4.23) 


| 

where  Spq  is  the  scattering  function  for  the  capillary  waves.  The 
scattered  fields  in  the  reference  frame  become 


(4.24) 


(h-h')  £h'+  (h.v') 


(4.25) 


It  is  clear  from  (4.13)  that  to  compute  the  scattering  coefficient  for 
horixontal  polarization  the  local  vertically-polarized  scattered  field 
is  needed  and  vice  versa.  Furthermore,  there  is  a  term  involving  the 
product  of  £ fo*  «o«l  2fv' 

Special  cases  where  S^'y',  ^v’h'  are  zero  and  tilt  angles  are 
small  have  been  treated  by  Wright  [1968],  and  cases  where  the  tilt 
angle  is  restricted  to  the  plane  perpendicular  to  the  plane  of 
incidence  have  been  studied  by  Valenzuela  [1968].  Note  that  Valenzuela 
did  not  average  over  the  tilt  angles  as  Wright  did.  Thus,  his  results 
show  the  effect  of  polarization  changes.  The  study  of  these  special 
cases  indicates  that  the  polarization  effects  due  to  tilting  are 
unimportant  for  vertical  polarization.  The  same  statement  can  be  made 
for  horizontal  polarization  when  the  incidence  angle  is  restricted  to 
less  than  about  70°.  Hence,  when  averaging  is  included,  the  scattering 
coefficient  for  polarization  scattering  is  approximately 
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(4.26) 


where  is  given  (4.1),  &  '  is  the  local  angle  of  incidence, 
Po(Zx',ZyO  is  tlie  slope  distribution  of  the  large-scale  waves  as 
▼ieved  at  an  incidence  angle, &,  and  is  defined  in  the  priae  coordinate 
whose  x'  axis  is  parallel  to  the  wind  direction.  Zx',Zy'are  partial 
derivatives  of  Z.  It  is  assuaed  that  the  plane  of  incidence  is  the  x-z 
plane  and  that  t&e  angle  between  tne  x  axis  and  the  x'  axis  is  f , 
(Figure  4.1)  so  that  an  upwind  observation  occurs  when^**0.  The 
relations  between  the  priaed  and  oaprined  coordinates  are 


Zy  Coi(f)  -  2*  u-2,) 

The  cross-polarized  scattering  coefficient,  including  polarization 
and  averaging  effects  for  saall  tilts,  is  approx iaat el y 

QO  OO 

<r  gu  J  J  Cf'Mf)  fa.  d £(z. 

where 

2 

9)  =  ?  !zy  (,ivv- ^Yssi  ff’j  wfaf)'*-™ 

and^  hv^**^  i*  *iT«a  by  (4.1). 

The  relation  between  tne  slope  and  density  function  P^(Zx>,Zyf) 
and  function  Pp(Zx',Zy>)  defined  by  Cox  and  Nunk  [19S4]  is 
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Diagram  of  the  relation  between 
(2,X,Y)  and  (Z,X',Y')  coordinates. 
[Chan  and  Fung,  1977] 


%  (?,-  >  ?y)  =  2y')  <4-31’ 

Hie  fora  of  P(Zx>,Zy>)  has  been  found  to  bo  a  Grom-Charl ier  tar lea 
(Cox  and  Monk  (1954),  Longue t-Higgins  (1963),  Pierson  ( 1975)1 . 
However ,  the  values  of  the  slope-variance  parameters  in  the  function 
require  further  studies.  The  slope  variances,  0'a  ,  •  given  by 

Cox  and  Munfc  [1954]  are  not  restricted  to  large-scale  waves.  It  is  not 
clear  at  the  present  tine  how  the  slope  variances  for  the  large-scale 
wave  alone  can  be  obtained  at  a  given  nicrowave  frequency.  It  has  been 
suggested  that  one  possible  way  to  coaipute  these  slope  variances  is  to 
integrate  the  low-frequency  portion  of  the  slope  spectrum  of  the  sea 
surface  [Chan  and  Fung,  1977 J.  If  we  assume  that  some  suitable  slope 
variance  parameters  can  be  found,  then  sll  the  quantities  in  the 
scattering  coefficient  have  been  discussed  except  the  sea  spectrum  for 
the  sea  spectrum  1(0  ,f)  to  be  discussed  in  the  following  section. 
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4.3.2 


The  high-frequency  portion  of  the  sea  spectrua  has  undergone 
intense  study  in  recent  years  [Pierson  (1973) ,  Mitsuyasu  and  Honda 
(1974),  Afarani  (1973)1.  The  general  fora  of  the  directional  sea 
spectrua  is 


S(K,<p)s  SLK)(a.*a,c0s(.z9))  (4.32) 

when  higher  order  terns  in  are  neglected  [Pierson,  1973] .  The 
quantity  determined  froa  measurement  is  S(K) .  One  fora  for  S(K) 
proposed  by  Pierson  [1975]  for  the  fully-developed  sea  is  as  follows: 

[The  large-K  regions  are  baaed  on  wave-tank  aeasureaents.  It  should  be 
noted  that  further  revision  of  S(I)  is  being  conteaplated  [Pierson, 
19781.] 
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Kg  oan  be  found  nuaerically  by  setting  84(^4)  equal  to  Ss(*4  )  . 

Other 
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paraaeters  are  as  follows: 


u*  »  friction  velocity  n*  >  u*B 
I.  -  (13.1769) 

q  -  log[Sg(t2)/S4(E3)]/log[K2^3] 

Pj  *  5.0  -  log  u* 

ZQ  -  0.684/a*  +4.28  x  10“5u*2  -  4.43  x  10-2 
u(u*)  “  u*/0.4  Ln(Z/Zg)  cm/ sec. 
a  »  4.05  x  10"3 
g  *  980  cm/ sec. 
u*m  *  12  cm/ sec. 

Note  that  this  sea  spectram  shows  that  the  capillary  waves  (i.e., 
^4(1))  grow  with  the  wind  and  that  the  higher  tne  £  valae  the  faster 
the  growth.  The  parameters  ag  and  ai  have  been  derived  by  Chan  and 
Fnng  [1977]  and  the  relation  between  1(0',  $)  in  (2)  and  (3)  and  S(£, 
has  been  found  to  be: 

*  4-  Z 

where  I  ■  ltsing’,  V  *  ?^Ct)  ^«T  are  tne  slope  variances  of 

the  sea  surface  along  the  crosswind  and  upwind  directions.  A  possible 

way  of  estimating  V  is  to  use  the  slope  variance  parameters  given  by 

Cox  and  Mnnh  [1954]  for  their  clean  sea  model.  Note  that  1(0  ,f)  is 

resrioted  to  that  portion  of  the  sea  spectrum  which  satisfies  the 

assumptions  in  the  perturbation  theory.  Thus,  the  value  of  £  has  a 

lower  bound  which  is  a  function  of  the  electromagnetic  wavelength. 

With  the  sea  spectrum  given  in  the  preceding  section,  the 

sea-scatter  model  given  by  (4.14)  can  be  evaluated.  Figure  4.2  shows  a 

.  o 

comparison  of  the  theoretical  0  versus  azimuth  angle  curves  with 
experimental  data  taken  on  two  different  AAFE  RAD SCAT  flights  [Jones, 
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Schroeder  and  Mitchell,  1977,  1978].  Figure  4.3  illustrates  the 
relation  between  scattering  coefficient  and  wind  speed  at  tnree  angles 
for  horizontal  polarization  and  13.9  GHz  frequency,  along  with 
experimental  points  obtained  with  the  AAFE  RADSCAT  instrument.  The 
lines  on  the  figures  are  theoretical  values  based  on  the  CT'  model  and 
spectrum  of  Section  4.3.2. 

Our  following  analysis  for  determination  of  theoretical  scattering 

coef f icient<f*  ,  is  based  on  having  instantaneous  slope.  One  can 
Htiiofy  , 

derive  instantaneous  local  angle  of  incidence,  @  ,  and  also  roughness 

spectrum  W(  &  ,  )  and  directly  calculate  the  instantaneous 


Figure  4.2 

Comparison  of  theoretical  results  with  a  sample  AAFE 
run  [Chan  and  Fung,  1977]. 
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oa  these  modulations  [Wright,  1977]  which  ney  also  induce  transfer  of 
substantial  energy  between  short  capillary  and  long  large-scale  waves. 

An  FM  radar,  such  as  that  used  by  K.O.,  is  a  self-contained 
two-scale  wave  probe.  The  measure  of  the  large-scale  waves  (several 
aeters  in  length)  is  the  instantaneous  wave  height,  which  is  obtained 
froa  tiae-delay  of  the  received  signal  (slant  range)  (Section  5.3).  At 
appropriately  chosen  viewing  angles  the  primary  scatterers  are 
capillary  waves,  a  few  centimeters  in  length.  Hence  the  amplitude 
aodnlation  of  the  received  signal  is  a  measure  of  the  modulation  of  the 
capillary  waves.  Two  inportant  sources  of  modulation  of  scattered 
microwave  power  have  been  noted:  modulation  of  the  orientation  of  the 
capillary  waves,  called  tilting,  and  modulation  of  the  amplitude  of  the 
capillary  waves,  ezeaplified  by  the  straining  of  the  short  capillary 

waves  by  the  horizontal  component  of  orbital  velocity  of  the  large 
waves. 

If  one  assuaes  that  the  capillary  wave  amplitudes  are  uniform  over 
the  large-scale  waves,  the  modulation  of  the  received  signal  is 
completely  governed  by  the  variation  of  the  local  angle  of  incidence 
(tilting  effect)  caused  by  slope  changes.  Note  that  tne  instantaneous 
large-scale  wave  slope  can  easily  be  derived  froa  the  instantaneous 
wave  height  time  profile,  provided  a  long-crested  assumption  is  made. 

If  one  knows  the  spectrum  of  the  capillary  waves  S(K)  and  the 
slope  of  the  large-scale  waves,  the  scattering  coefficient  can  also  be 
established  (Section  6.2).  Since  the  capillary  wave  amplitudes  are  not 
the  aaae  at  all  points  on  the  large-scale  waves,  due  to  modulation,  one 
of  the  major  goals  is  to  establish  their  distribution.  This  task  can 
be  done  by  coaparing  the  radar  signals  calculated  as  mentioned  above 


with  the  actual  radar  signal  observed.  The  difference  must  be  due  to  a 


change  in  the  amplitude  of  the  capillary  waves  at  tne  Bragg-resonant 
wavelength.  Then  the  variation  of  this  amplitude  difference  is 
cross-correlated  with  the  height  or  slope  of  the  large-scale  waves,  an 
estimate  of  the  phase  relations  of  the  peak  of  the  distribution  of 
'apparent'  capillary  wave  amplitude  (with  the  crest  of  the  wave  or 
point  of  SMzimnm  slope)  can  be  obtained  (Section  8.0). 

The  term  ’apparent’  is  used  because  interactions  between  the 
larger  wavelength  ocean  waves  and  the  capillary  wave  can  result  in 
stretching  or  compressing  a  persistent  train  of  capillary  waves,  so 
that  the  wavelength  associated  with  a  particular  pair  of  wave  crests 
may  be  different,  depending  upon  their  position  on  the  large-scale 
wave.  Since  the  Bragg  resonant  condition  for  the  radar  can  not  follow 

these  expansions  and  contractions  of  the  wavelengths  for  a  wave  train, 
the  signal  observed  might  actually  come  from  different  wavelets, 
depending  upon  their  position  on  the  large-scale  waves,  even  if  the 
local  angle  of  incidence  were  to  remain  the  same. 


4.4.2  Definition  of  the  Modulation  Index 

Note  that  the  scattering  coefficient,  Ol*.  ,  can  be  resolved  into 

total 

the  two  components:  contribution  of  backscattering  due  to  orientation 


aodulation  (tilting) 


<r° 


f//r 


,  and  contribution  of  backscattering  due  to 


amplitude  modulation,  (J*  .  . 


Hence 


CT'  "  —  ^  J  £7"  °  (4.41) 

total  ~  tlfr  n»md 


but 


(4.42) 
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(4.43) 


(4.44) 


(4.45) 


vh«r«  |i  is  called  the  'modulation  index.' 

Hie  concept  of  the  modulation  index  is  fairly  straightforward.  It 
is  simply  the  ratio  of  the  contribution  of  amplitude  modulation  to  the 
contribution  of  tilting  to  the  baekscatter.  The  phase  relations  of  the 
modulation  index  with  respect  to  the  crest  and  slope  are  discussed  in 
Section  8.0. 


5.0  EIPEBTMFNTAL  SET-OP.  MEASUREMENT  TECHNIQUES  AND  'EVENT*  SOLUTION 
5.1  Experimental  Set-Up  and  System  Description 

The  University  of  Kansas  system  was  an  FM  radar  operated  over  the 
band  from  8-18  GHz,  Ku-band.  The  block  diagram  of  the  system  is  shown 
in  Fignre  5.1.  The  platform  near  the  Dutch  coast  was  instrumented  with 
both  microwave  and  ocean-wave  equipment  for  measurement  of  waves, 
currents  and  tides  (Figures  5.2a  and  5.2b). 


5.2  Basis  of  Experiment 

The  microwave  antenna  was  mounted  just  below  the  helicopter  deck 
of  the  platform  (Figures  5.3a  and  5.3b).  It  illuminated  a  portion  of 
the  ocean  surface  that  was  small  in  comparison  with  the  wavelength  of 
any  modulating  (large-scale)  ocean  wave  of  interest.  The  short 
capillary  waves,  the  predominant  microwave  scatterers  are  advected 
about  by  the  large  waves.  The  backscatter  field  caused  by 
Bragg-resonant  phenomena  was  measured  at  the  receiver.  From  the 
frequency  difference  between  transmitted  and  received  signals,  the 
slant  range  between  the  radar  and  the  centroid  of  the  illuminated  area 
could  be  found  and  converted  to  wave  height.  From  the  signal  return 
and  the  instantaneous  area  of  the  observed  spot,  the  scattering 
coefficient  could  be  determined.  The  result  of  these  conversions  would 
be  a  sample  time  profile  of  wave  height  in  the  center  of  the 
scattering-measurement  area  which  is  comparable  with  that  obtained  by  a 
wave  gauge  at  the  tower  and  a  simultaneous  time  history  of  the 
scattering  coefficient.  Brief  remarks  on  radar  equations  and 
calculation  of  scattering  coefficients  and  range  follow.  Note  that  the 
simplification  of  the  radar  equation  when  the  radar  parameters  remain 
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west -side 


south-side 


1.  poles  with  oceanographic  sensors 
A.  place  of  fM/CW  measuring  radar  in  1977 
8.  ditto  in  1973 

C.  place  of  french  RAftSES  system 
0.  place  of  8  mm  navigation  radar 

L.  luneperg  lens,  protruding  from  lielicooter  deck  (calibration) 


Figure  5«2a 

Measurement  platform  Noordwijk 


50 


figure  5.3 

Illustration  of  antenna  mounted  at  tie  i 
of  tde  Noordwi ik  Tower 


tially  unchanged  caa  ba  represented  by 


Pr-  ftGtGr 
Lkn)3  Z* 


(5.1) 


where : 

Pr  ■  received  power 

Pt  “  transaitted  power 

Gg  »  traasait  antenna  gaia 

Gr  ■  receive  aateaaa  gaia 

'X  *  waveleagth 

©  ■  scatteriag  coefficient  (scattering  cross-section  per 
oait  area) 

Aill  *  illaaiaated  area 

R  *  slant  range  to  centroid  of  observation  spot 

The  retara  power  is  aeasared  by  ase  of  a  square-law  detector.  The 
power  aeasared  by  the  detector  aay  be  related  to  tae  returned  power  at 
the  receive-anteana  terminal  through  aa  unknown  constant  Kt,  which 
represents  the  effects  of  attenuation  and  conversion  losses  between 
antenna  and  detector.  Pdt»  the  detector  power,  is  therefore  related  to 
the  returned  power  by 


P  ft  GtGr  'f/i’i- "I 
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(5.2) 


:  'JS&JllE&SCi.'  .a. A*, sn '  AMMiditaM  a  n 


Iaaediately  before  end  after  recording  the  return  from  the  ocean, 
a  coaxial  delay  line  of  loss  L  is  switched  into  the  circuit  at  tne 
antenna  ports  to  replace  the  path  through  the  transnit  antenna  to  the 
surface  and  bach  through  the  receive  antenna. 

The  power  (through  the  delay  line)  after  square-law  detection  is 
given  by 

9  =  Kl  ft  L 

PIT  * 


Taking  the  ratio  of  equations  (5.2)  and  (5.3) 

.  f  y  (T  */4/lL  J 

%LT  ^  (4/7  )3 /f/ L  J 


To  conplete  calibration  of  the  system,  the  returns  fron  a  standard 
target  of  a  known  radar  cross-section  sre  measured.  The  Luneberg  lens 
of  the  Dutch  group  was  used  as  a  standard  target.  The  returned  power 
of  standard  radar  targets  is  given  by 


r  'X  & 

Ums'  c  1 -(tin)3  . 


(5.5) 


where: 


Kc  “  system  constant  k  at  time  of  calibration 
“  range  to  standard  radar  target 
(Tsrt  *  scattering  coefficient  of  standard  target 

The  power  detected  through  the  coaxial  delay  line  at  the  time  of 
calibration  using  the  standard  radar  target  is  given  by 

P  =  fcc  Pi  £  <*•«> 


The  ratio  of  equations  (5.5)  and  (S.6)  is  given  by: 


L 


(5.7) 


Combining  equation  (5.4)  with  (5.7)  and  solving  for 

%f*s 

Note  that  the  radar  cross-section  of  the  Luneberg  lens  was 
supplied  by  the  manufacturer,  but  a  recalibration  was  performed  at  the 
University  of  Kansas  [Kim,  1982]. 

The  illuminated  area  was  calculated  using  the  geometry  shown  in 
Figure  5.4.  The  projection  of  the  radar  beam  as  seen  on  the  flat 
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P 

'per 


surface  is  a  skewed  ellipse.  The  area  of  the  ellipse  was  calculated 


JL(Sf  )Ut  ) 

*■  »*ij  '  \  *ris) 


where  Max^s  it  the  major  axis  and  naX£S  is  the  minor  axis. 

Froa  the  geometry  of  Figure  5.4.  the  expressions  for  the  aajor  and 
ainor  axes  were  derived. 

M  -  Cos  * 

ent'$ 


where: 

Maz^s  «  aajor  axis 
aaX£S  *  ainor  axis 
R  *  range  to  center  of  target 
9  *  pointing  angle  of  antennas  off  vertical 
^  E  effective  gain-product  (for  receiving  and  transmitting 
antennas)  beaanridth  in  elevation  plane 
At  »  effective  gain-product  beaawidth  in  the  azimuthal  plane 

For  a  given  target  at  range  R.  the  time  that  returned  signals  are 
delayed  is 


c. 


(5.11) 
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An  illustration  of  the  frequency  relationship  between  tne  transmitted 
and  received  signal  is  shown  in  Figure  5.5.  It  can  be  shown  that 


'5.12/ 


so 

*(*/)& 


(5.13) 


where : 

c  -  speed  of  light 
Fjjr  »  intermediate  frequency 
Fjj  -  sweep  rate 
^  f  *  sweep  frequency  width 

In  this  system,  Fjp  was  fixed  st  50  kHz  and  Af  was  fixed  at  1  GHz. 
Fy  is  adjusted  by  a  tracking  loop  to  center  tne  target  return  at  50 
kHz. 

5 .3  Measurement  Techniques 

Measurements  were  made  at  angles  of  incidence,  Q  .  from  10c  to 
O 

80  (Figure  5.6),  with  both  vertical  and  horizontal  polarizations. 
Frequencies  for  the  measurements  were  at  1  GHz  intervals  from  9  to  17 
GHz.  Measurements  were  made  at  various  angles  relative  to  tne  wind 
direction  with  a  large  number  of  measurements  in  the  upwind  and 
downwind  direction,  a  smaller  number  of  measurements  in  the  cross-wind 
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direction  end  e  few  measurements  in  which  the  azimuth  angle  was  varied 


during  collection  of  a  single  set  of  data. 

Sampling  rate  was  1  Hz  (one  sample  per  second)  ao  tne  maximum 
frequency  one  expected  to  see  was  0.5  Hz.  For  both  convenience  and 
assurance  that  aliasing  would  be  prevented,  the  wave  height  time  series 
were  low-pass  filtered  to  a  cutoff  frequency  of  0.35  Hz.  This  should 
not  eliminate  significant  information  in  the  wave  height  time  series, 
since  one  expects  about  95  percent  of  signal  energy  to  be  distributed 
between  0.1  and  0.2  Hz  in  normal  climate  conditions  for  wind  speed  less 
than  20  m/sec.  Each  set  of  experiments  usually  took  10-15  minntes.  At 
the  start  of  a  set,  the  wind  speed,  polarization  and  azimnth  were 
recorded  and,  in  each  set,  the  instantaneons  slant  range  (actually  Fjj) , 
wave  gauge  record,  and  power  return  were  recorded  on  the  cassette  tape. 

Examination  of  the  reduced  data  showed  that  a  very  large  number  of 
successful  measurements  was  made,  but  that  many  of  the  measurements 
contained  errors  due  to  such  causes  as  lack  of  proper  tracking  of  the 
illuminated  area  by  the  range-tracker  circuitry  and  poor 
signal-to-noise  ratio. 

Although  measurements  were  made  at  frequencies  separated  by  only  1 
GHz,  many  of  the  data  sets  were  obtained  at  only  two  frequencies  —  10 
and  15  GHz.  This  was  done  to  permit  longer  sampling  intervals  than 
would  have  been  possible  if  both  angle  of  incidence  and  frequency  had 
been  changed  over  the  entire  range  of  parameters.  Because  of  the 
variation  of  the  Bragg-resonant  wavelength  with  angle  of  incidence,  a 
nearly  complete  set  of  the  possible  Bragg-resonant  conditions  was 
obtained  with  the  two-frequency  measurements. 


5.4  'Specular  Events'  Observation 


Interesting  phenomena  that  were  observed  in  the  time  history  of 
radar  return  records  were  designated  'specular  events'.  Study  of  the 
time  history  of  radar  returns  indicates  that  aost  of  the  time  the 
signal  fades  about  a  mean  level  associated  with  the  local  angle  of 
incidence.  Occasionally  certain  'events'  occur  in  which  the  signal  is 
auch  larger.  This  has  been  noted  particularly  at  angles  of  incidence 
near  grazing  [Kalmykov  and  Pustovoytenko.  1976J,  but  has  also  been 
observed  at  other  angles  of  incidence.  An  example  of  this  type  of 
behavior  is  shown  in  Fignre  5.7  that  contains  six  noncontiguous 
128-second  samples  of  the  time  series  of  scattering  coefficient. 
Notice  the  generally  small  fluctuation  about  the  varying  mean  level 
with  the  occasional  events,  where  the  signal  increases  greatly.  No 
corresponding  events  occur  for  decreasing  signal,  so  these  large 
signals  appear  to  be  caused  by  different  phenomena  than  the  normal 
fading.  These  phenomena  were  observed  by  Kwoh  and  Lake  [1981]  in  a 
wave  tank,  where  they  identified  two  classes  of  events.  In  one  class, 
both  polarizations  gave  the  same  return,  and  the  laser  sensor  indicated 
that  specular  reflection  was  possible.  In  the  other  class,  the 
vertically  polarized  returns  were  larger  than  the  horizontal  returns. 
These  were  attributed  to  edge  diffraction  from  sharp  edges  where  waves 
approach  the  breaking  point. 

The  capillary-wave  return  power  is  proportional  to  the  square  of 
the  number  of  scatterers  within  the  footprint,  because  their 
backscattered  fields  add  coherently,  whereas  the  mean  return  from  the 
events  is  proportional  to  the  average  number  of  events  per  footprint. 
Since  the  events  are  randomly  located,  their  backscattered  fields  add 
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SigmaO  (dB) 


noncoherent ly.  The  significant  differences  which  have  been  reported 
between  the  scattering  coefficient  of  the  ocean  aeasnred  by  radars  in 
aircraft  and  spacecraft  and  those  aeasnred  by  radars  on  towers  aay  be 
explainable  with  this  theory.  The  larger  footprints  Bean  a  larger 
capillary  contribution  relative  to  the  event  contribution.  Therefore, 
the  specular  and  other  (perhaps  edge  diffraction)  events  cause  a  auch 
larger  aodification  to  the  Bragg-scatter  signal  strength  for  the  low 
altitude  (where  the  footprint  is  saaller)  than  for  the  high  altitude 
(where  the  footprint  is  large).  The  derivation  of  this  .conclusion  is 
shown  in  Appendix  B. 
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6.0  THE  BASIS  OF  THE  CALCULATIONS 
6.1  Slope  Calculation 


The  height  of  the  large-scale  gravity  waves  for  a  finite  record 
length  can  be  written  as  follows  (Section  4.1.5): 


(6.1a) 


=■  2-  C*i  (jtCK  +”**+&)  <6- 


lb) 


where : 

aj  »  the  amplitude  of  the  i£k  component  of  the  wave 
U>i  *  angular  frequency  of  the  i£k  component 
li  -  i£*  wave  number 

The  water  at  Noordwijh  is  about  20  m  deep.  This  is  close  enough 
to  the  condition  for  use  of  a  deep-water  dispersion  relation  (Section 
4.1.4)  that  such  a  relation  can  be  used  with  little  error.  Here  we 
used  the  deep-water  form: 


Z 


where  g  is  the  gravitational  constant  9. SI  m/sec^. 

Another  simplifying  assumption  made  here  is  that  tne  wave  is 
sufficiently  long-crested  that  its  slope  may  be  considered  to  be  in  the 
direction  of  wave  travel  in  the  -z  direction.  Making  this  assumption, 
slope  s(x,t)  is  given  by: 


6(x,t) 


s  —  4*  /^'c  fa  ?* ,3) 
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6.2 


Scattering  Coefficient  Calculation 
Tire  Bragg-resonant  wavelength  is  a  function  of  the  instantaneous 
angle  of  incidence.  The  angle  of  incidence,  normally  assumed  to  be  the 
pointing  angle  for  a  radar,  is  modified  by  the  slope  of  the 

underlying,  large-scale  wave  (Figure  S.7).  Thus  we  may  write  for  the 
Bragg-resonant  condition: 


(6.7) 


where  k  is  the  radar  wave  number  and  ^  is  the  pointing  angle, 
measured  from  vertical  for  the  radar  antenna. 

Let  the  local  angle  of  incidence  be  V  ,  where 


-  * 

so  the  instantaneous  theoretical 
would  be  (Section  4.3.1): 


(6.8) 


scattering  coefficient  o' 9). 

fh  Cm*y 


Z<r'w(0',?) 

rp 


(6.9) 


where  for  horizontal  polarization 


C.sV') 


(6.10a) 


Rji  being  the  Fresnel  reflection  coefficient  for  horizontal 
polarization: 
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(6.10b) 


s  Cos  £_  - 

c®>  & 


For  vertical  polarization 


«  _  £,  Code  ) ->■(£'- if )  T* 


(6.11a) 


8^  and  Ty  being  the  Fresnel  reflection  and  transmission  coefficients. 


For  vertical  polarization: 


(6.11b) 

A-  C 6  SC-  ’)Oji>  +///Vf<J/-Jf  V 


fysr  1+ 


(6.11c) 


In  k  i*  the  radar  wave  number  in  air,  k*  is  the  radar  wave 

number  in  sea  water,  and  9  is  the  local  angle  of  incidence. 

is  the  normalized  anisotropic  sea  spectrum,  and  O',  is 
the  variance  of  the  small-scale,  capillary  waves  as  shown  by  Fung  and 
Chan  [1977] . 


<r, V( *;  f)a  $00 [uz (&&) 


(6.12) 


Z  z  z 

where  R  »  ^"cZ/^UT  i«  the  ratio  of  slope  variance  at  crosswind,  <%t 
,  to  slope  variance  at  upwind,  .  and  K  *  2ksin(i^/ 

Determination  of  S(I) .  the  spectrum  of  the  small  capillary  wave  is 
based  on  a  modification  of  Pierson  and  Stacy's  [19731  sea  spectrum 
(Section  4.3.2).  The  region  of  their  spectral  model  that  we  are 
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interested  in  is  the  cspillary  ref ion,  84(1),  where 


6(k)=o.43t  iw  . 

The  paraaeters  which  have  been  used  in  this  model  have  been  described 
in  Section  4.3.2.  Even  though  the  mode)  looks  very  cumbers one,  the 

g.  r 

instantaneons  sea  spectrum  S ( K ( &  (t))),  based  on  instantaneous  local 
angle  of  incidence  and  windspeed,  can  easily  be  determined. 

In  the  part  of  the  experiment  reported  here  tne  antenna  pointed  to 
upwind  (  $  ”0°),  downwind  (  ^*180°),  or  crosswind  (  ^**9U°) ,  so  the 
scattering  coefficient  would  be  considered  only  local-angle-of- 
incidence,  &  ,  dependent.  Hence 

It  is  worthwhile  mentioning  that  the  last  factor  in  the  numerator, 
[1+2(1-R/1+R) ] is  constant  for  a  specific,  given  windspeed  a,  so  one  can 
state 

cr '(<?’<»)  *l£lu orfsfotScie'lAv/ 

"  rr  '  ffl  /Stiff) 

—  ^  ( &lf),  L‘)/4(U)  (6.15b) 


where: 

A(u)  is  constant  at  a  given  windspeed 
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For  determination  of  instantaneous  scattering  coefficient, 

O'  ,  we  normalized  the  instantaneous  calculated  value  of 

Mi 

rr8(9'nu)  to  *  ,  where  &  is  mean  of  3  (&(&/**) 

is  mean  of  actual,  measured  scattering  coefficient  (Section 
5.2).  This  type  of  normalization,  here  called  theoretical  scattering 
coefficient  calibration,  not  only  takes  care  of  A(u) , [1+2(1-R/1+R) J  but 
also  compensates  for  any  other  scaling  factor  that  we  did  not  take  into 
account  (Section  5.4).  The  theoretical  scattering  coefficient  is 
therefore  defined  as  ,  with 


(6.16) 


6.3  Modulation  Index  Calculation 

As  discussed  in  Section  4.3,  if  no  amplitude  modulation  of  the 
capillaries  were  present,  the  only  source  of  modulation  in  backscatter 
from  the  sea  surface  is  tilting,  but  it  has  been  shown  very  clearly  by 
several  investigators  (Section  3.0)  that  the  amplitude  modulation  of 
the  capillary  waves  does  indeed  exist  and  has  very  important 
contributions  to  actual  scattering  coefficient.  Therefore  we  may  write 

cr'  =  CT  °(fU))  -ecT  *  {&’(*))  (617) 

ZoXof  Tift 
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where : 


<f  ’ 

/7rt*J 


is  instantaneous  scattering  coefficient 
contribution  due  to  tilting  effect, 
is  instantaneoas  scattering  coefficient 
contribution  due  to  capillary 
wave  modulation. 


Then 


»  / 

d'  ($ 

totaf 


'(*>)  =  o'  "(fft,)  (u  \ 

t,fr  ' 


(6.18) 


Note  that  f  B  (t)  )  is  the  same  as 

t,fr 

theoretical  scattering  coefficient  (Section  4.3), 


cr*  (&'(t) ) 

ffieo'y 


the 


<r‘  (9'tv)  =0"  (/+  ) 

t.W  rif.ry  0. .  (&,mj  J 


(6.19) 


Here  we  define  the  modulation  index 


m. 


o'  (t) 
y**  (t)=  — 

O'  °  (t) 
7ht»<y 


So,  from  (6.19) 


/*(t)  *  f 

rn..ry 


(6.20) 


ytu) 

amplitude 
backscatter 
waves . 


is  the  ratio  of  the  contribution  of  backscatter  due  to 
modulation  of  capillary  waves  to  the  contribution  of 
due  to  orientation  of  capillary  waves  over  tne  large-scale 
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6.4  Cross-Correlation  Calculation 

Cross-correlation  of  modulation  index,  with  wave  height 

t  h  ( 0  ;  t)  can  readily  show  the  phase  relationships  between  the 

peak  of  capillary  amplitude  and  the  crest  of  the  large-scale  wave 
(Section  8.0).  Analytically. 


J(t)s  1 7 1./*!# !><•>*+* )</*]/_  fc 


2  (6.21a) 


Numerically,  Section  8.0 

Zaj 


(6  21b) 

With  the  cross-correlation  of  modulation  ind /**({),  with  slope, 
s(o,t)  one  can  determine  the  phase  relationships  between  peak  of 
capillary  modulation  and  peak  of  the  slope.  Analytically 

Scd*  [ 

Numerically,  Section  8.0: 

2*  / 

/('»’)=  Ytf  ZT /"(n)  (6 ,22b) 


6.5  Modulation  Transfer  Function  Calculation 

Consider  the  backscattered  signal  to  a  microwave  radar  which 
illuminates  an  area  of  sea  surface,  at  least  one  dimension  of  which  is 
much  smaller  than  the  wavelength  of  the  dominant  surtace  wave.  The 
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backscattered  power  is  found  to  fluctuate  with  time  and  may  be 
represented  by  Fourier  series 

P(t)s  rf„  C, s (2 a* /£+&) 

o 

-  PoCci(%)  .  (6.23b) 

In  many  cases,  harmonics  of  the  ocean  wave  frequency  are  not  found 
in  the  return  power.  Furthermore,  it  is  found  that  the  amplitude  of 
the  fluctuation  in  the  power  increases  with  the  height  of  the  ocean 
wave  [Plant,  1980].  To  incorporate  these  observations  in.  ,  let 

‘P  —  P 

P,  =  P/Zl/t 

0.  f 

Pr>  s  o 

where  jp  is  mean  backscattered  power,  A  is  wave  amplitude  and  R  is  a 
coefficient  describing  the  power  modulations  with  these  definitions. 

f(t)  =  P[l+  /Z/JCoifafi/.t+Pj}  (6.2!, 

Let  the  surface  displacement  be 


(6.24a) 


(  6 . 2  4o  ) 


(6.24c) 


(6 .24d) 


2 


/ICos(l<*  -In/1.*) 


(6.26) 


Tie  wave  slope  is  given  by 

it  —  -  k (k/4  -  2n  /*t) 


We  nay  therefore  write  tie  magnitude  of  tie  wave  slope  as 


It  /  -  k/4=  ^*/c 

b*  1 


(6.28) 


where  C  is  the  phase  speed  of  the  ocean  wave. 


and 


CC»  - 


bt 


(6.29) 


is  the  magnitude  of  the  horizontal  component  of  orbital  velocity. 
Thus. 


p[t)&  P  [it Iw/Us  Cos(2fl/,'£+g>jj  (6. 


30) 


where 


l/wl-  l£l 
/< 


(6.31) 


If  we  write  this  in  complex  notation,  we  have 


c  2n/0 1 

Ptt)s  p[/rb?/(^)c  J 


(6.32) 


•  ; *’  '  /’s  *-v 
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where  a  is  now  a  complex  quantity  whose  phase  is  positive  if  it  leads 


the  wave  crest.  Since  a  aay  depend  on  ocean  wave  frequency,  it  is 
called  the  modulation  transfer  function.  It  is  a  dimensionless 
quantity  and  is  related  to  S.  which  has  dimensions  of  inverse  length. 


/W  = 


(6.33) 


The  modulation  transfer  function  may  be  evaluated  by 


cross-correlating 


with  either  wave  height  or  horizontal  component 


of  orbital  velocity.  In  our  case  we  evaluate  Ml'F  by  cross-correlating 
with  wave  height.  Correlating  , 


ftt)  4  h  W)  =•  J.  ?R/)*e 


(6.34) 


where  r  is  the  time  lag.  But, 


czp/>  'C' 


b(°}t)+  hi*!*)  -  -L/4  € 

c* 


(6.35) 


P(t)¥r  h  fait)  -  ? #  (hl**t)+h%,t)  ) 


Transforming  the  equation  to  Fourier  domain. 


PL  f)M*)  =  9  *tt)  |H(/;V 


(6.36) 


(6.37) 


so 


<<P  /#(/;/’> 

The  plots  of  R(J^  )j  MIF  and  f)  sad  ffifi  versos  frequency  are  shown 
ia  Figures  8.22  through  8.27  (Seotioa  8.0). 


6.6  Cohcreace  Function  Calculation 

The  last  quantity  that  we  have  calculated  is  coherence  function, 
defined  by 

z 


(6.39) 


C 7/; _  jztrtswt*]] 

The  coherence  function  has  a  maximum  of  one  at  any  given  frequency. 


Values  less  than  one  indicate  either:  (1)  noisy  signal,  (2)  m  and 
//(j^are  not  linearly  dependent,  or  (3)  P(/)  is  due  to  sources  in 
addition  to  Htf) .  The  discussion  about  coherence  function  is  in 
Section  8.0.  Figures  8.28  and  8.29  (Section  8.0  show  the  plot  of 
coherence  vs  frequency. 
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7.0  TRANSFORMATION  FROM  TEMPORAL  TO  SPATIAL  DOMAIN 
7 .1  Problem  Description 


Measurement  of  wave  height  was  made  with  both  the  wave  gauge  and 
the  FM  radar  capability  of  TRAMAS  (Section  S.3).  These  point 
measurements  of  time  history  are  interesting,  bnt  many  of  the  features 
of  the  waves  that  are  desired  are  spatial  rather  than  temporal,  i.e., 
distribution  of  capillary  wave  amplitude  over  large-scale  wave.  In  tne 
following  we  discuss  the  means  to  convert  tne  temporal  information  into 
spatial  information  on  the  assumption  that  waves  are  long-crested 
(Section  4.0),  and  that  observations  are  made  in  an  upwind  or  downwind 
direction.  If  we  know  the  spectrum  of  the  capillary  waves  S(K) ,  and 
the  slope  of  the  large-scale  waves,  as  a  function  of  distance,  the 
theoretical  scattering  coefficient  may  also  be  established  as  a 
function  of  distance.  Furthermore,  if  one  could  prove  that  the 
transformation  of  information  from  temporal  to  spatial  domain  is  a 
linear  process,  the  direct  conversion  of  the  scattering  coefficient  or 
any  other  quantity  of  interest  from  temporal  to  spatial  domain  would  be 
possible. 

Since  the  capillary  wave  amplitudes  are  not  the  same  at  all  points 
on  the  large-scale  wave  (Section  3.0),  one  of  the  major  goals  is  to 
establish  their  distribution.  To  establish  the  spatial  distribution  of 
capillary  wave  amplitude  over  the  large-scale  wave,  one  must  convert 
the  waveheight-time  profile  to  a  waveheight-distance  profile.  Because 
of  the  nonlinear  relationship  between  tne  wave  number  K  and  the  angular 
frequency  ,  this  is  not  as  trivial  a  task  as  one  might  suppose.  In 
the  following,  we  undertake  to  lay  the  basis  for  tne  calculation. 


r/'Vi a  -> aSMMfthtS BiSs tew i  wfc 
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Consider  s  wave  traveling  in  the  -x  direction.  It  nay  be 
represented  by  a  series  of  the  fora  (Section  4.1.5): 


k  C-x.t)  =  %  &<:  ZeS  t  Cpt  ) 


(7.1a) 


where  b  j  is  a  complex  amplitude. 


(7.1b) 


W  -  *  ^ 


(7.2a) 


la  deep  water  the  relation  between  K  and  <0  is  (Section  4.1.4) 


(7.2b) 


In  the  radar  or  wave-gauge  measurement,  we  observe  the  height 
h (0,t)  where  it  is  assumed  that  tne  measurement  is  made  at  x*0  (Section 
5.3).  To  express  this  in  terms  of  Fourier  series  for  a  measurement  of 
duration  T,  we  may  write 


9i  *  J  A  f©,o 


e  srff 


(7.3) 


where 


*  ■  ijMarisiean  ^ 1 


A  PC 

T 


(7.4) 


CJt  * 


Recovering  the  tine  variation  of  height  from  the  Fourier  series 
then  consists  merely  of  evaluating 


Ji  ( o,t ) 


(7.5) 


If  somehow  we  could  observe,  at  a  single  instant,  the  height  as  a 
function  of  distance  h(x,0)  we  would  obtain  tne  Fourier  coefficient  for 
this  function  by 


where 


: 

-7a 


(7.6) 


X 


(7.7) 


We  could  regain  the  function  of  height  by 


h  (-X,  o) 


£*Zl>r 


e 


(7.8a) 


where 


lc» 


(7.8b) 


also 


78 


hC't.O') 


(7.8c) 


Note,  however,  that  the  terms  in  (7.8)  are  uniformly  spaced  in  K, 
whereas  the  terms  in  (7.5)  are  uniformly  spaced  in  .  The  Fourier 
coefficients  are  different  even  tnough  the  values  and  bn  are  samples 
of  the  same  continuous  spectra.  It  is  this  difference  that  causes 
problems  in  using  the  coefficients  obtained  by  (7.3)  to  evaluate 
(7.8a). 

The  slope  might  be  obtained  easily  by  differentiating  the  series 
of  (7.1),  that  is 


3- 


AC*,*) 

•c 


SO^it  +  ICcX.) 

e 


(7.9) 


3  A  (<»>*)  "dhiX.,0) 

so  either  or  '9’X.  can  be  recovered  by 

multiplying  terms  in  the  appropriate  series  by  jK. 


7.3  Continuous  Form 

It  is  instructive  to  consider  what  happens  with  continuous 
frequency  and  C  functions  rather  than  Fourier  series,  particularly 
since  the  FFT  algorithms,  although  actually  involving  Fourier  series, 
are  addressed  as  if  they  were  for  continuous  functions. 

In  continuous  form  we  may  write  the  time  series  as 


-  00 


(7.10) 


or  the  spatial  series  as 
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(7.11) 


Ji  Oe,o) 


oc 


o) 


e 


where 

<*>  *  f*f 


so 


(7.12) 


The  two  integrals  (7.10)  and  (7.12)  obviously  are  not  tbe  same 
because  K  tad  uj  ere  related  by  uj  -  Note  that  H(u)  and 

H(  vTg)  are  generally  not  equal. 

One  can  also  find  the  slope  in  either  the  time  or  tne  frequency 
doaain  by  differentiating  the  appropriate  functions  obtaining 

fu >t 


3L  A  (o,t)  =•  f  e  fAji 

~drx.  J  jf  X1T 


(7.13) 


. I f  K  rt CJSj )  F* -to 


(7.14) 


7.4  Numerical  Evaluation 

The  difficulty  with  converting  the  formal  solution  into  a 
numerical  one  has  to  do  with  the  way  in  which  the  FFT  provides  the 
values  for  the  function  H(t*>).  Clearly  the  inverse  FFT  requires 
samples  of  H(  (VTg)  spaced  uniformly  in  K,  but  the  uniform  spacing  in  u> 
means  that  the  values  are  non-uniformly  spaced  in  K.  Thus,  one  cannot 
directly  feed  the  outputs  of  the  FFT  of  h(0,t)  into  an  inverse  FFT  in  K 
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to  got  h(x, 0) . 


If  the  function  were  smoothly  varying  in  both  amplitude  and  phase, 
one  coaid  ase  an  interpolation  routine  to  obtain  values  of  the  function 
at  uniformly  spaced  points  on  the  E  axis,  but  unfortunately  the  values 
of  the  Fourier  coefficients  fluctuate  widely.  So  the  interpolation  is 
not  possible. 

The  straightforward  way  to  evaluate  the  function  is  siaiply  to  sum 
the  series  of  (7.1a)  with  t=0  and  E  given  by  (7.2b).  This  can  be  done 
in  the  computer,  but  not  with  the  inverse  FFT  routine.  Hence,  the 
nonnniformly  spaced  trigonometric  functions  themselves  must  be 
evaluated,  which  is  time-consuming.  For  handling  masses  of  data,  this 
means  that  a  table  must  be  generated  of  the  appropriate  trigonometric 
values,  so  that  multiplication  can  be  substituted  for  the  repeated 
evaluation  of  numerous  trigonometric  functions  by  selecting  a  standard 
record  length,  such  as  128  or  2S6  samples. 

Before  considering  the  setting-up  of  such  a  table,  we  must 
consider  the  apropriate  ranges  of  X  and  T.  We  can  write  the  phase  for 
a  component  of  the  surface  as 


If  we  are  to  use  values  of  X  and  T  that  correspond  to  comparable 
maximum  phase  shifts,  we  have 


T-  2^4.  x. /j 


(7.16a) 


where  T  is  the  sample  length  in  time  domain  and  fB  is  the  maximum 
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frequency  need.  Since  we  start  with  a  tine  record  of  limited  duration, 
X  must  be  expressed  in  terms  of 


X= 


(7.16b) 


which  means  that  we  shonld  not  try  to  reproduce  h(z,0)  over  a  larger 
range  than  this.  That  is,  we  must  keep 

X  4i  <7-i7) 


Reasonable  record  lengths  to  us~  with  the  MARSEN  data  are  64,  128 

and  sometimes  256  seconds.  If  we  continue  to  limit  the  maximum 
frequency  used  to  .35  Hz,  this  means  that  x  should  be  less  than  225, 
570  or  1140  m,  respectively.  The  predominant  wavelength  is  of  the 
order  of  54.0  m,  so  a  64-second  record  would  only  have  5  cycles  of  this 
wave  and  oven  fewer  cycles  of  the  low-frequency  components,  so 
correlation  achieved  with  individual  records  would  be  very  noisy. 
Accordingly,  the  minimum  record  length  that  shonld  be  used  is  128 
seconds.  An  analytical  experiment  was  perlormed  to  determine  the 
validity  of  this  procedure. 


I 

! 

f 


8.0 


A  very  large  data  set  (about  115  cassette  tapes)  was  collected 
during  the  experiment.  The  raw  cassette-tape  data  were  converted  into 
files  in  the  University  of  Kansas  Honeywell  computer  system.  The  files 
were  edited  to  insert  headers  ant?  eliminate  bad  runs.  Software  was 
written  to  produce  scattering  coefficients,  radar  wave  heights  and 
other  quantities  of  interest.  Preliminary  calculations  have  been  made 
on  all  the  basic  quantities. 

Figures  8-.1  and  8.2  illustrate  some  processed  'raw  data'  records. 
Figure  8.1  shows  a  sample  set  of  versus  time,  and  Figure  8.2  shows 
the  measured  wave  height  both  from  slant-range  measurement  and  wave 
staff.  The  power  spectrum  for  the  wave  height  profile  has  been 
calculated  both  from  radar  and  wave-staff  measurements.  Figures  8.3 
and  8.4  illustrate  samples  of  these  spectra.  They  agree  very  well. 
Both  have  a  single  peak  at  about  0.17  Hz.  This  agreement  gave  some 
confidence  for  further  analysis.  The  sampling  rate  of  the  experiment 
was  1.0  Hz.  For  slope  determination,  wave  height  data  set  has  been 
low— pass  filtered  with  maximum  frequency  of  0.35  Hz.  This  reduces  the 
possibility  of  aliasing.  Figure  8.5  shows  a  sample  histogram  of 
calculated  slope  for  a  wind  of  12.5  m/sec,  obtained  from  a  filtered 
wave— height  measurement.  The  standard  deviation  is  abont  10  degrees. 
For  wind  speed  cases,  say,  between  8  and  10  m/sec,  one  would  not  expect 
such  a  high  value  of  standard  deviation.  Figure  8.6a  illustrates  a 
sample  set  of  measured  scattering  coefficients  with  an  average 
windspeed  8.61  m/sec.  Note  that  any  point  which  is  greater  tnan  6.0  dB 
above  the  mean  is  called  specular  ’event'  (Section  5.4).  Figure  8.6b 
illustrates  the  same  sample  set  of  scattering  coefficients,  replacing 
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average  of 


tie  specular  events  by  a  seven-point  triangularly— weighted 
tie  adjacent  samples.  Tie  expected  dependence  of  scattering 
coefficient  on  tie  windspeed  las  been  observed,  and  is  illustrated  in 
Figure  8.7.  Tlis  dependence  is  due  to  tie  fact  that  tie  capillary  wave 
amplitude  depends  on  windspeed  (Section  4.3.2)  and  tie  amplitude  of  tie 
capillary  spectrum  directly  governs  tie  backscattered  power  (Section 
4.3.1) . 

Two  samples  of  measured  scattering  coefficient  for  W-  and 
HH-polarization  as  a  function  of  incidence  are  shown  in  Figure  8.8. 
Tie  scattering  coefficient  is  inversely  proportional  to  tie  incidence 
angle.  If  tie  effect  of  the  large-scale  wave  slope  on  capillary-wave 
amplitude  could  be  neglected,  tie  Bragg-wave  number,  K,  would  be 
completely  governed  by  the  angle  of  incidence.  Figure  8.9  shows  a 
sample  preliminary  estimate  of  tie  measured  spectrum  in  the 
Mitsuyasu-Honda  range,  S^K),  based  on  average  incidence  angle  (the 
pointing  angle).  The  measured  values  are  compared  with  theoretical  K 
spectra  (lines  on  the  figure).  Note  tnat  tne  K-spectrum  decreases 
rapidly  with  increasing  Bragg-wave  number,  K.  Since  the  scattering 
coefficient  is  directly  governed  by  K-spectrum,  one  would  expect  the 
smaller  value  of  scattering  coefficient  for  larger  incidence  angle. 
Note  that  Bragg-wave  number  K  is  proportional  to  angle  of  incidence. 

One  of  the  major  goals  of  this  analysis  is  to  estimate  the 
modulation  of  the  capillary  wives  over  the  large-scale  waves,  both 
temporally  and  spatially.  If  the  capilary  wave  spectrum  S(K)  and  the 
instantaneous  slope  of  the  large-scale  wave  are  known,  one  can 
calculate  the  theoretical  instantaneous  scattering  coefficient  based  on 
uniform  capillary  distribution  (Section  4.3).  The  RMS  value  of 
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SigmaO  (dS)  SigmaO  (dB)  SigmaO  (dB) 


AD' A 129  690  INTERACTION  OF  THE  RADAR  WAVES  WITH  THE  CAPILLARY  WAVES  %/  \ 
ON  THE  OCEAN.. (U)  KANSAS  UNIV/CENTER  FOR  RESEARCH  INC 
LAWRENCE  REMOTE  SENSING  L-.  S  BARKESHLI  ET  AL .  MAY  83 
UNCLASSIFIED  CRINC/RSl-TR-4 19- 1  N00014-79'C-0533  F/G  17/9  NL 


Radar  Wave 


Avm.  Of  5,  64Pt.  Svbr-unm  &  5Pt,  Moving  A 


Figure  8.3 

Normalized  wave-height  spectrum  (radar); 
single  peak  at  0.17  Hz. 


Avm.  OF  5,  64Pt.  Subrvns  &  5Pt.  Moving  Ave. 
Figure  8.4 

Normalized  wave-height  spectrum  (wave  staff); 
single  peak  at  0.16  Hz. 
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Figure  8.5 

Sample  slope  distribution  obtained  from  instantaneous  slopes 
derived  from  radar  wave-height  measurement.  Standard 
deviation  10®.  Wind  speed  12.5  m/s.  Scale  on  Figure:  Number 
of  observations  (upper  list  of  numbers). 

Slope  in  degrees  (lower  list  of  numbers). 


SigmaO  (dB) 


Time  Profile  of  SigmaO  in  dB 

U=*8. 69(m/sec),Theta*=40(Deg.)^zim.  «=0. 0( Oeg.),Freq.  *»  1 4.5( GHz),  W  Pol. 


SigmaO  (dB)  SigmoO  (dB)  SigmaO  (dB) 


352  384  416 

Time  (Sec.) 


544  576  608  640  672  704 

Time  (Sec.) 


Figure  8.6b 

Sample  of  a  (dB)  replacing  the  specular  events  by  7  point 
average  of  an  adjacent  sample.  Each  128-sec  record  is 
from  a  different  time. 
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SigmaO  (  dB  ) 


' 


SigmaO  vs.  log  windspeed  with  least  squares  fit. 
W  polarzation.  Inc.  angle  =*  40  deg.  Downwind. 


Figure  8.7 

Sample  preliminary  estimate  of  wind  speed 
response  of  o°  from  3  days  of  data 
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Incident  angle  (  deg.  ) 


Figure  8.8 

A  sample  preliminary  estimate  of  o°  vs  angle  of  incidence 
based  on  10  GHz  calibration  used  at  15  GHz. 


K-Spectrum  for 


K 


Figure  8.9 

Sample  preliminary  estimate  of  measured  Mi tsuyasu-Honda 
range  of  wave  spectrum.  Some  of  scatter  may  result  from 
use  of  10  GHz  reference  calibration 
for  both  10  GHz  and  15  GHz  data. 


Modulation  index  |i  or^| f*  (Section  4.4.2)  can  show  the  nodule tion 
effect  and  tilting  effect  on  total  aoattering  coefficient.  Seeing  the 
cross-correlation  between  the  ji(t)  and  slope  one  can  estimate  the  phase 
relation  of  the  peak  of  the  modulation  with  respect  to  peak  of  the 
slope.  The  phase  relation  of  the  peak  of  the  aodulation  with  respect 
to  crest  of  the  ware  can  be  found  by  the  same  means. 

The  ensemble-average  RMS  value  of  aodulation  index  has  been 
calculated  both  for  HH—  and  W— polarization.  The  average  RMS 
modulation  index  is  0.876  with  standard  deviation  of  0.34  for  HH-case 
and  is  0.673  with  standard  deviation  of  0.044  for  W-case. 

Figures  8.10a  and  8.11  illustrate  the  autocorrelation  of  wave 
height  and  slope.  The  wave  height  autocorrelations  show  that  the  wave 
height  is  a  quasi-per iodic  type  of  function  with  a  period  between  3.S 
to  6.3  seconds,  which  corresponds  to  a  dominant  wave  with  frequency  of 
about  0.13  -  0.17  Hx.  The  slope  autocorrelation  function  shows  that 
the  slope  is  the  same  type  of  the  quasi-periodic  function,  but  with 
shorter  periods,  about  3.3  seconds,  which  corresponds  tc  a  dominant 
frequency  of  0.29  Hz.  The  higher  dominant  frequency  in  slope  is 
expected  because  the  slope  spectrum  components  are  calculated  by 
multiplication  of  wave-height  spectrum  components  with  the  square  of 
corresponding  angular  frequency,^!,  divided  by  g,  the  gravitational 
acceleration.  This  would  definitely  shift  the  spectrum  peak  to  the 
higher  frequency.  Figure  8.10b  shows  the  cross-correlation  of  slope 
and  wave  height  and  cross-correlation  of  theoretical  scattering 
coefficient  and  slope.  As  is  shown,  the  slope  leads  the  wave  height  by 
about  1.0  second.  The  cross-correlation  of  theoretical  scattering 
coefficient  and  slope  has  a  maximum  at  zero  lag,  as  it  should. 
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MJTOCOOL  (MHTI 


-8.0-7.0-6.0-5.0-4.0-3.0-2.0-1.0  0.0  1.0  2.0  3.0  4.0  S.O  6.0  7.0  8.0 


TIME  Lag.  (Sec.) 


Figure  8.10a 

Temporal  autocorrelation  of  wave-height.  Dominant  period  is  5.8  sec. 
Temporal  autocorrelation  of  slope.  Dominant  period  is  3.5  sec. 


Ave.  of  10  Runs  Sept.21st/70  (SPCLR  OUT  6.  dB  GT  MEAN)  -  (FILT  .35  Hz) 


U-S . 61  («/Sec . ] .  Theta-40 .  (Deg . ) .  AZIM-0 . 0  (Deg . ) .  Fr eq-14 . 5  (GHz) .  Pol . -HH 


TIME  Lag.  (Sec.) 


TIME  Lag .  (Sec . ) 


Figure  8. 10b 

Cross-correlation  of  slope  and  wave-height.  Slope  leads  the 
wave-height  by  about  1.0  sec.  Cross-correlation  of  theoretical 
0°  and  slope.  Maximum  peak  at  zero  lag. 


Ave  of  5  Runs  Sept.  21st/79 


-  (FILT  .35  Hz) 


-8.0-7.0-6.0-5.0-4.0  -3.0  -2.0  -1.0  0.0  1.0  2.0  3.0  4.0  5.0  B.O  7.0  B.O 


TIME  Lag  .  (Sec  .) 


TIME  Lag .  (Sec . ) 


Figure  8.11 

Autocorrelation  of  wave-height.  Dominant  period  is  about  5-8  sec. 
Autocorrelation  of  slope.  Dominant  period  is  about  3-5  sec. 
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Figures  8.12  and  8.13  illustrate  the  cross-correlations  of  the 
seatteriag  coefficient  and  the  wave  height  and  also  tne  scattering 
coefficient  and  the  slope,  for  both  HH-  and  W-polarization  cases.  The 
values  of  cross-correlation  are  not  large,  since  the  graphs  are  the 
result  of  averaging  ten  runs  for  HH-cases  and  five  runs  for  W-cases. 
Even  though  auch  larger  correlations  have  been  observed  for  single 
cases  at  different  lags,  the  aost  significant  result  is  that  consistent 
values  were  observed  for  the  two  polarizations,  based  on  observations 
at  different  tines.  The  scattering  coefficient  leads  the  slope  by 
about  0.73  seconds  or  about  1/8  of  the  doainant  large-scale  wave  period 
(4 5°)  or  1/4.7  of  the  doainant  slope  period,  77°.  The  scattering 
coefficient  leads  the  wave  height  even  wore,  about  1.7  seconds  or  about 
104°  ahead  of  the  crest  of  the  doainant  wave.  Another  peak,  saaller 
than  the  doainant  one.  can  be  observed  for  both  cross-correlation 
functions.  The  phase  difference  between  tne  two  peaks  is  about  a  full 
period  of  the  doainant  slope  quasi-  periodic  coaponent,  3.3  seconds. 
This  obviously  shows  that  the  scattering  coefficient  is  doainated  by 
the  slope  effect,  which  one  would  ezpect  froa  the  theory  (Section 
4.3.1).  Note  that  the  cross-correlation  function  of  theoretical 
scattering  coefficient  and  slope  peaks  to  aaziaua  at  zero  lag  (Figure 
8.10b),  since  aaziaua  theoretical  scattering  coefficient  occurs  at 
aaziaua  slope.  Ve  believe  the  'unezpected'  shift  of  the 
cross-correlation  of  scattering  coefficient  and  slope,  1/8  of  the 
doainant  large-scale  period,  is  due  to  the  fact  that  the  capillary 
waves  are  not  uniforaly  distributed  over  the  large-scale  wave. 

The  aodulation  indez  |i(t)  is  called  ^  (J-'  in  the  figures.  For 
both  HH  and  W  cases.  Figures  8.14  and  8.13  show  that  A^leads  the 
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CRSCOflL  (SIGOCSLOPEI  CHSCOHL  (SIGOKHHT) 


Ave.  of  10  Runs  Sept.21st/79  (SPCLR  OUT  6.  dB  GT  MEAN)  -  (FILT  .35  Hz) 


U-8.81  (BI/S^C.J.  Thet«-40.  (Oeg.i,  AZIM-0.0  Frtq-14 .  S  (GHz) .  Pol.-HH 


TIME  Lag .  (Sec . ) 


0. 

0. 

0. 

-0. 

-0. 

-8.0-7.0-8.0-5.0-4.0  -3.0  -2.0  -1.0  0.0  1.0  2.0  3.0  4.0  5.0  8.0  7.0  8.0 

TIME  Lag .  (Sec . ) 

Figure  8. 12 

Cross-correlation  of  0°  and  wave-height.  0°  leads  the 
wave-height  by  about  1.7  sec  (104®).  Cross-correlation  of 
o°  and  slope.  o°  leads  the  slope  by  about  0.75  sec  (90°). 
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Ave.  of  S  Runs  Sept.21st/79  (SPCLR  CUT  6.  dB  GT  MEAN)  -  (FILT  .35  Hz) 


U-e .  89  ta/See . ) .  Tbete-40 .  (Deg . ) .  AZIM-0 . 0  (Deg . ) .  Freq-14 . 5  (GHz] .  Pol  ,-W 


TIME  Lag.  (Sec.) 


TIME  Lag .  (Sec  . ) 


Figure  8. 1 3 

Cross-correlation  of  o°  and  wave-height.  0°  leads  the 
wave-height  by  about  1.7  sec  (104®).  Cross-correlation  of 
o°  and  slope.  0°  leads  the  slope  by  about  0.75  sec  (90°). 


jlopi  by  «boot  1.6  seconds  or  about  1/3.6  of  tbo  do* inant  period  of  tba 
large-scale  van  (about  1/2.19  of  tba  dominant  period  of  tba  alope 
(about  163°)).  Tba  cross-  correlation  of  Aff"  and  ware  height,  Pigurea 
8.14  aud  8.15,  above  that  ^5"  laada  tba  wave  height  by  about  2.5 
aeeoade  vbieb  ia  about  1/2.35  of  doainant  period  of  tba  large-scale 
wave  or  about  153  dagreaa  ahead  of  tba  creat  of  tba  wave.  Note  tuat 
the  croaa-corralation  of  tba  and  wave  height  for  HH  ia  stronger 
than  for  W,  but  tba  locations  of  tba  peaks  are  reasonably  consistent. 

Since  tba  spatial  cross-correlation  functions  should  be  studied 
for  an  estiaation  of  the  spatial  distribution  of  the  capillary  waves 
over  tba  large-scale  waves,  a  set  of  spatial  profiles  has  been 
determined  using  the  aetbod  described  in  Chapter  7.  Figure  8.1 6  shows 
the  profiles  of  the  wave  height  for  a  saaple  record  in  the  teaporal  and 
spatial  doaains  and  Figure  8.17  shows  the  coaparable  slope  profiles. 
These  spatial  slope  profiles  were  need  to  compute  a  theoretical 
scattering  coefficient  profile  ia  the  spatial  doaain  as  shown  in  Figure 
8.18.  The  relationship  between  tne  scattering  coefficient  and  the  wave 
height  is  very  nonlinear.  Consequently,  the  assumptions  of  Chapter  7 
regarding  transformation  froa  teaporal  to  spatial  doaain  probably  are 
not  sufficiently  good  to  justify  a  direct  transforaation  of  the 
scattering  coefficient  to  the  spatial  doaain.  Accordingly,  a 
transforaation  has  been  aade  of  the  values  of  A9*  .  For  relatively 
saall  values  of  69"  the  relation  should  be  sufficiently  close  to 
linear  that  the  transforaation  has  a  reasonable  chance  of  success.  A 
saaple  is  shown  in  Figure  8.19.  For  larger  values  of69"  the  values 
presented  in  the  spatial  doaain  probably  are  not  very  accurate.  A 
simple  set  of  spatial  autocorrelations  is  presented  in  Figure  8.20  for 


CRSCOFS-  (Delta  S16CSL0PE)  CRSCORL  (Delta  SI6CWHTI 


Ave.  o£  10  Runs  Sept.21st/79  (SPCLR  CUT  6.  dB  GT  MEAN)  -  (FILT  .35  Hz) 


U-A.81  (m/See . ) .  Theta-40 .  (Deg .  J .  AZI.M-0 . 0  (Deg . ) .  Freq-i 4 . 5  (GHz] .  Pol . -HH 


TIME  Lag .  (Sec . ) 


TIME  Lag .  (Sec . ) 


Figure  8. 14 

Cross-correlation  of  A o  and  wave-height.  Ao  leads  the  wave-height 
by  about  2.5  sec  (153°).  Cross-correlation  of  Ao  and  slope. 

Ao  leads  the  slope  by  about  1.6  sec  (165°). 
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CRSCOflL  (Delta  SIGCSLOPE)  CRSCORL  (Delta  S1GCWHT) 


Ave.  of  5  Runs  Sept.  21st/79  (SPCLR  OUT  6.  dB  GT  MEAN)  -  (JILT  .35  Hz) 


U-6 . 88  (a/Sec . ) .  Theta-40 .  (Deg . ) .  AZIM-0 . 0  (Deg . ) .  Freq-14 . 5  (GHz) .  Pol .  -W 


TIME  Lag .  (Sec . ) 


•TIME  Lag.  (Sec.) 


Figure  8.15 

Cross-correlation  of  Aa  and  wave-height.  Aa  leads  the  wave-height 
by  about  2.5  sec  (153°).  Cross-correlation  of  Aa  and  slope. 

Aa  leads  the  slope  by  about  1.6  sec  (165°). 
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ware  height  and  slope.  It  seeas  the  doainant  spetiel  wave  length  is 
aboat  70  a  and  the  doainant  slope  ware  length  is  about  42  a.  Figure 
8.21  shows  the  spatial  eross-eorrelationa  between  AO*  and  ware  height 
and  slope  for  the  saae  data  set. 

Spatial  analysis  has  just  started  and  requires  considerably  aore 
study  before  conclusions  can  be  drawn.  Hence,  the  abore  presentations 
of  results  aust  be  considered  as  rery  preliminary . 

The  final  quantity  of  interest  is  the  modulation  transfer 
function,  MTF  (Section  6.5).  Figures  8.22  and  8.23  illustrate  the  MIF 
for  HH  and  W  polarization  cases.  Both  tne  Alpers  fora  B(f)  and  the 
NHL  fora  a(f)  are  shown.  Below  0.35  Hz  the  MIFs  hare  aore  or  less  tne 
saae  characteristics.  Rtf)  has  a  aininua  at  about  0.1  Hz,  and  a(f) 
decreases  a»notonieally.  More  saoo thing  and  aore  sample  areraging 
needs  to  be  done  for  aore  stable,  saooth  characteristics.  R(f)  abore 
0.35  Hz  has  aore  noiselike  characteristea  because  tne  ware-height  power 
speetrua,  which  is  in  the  denominator  of  the  MIF  expression,  would  not 
hare  proper  inforaation  beyoad  0.35  Hz  (Section  5.3).  For  both  W  and 
HH  polarization  cases  shown  in  Figures  8.24  and  8.25,  the  power  spectra 
hare  siailar  characteristics,  except  for  soae  peaks  at  0.22  Hz  and  0.35 
Hz  for  W  cases.  Note  that  for  W  only  5  saaple  records  were  areraged 
whereas  for  HH  10  saaples  hare  been  areraged.  Hence,  one  would  expect 
aore  saooth  and  stable  characteristics  in  power  spectra  and  MIF  for  HH 
case.  The  ware  height  power  spectra  of  both  cases  hare  been  plotted  in 
Figures  8.26  and  8.27.  Both  hare  a  single  peak  around  0.15  Hz. 

The  arerage  coherence  for  function  for  radar  power  return  spectra 
and  ware  height  spectra  hare  been  found  both  for  HH  and  W  polarization 
and  plotted  in  Figures  8.28  and  8.29.  The  arerage  ralue  of  coherence 


105 


WHT  (Meters)  WHT  (Metei$ 


Sample  Of  Wave  Height  Profile  In  TemporaldcSpatial  Domains 
U=*8. 61  (m/sec), Theta^40(Deg.)rAzim.—0.0(Deg.),Freq.—1 4.5(GHz),HH  Pol. 
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Figure  8.16 

Sample  of  wave-height  profiles  in  temporal  and  spatial  domains 


Sample  Of  Slope  Profile  In  Temporal&Spatial  Domains 
Uz*8. 6 1  (m/i sec),  Theta*‘40(Deg.)*Azim.  =*0. 0(Deg.),Freq. *  1 4.5(GHz),HH  Pol. 


Figure  8.17 

Sample  of  slope  profile  in  temporal  and  spatial  domains 


SigmaO  Theory  (dB)  SigmoO  Theory  (dB) 


Sample  Of  SigmaO  Theory  Profile  In  TemporaldcSpatial  Domains 
U=3.61  (m/sec),  Theta= 40 ( Deg.  ),Azim. =0.0(Deg.),f req. =14.5( GHz),HH  Pol. 


Figure  8. 1 8 

Sample  of  a  theory  profile  in  temporal  and  spatial  domains 
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Sample  Of  Delta  Sigma  Profile  In  Temporal&Spatial  Domains 
U-8. 61  (m/sec), Theta=40(Deg.),Azim.=0.0(Deg.),Freq.  =  l 4.5(GHz),HH  Pol. 


AUTOCOW.  (SLOPE)  AUTOCOPL  (NHT) 


Ave.  of  10  Runs  Sept.  2lst/79  -  (FILT  .35  Hz) 


U-e.Bl  («/S*c .  1 .  Theta-40 .  (Deg.) .  AZIM-0.0  (Deg.) .  Freq-14.5  (GHz) .  Pol.-HH 


Figure  8.20 

Spatial  autocorrelation  of  wave-height.  The  dominant 
period  is  about  72  meters. 

Spatial  autocorrelation  of  slope.  The  dominant 
period  is  about  42  meters. 
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Figure  8.21 

Cross-correlation  between  Aa  and  wave  height  and  slope. 
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Sept.  21st,  79,  Ave.  Of  10,  128pts.  Subruns  S  7pt.  Moving  Ave. 
11=9.61  (m/sec),  Theta=40  (Deg.) ,  Azim.=0.0  (Deg.j ,  Freq=14.5  (GHz) ,  HH  Po 

Spclr  Out  6.  dB  GT.  Mean 
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Figure  8.22 

Normalized  modulation  transfer  functions  R(f)  (dash  line) 
and  m(f )  (solid  line). 


Sept.  21st,  79,  Ave.  Of  5,  128pts.  Subruns  S  7pt.  Moving  Ave. 
U=B.69(m/sec),  Theta=40  (Deg.) ,  Azim.=0.0  (Deg.) ,  Freq=14.5  (GHz) ,  vv  Pol. 

Spclr  Out  6.  dB  GT.  Mean 
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Figure  8.23 

Normalized  modulation  transfer  functions  R ( f )  (dash  line)  and 
m(f)  (solid  line).  Note  it  is  not  as  smooth  and 
stable  as  the  HH  MTF. 


Sept.  21st,  79,  Ave.  Of  10,  12Bpts.  Subruns  S  7pt.  Moving  Ave. 
U=8.61  (m/sec),  Theta-40  (Deg.),  Azim.sO.  0  (Deg.),  Freq=14.5(GHz),  HH  Pol 

Spclr  Out  6.  dd  GT.  Mean 
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Figure  8.24 

Normalized  power  return  power  spectra. 


Sept.  21st,  79,  Ave.  Of  5,  128pts.  Subruns  8  7pt.  Moving  Ave. 
U*8. 69 (m/sec) ,  Theta*40(DegJ,  Aziw.=0.0(DegJ.  Freq=14.5(GHz),  vv  Poi 

Spclr  Out  6.  dB  GT.  Mean 


Freq.  (Hz) 


Figure  8.25 

Normalized  radar  power  return  spectra. 


Sept.  21st,  79,  Ave.  Of  10,  128pts.  Subruns  S  7pt.  Moving  Ave. 
U*8.61  (m/sec),  Thetas40(Deg.) ,  Azim.=0.0(Deg.) ,  Freq=14.5(GHz).HH  Pol 


Freq.  (Hz) 


Figure  8.26 

Normalized  wave-height  power  spectra; 
single  peak  about  0.15  Hz. 
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Sept.  21st,  79,  Ave.  Of  5,  12Bpts.  Subruns  S  7pt.  Moving  Ave. 
U*8. 69 (to/sec) ,  Tbeta^OfDegJ.Aziio.^O.OfDegJ.Freq^U.SfGHz).  vv  Pol 


Freq.  (Hz) 


Figure  8.27 

Normalized  wave-height  power  spectra 
single  peak  about  0.14  Hz. 


Coherence 


over  the  frequency  for  HE  is  0.321  with  standard  deviation  of  0.013  and 
for  W  is  0.72  with  standard  dewiaton  0.02.  The  coherence  functions 
show  nore  dependence  between  radar  power  return  spectra  and  ware  height 
for  W  than  for  HH.  Further  study  is  needed  to  explain  the  difference. 
Plant  [1980]  stated  that  a  useful,  though  arbitrary,  convention  in 
analysing  the  HIF  is  to  use  only  data  for  which  the  coherence  is 
greater  than  0.3.  One  result  of  a  low  value  of  coherence  function  is 
that  the  corresponding  values  of  MTF  have  large  variances. 


.4  vc.  Of  10  Ul in  Scpi.  PI  si/ 1979 

'J-  8.81  (m..  /s),Thcia--40.0(Dcg.),A:am ■  0.0(Dcg.),rrcq--1  4 .5(Ghz),Po!r-H! ! 


Freq.(Hz.) 


Figure  8.28 

Coherence  function  for  horizontal  polarization 
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Ave.  Of  5  Run  Sept .  21st/1979 

U=8.69(m/s),Theta=40.0(Deg.),Azim=0.0(Deg.),Freq=1 4.5(Ghz),Pol=W 


Figure  8.29 

Coherence  function  for  vertical  polarization 
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Methods  have  been  developed  for  analysis  of  the  University  of 
Kansas  aeasureaents  of  radar  cross-section  and  wave  height  as  a 
function  of  tiae  that  were  aade  as  part  of  Project  MARSEN  from  the 
Noordwijk  pi a tf ora  in  the  North  Sea  off  the  Dutch  coast  during 
Septeaber  through  Noveaber  1979.  The  radar  spectroaeter  was  unique  in 
two  ways: 

1.  It  could  operate  over  a  band  of  frequencies  froa  8  to  18  GHz. 
whereas  aost  previous  aeasureaents  have  been  made  at  single 
frequencies  or  groups  of  widely  separated  spot  frequencies. 

2.  The  range-tracker  in  the  systea  pemitted  determination  of  the 
instantaneous  wave  height  at  the  centroid  of  the  radar  footprint, 
whereas  previous  aeasureaents  have  depended  upon  wave  gauges 
located  soae  distance  froa  the  radar  footprint  or  upon  integration 
of  the  velocity  obtained  by  Doppler-frequency  measurement  to 
obtain  the  instantaneous  wave  height. 

The  first  capability  permits  observations  over  a  range  of 
Bragg-resonant  ocean  wave  numbers  (K)  not  previously  studied.  The 
second  permits  development  of  cross-correlations  between  radar 
cross-section  and  wave  height  not  previously  possible. 

The  overall  experiaent  provides  information  that  permits 
determination  of  the  apparent  wave  spectrua  S(K)  in  the  capillary 
region,  as  well  as  eapirical  observations  of  the  radar  scattering 
coefficient  over  a  wide  range  of  controlled  conditions.  These  studies 
are  the  subject  of  other  reports  in  this  series.  This  report 
concentrates  on  the  aethodology  and  preliainary  results  of  the  study  of 


> 

I, 


! 
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modulation  of  the  radar  cro»*-*ection,  and  consequently  of  the 
Bragg-resonant  capillary-wave  components,  associated  with  the  strncture 
of  the  underlying  larger  waves. 

Because  of  the  large  number  of  parameters  that  could  be  varied 
during  the  experiment  (frequency,  polarization,  angle  of  incidence, 
look  direction  relative  to  the  wind),  the  individual  observation  runs 
were  relatively  short.  This  means  that  the  statistics  of  this 
experiment  are  somewhat  less  stable  than  those  for  experiments  where 
fewer  parsmeters  had  to  be  varied.  Typical  data  runs  varied  from  3  to 
7  minutes.  However,  on  many  days  observations  with  the  same  radar 
parameters  were  repeated  often  enough  to  improve  the  stability  of  the 
statistical  results.  Because  of  the  relatively  short  duration  of  the 
individual  runs,  sisplif ication  of  data  processing  (involving  numerous 
FFT  operations)  was  achieved  by  standardizing  on  128-second  sequences. 
For  the  case  analyzed  as  an  example  in  this  report  (21  September  1979, 
40°  angle  of  incidence)  this  meant  that  averages  could  be  obtained  for 
3  vertically  polarized  (W)  and  10  horizontally  polarized  (HH) 
128-second  samples. 

Algorithms  were  developed  to  compute  the  modulation  transfer 
functions  (MTFs)  used  by  other  investigators,  so  that  comparisons  are 
facilitated.  Because  the  instantaneous  wave  heights  were  measured,  it 
is  possible  to  calculate  the  instantaneous  slopes  for  the  underlying 
larger  waves,  on  the  assumption  that  these  waves  are  long-crested.  One 
can  then  calculate  the  scattering  coefficient  variations  that  would  be 
observed  if  there  were  no  modulation  of  the  capillary  waves  by  using 
the  capillary  wave  spectrum  S(K)  in  the  theoretical  expression  for 
scattering — and  the  S(K)  used  can  be  that  derived  from  the  experiments 

121 


ami-. 


themselves.  The  difference  of  the  observed  radar  scattering 
coefficient  and  the  one  that  would  be  observed  without  modulation  is 
therefore  a  measure  of  the  modulation,  here  called  the  modulation 
index. 

This  instantaneous  modulation  index  is  then,  in  our  analysis, 
cross-correlated  with  the  instantaneous  record  of  wave  height  and  that 
of  wave  slope  to  obtain  new  measures  of  the  variation  of  the  modulation 
over  the  underlying  larger  waves  (previously  only  frequency-domain 
measures  have  been  used) .  Furthermore,  the  RMS  value  of  this  index 
represents  a  new  measure  of  the  size  of  the  modulation  of  the  capillary 


waves. 

Preliminary 

interpretations 

for 

the 

correlation  with 

wave 

height 

have  been 

made  here  for 

the 

sample 

case 

studied. 

The 

interpretation  of 

the  correlation 

with 

wave 

slope 

remains 

to  be 

developed. 

In  the  course  of  the  analysis  it  became  obvious  that  ''specular 
events'*  or  *'sea  spikes''  were  present  in  most  of  the  measurements.  A 
technique  was  developed  for  removing  these  from  the  data  used  in 
analysis  of  the  modulation  phenomenon,  but  a  study  of  these  anomalies 
themselves  will  be  important  as  the  data  analysis  proceeds. 

The  observations  used  in  this  preliminary  study  were  for  the  radar 
looking  in  the  upwind  direction.  The  same  methods  can  be  used,  for  the 
downwind  direction.  For  the  crosswind  direction  the  analysis  can  be 
simplified,  since  the  long-crested  assumption  (and  the  available  point 
measurement  of  wave  height)  precludes  assuming  any  slopes  relative  to 
the  radar  look  direction.  Crosswind  data  analyses  will  be  presented  in 
future  reports  in  the  series. 

An  algorithm  has  been  developed  to  convert  the  time  histories  of 
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the  wave*  into  spatial  dasoriptions  of  appropriate  parts  of  the  wave 
trains.  This  technique  can  be  need  readily  with  the  wave-height  and 
wave-slope  records,  but  its  use  with  the  radar  cross-section  has  more 
doubtful  validity.  The  wave  height  can  be  described  by  a  series  with 
terms  of  the  form  •*%*))  ,  but  the  radar  signal  is  nonlinearly 
related  to  this  series.  Further  study  is  required  to  determine  whether 
the  difference  between  the  observed  radar  signal  and  that  due  to 
uniformly  distributed  capillary  waves  is  sufficiently  small  that  a 
linear  approximation  can  be  made  that  would  validate  the  spatial 
Fourier-like  series  for  this  quantity. 

The  remainder  of  the  discussion  in  this  section  relates  to  the 
preliminary  results  obtained  by  analysis  of  the  sample  data  set  for  21 
September  at  40°  angle  of  incidence. 

9. 2  Discussion  of  Results  for  Sample  Data  Set 

The  preliminary  analysis  of  the  Noordwijk  1979  data  has  shown  that 
the  measured  scattering  coefficient  leads  the  wave  height  by  about  1.8 
sec.  The  full  period  of  the  dominant  wave  is  about  5.8  sec  so  1.8  sec 
corresponds  to  about  104°.  The  measured  scattering  coefficient  leads 
the  maximum  slope  by  1.0  sec.  Note  that  the  full  period  of  the 
dominant  slope  component  is  about  3.5  sec,  so  1.0  sec  corresponds  to 
103°  for  both  Hfl-  and  W-polarizations. 

These  phase  differences,  especially  the  phase  difference  of  the 
scattering  coefficient  with  respect  to  slope,  brought  up  the  idea  of 
modulation  index,  p(t)  or  A 5"  ,  the  ratio  of  total  scattering 
coefficient  (based  on  the  instantaneous  footprint  of  the  radar  antenna 
beam)  to  the  theoretical  scattering  coefficient  that  would  be  found  if 


capillary  wave  amplitude  were  uniform  over  the  large-scale  waves.  The 
theoretical  value  is  also  based  on  84(C),  where  C  is  the  instantaneous 
Bragg-resonant  wave  number.  The  RMS  value  of  modulation  index  is  0.876 
with  a  standard  deviation  of  0.349  for  HB-polarization  and  0.623  with 
standard  deviation  of  0.044  for  W-polarization. 

The  apparent  phase  relations  between  the  modulation  and  the 
dominant  wave  and  slope  components  are  illustrated  in  Figure  9.1.  An 
apparent  discrepancy  exists  between  the  relations  with  the  dominant 
''slope  wave*'  and  the  dominant  ' 'waveheight  wave'*  because  of  the 
higher  frequency  of  the  dominant  slope  (associated  with  the  larger 
influence  of  the  high-frequency  components  of  the  wave  on  the  slope 
than  on  the  height) .  Interpretation  of  this  point  remains  to  be 
established. 

Figure  9.1a  shows  that  the  largest  modulation  appears  62°  back 
from  the  crest  of  the  wave,  on  the  upwind  side,  but  a  relatively  strong 
modulation  peak  appears  to  be  153°  ahead  of  the  crest.  These  do  not 
appear  to  be  artifacts  of  the  data  handling  procedure,  so  one  must 
assume  that  they  are  real,  pending  analysis  of  other  data  sets. 

Figure  9.1b  is  a  similar  illustration  indicating  the  location  of 
the  modulation  peaks  relative  to  the  dominant  slope  component.  These 
appear  to  be  symmetrical  and  near  the  point  of  minimum  slope.  This  is 
consistent  with  the  location  of  the  maximum  modulation  on  the  back  face 
of  the  dominant  wave,  but  since  the  dominant  wave  and  slope  have 
different  periods,  one  must  use  caution  in  this  interpretation. 

As  indicated,  the  interpretation  of  the  locations  of  the 
modulation  peaks  is  confused  because  one  cannot  simply  use  the  ideas 
that  come  from  a  very-narrow-band  process.  Furthermore,  the  1-second 
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Wind  Direction 


figure  9.1 

Graphic  illustration  of  the  nature  of  the  results  of  the  correlation  study. 
The  locations  of  the  maximum  correlations  between  apparent  capillary  ampli¬ 
tude  (from  Ao)  and  wave  height  on  the  dominant  component  are  shown  in  (a) 
and  the  locations  relative  to  the  dominant  slope  component  are  shown  in  (b). 
The  apparent  locations  of  maximum  capillary-wave  amplitude  are  sketched  on 
large  sinusoids  representing  dominant  large-wave  components.  The  sketches 
below  show  more  quantitatively  the  apparent  distribution  of  capillary  ampli¬ 
tudes,  but  these  sketches  should  not  be  considered  accurate  as  to  amplitude 
they  are  intended  only  to  illustrate  the  nature  of  the  waves. 


sampling  interval  used  to  measure  the  wave  height  may  have  been  too 
great  to  give  a  good  estimate  of  the  slope  spectrum,  i.e.,  the  latter 
may  be  corrupted  by  aliasing  since  its  components  are  multiplied  by 
2,  thereby  enhancing  the  effect  of  the  higher-frequency  components  of 
the  measured  spectrum  that  are  more  subject  to  corruption  by  aliasng. 
Clearly  more  study  is  needed. 

The  study  of  modulation  in  the  spatial  domain  has  just  started, 
and  more  attention  to  the  basic  assumptions  underlying  it  should  be 
given  before  quantitative  results  are  presented  in  the  spatial  domain. 


9.3  Discussion 

9.3.1  The  Validity  of  the  Analysis 

The  validity  of  this  analysis  is  directly  related  to  the  goodness 
of  the  theoretical  scattering  coefficient  estimation.  Two  different 
problems  arise  in  estimating  the  theoretical  scattering  coefficient: 
first,  the  validity  of  the  estimate  of  the  K-spectrum  of  the  capillary 
wave  region  84(C) ,  and  second  the  effect  the  computer  calibration 
factor  which  differs  from  run  to  run  (Section  6.2)  has  on  the  true 
value  of  the  theoretical  scattering  coefficient.  The  $4(0  which  has 
been  used  in  preliminary  analysis  was  based  on  an  average  of  three  days 
of  data  in  the  windspeed  range  of  8.0  -  10.0  m/sec.  The  average 
calibration  factor  (ratio  of  the  mean  of  the  measured  scattering 
coefficient  to  the  mean  of  the  theoretical  scattering  coefficient) 
is  17.3  with  standard  deviation  of  4.8  for  HH-polarization  and  3.04 
with  standard  deviation  of  0.66  for  W-polarization. 
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9.3.2  Farther  Study  Needed 


The  preliminary  results  have  shown  that  the  idea  of  modulation 
index  can  describe  the  phenomenon  of  capillary  modulation  over  the 
large-scale  wave.  For  further  analysis  a  more  accurate  model  of  the 
K-spectrum  S^(L)  is  needed.  Note  that  in  preliminary  analysis,  each 
set  of  theoretical  scattering  coefficient  profile  has  been  normalized 
to  the  ratio  of  the  average  of  measured  scattering  coefficient  to  the 
average  of  theoretical  scattering  coefficient.  This  should  be  checked 
to  make  sure  that  this  is  the  best  approach  for  further  analysis. 

The  next  phase  of  the  analysis  is  to  test  the  idea  of  modulation 
index  for  different  conditions  (wind  speed  and  pointing  angle)  to  find 
out  how  the  modulation  index  varies  as  a  function  of  different 
parameters.  The  third  phase  of  the  analysis  is  to  see  if  the 
modulation  index  can  be  expressed  as  a  'simple*  function.  At  the  end, 
one  not  only  can  estimate  the  amount  of  modulation  and  the  location  of 
the  peaks  of  the  modulation  with  respect  to  the  slope,  but  also  can 
give  a  idea  about  the  'dominant*  shape  of  this  modulation  as  a  function 
of  wind  speed. 

Comparison  with  cross-wind  observations  should  be  especially 
interesting.  The  slope,  under  the  long-crested-wave  assmption  is 
always  zero  at  crosswind.  Viewing  the  capillary  wave  components  normal 
to  the  predominant  waves  may  also  provide  useful  insights  into  the 
capillary  distribution  process. 
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APPEND  EX  A 


The  cross-correlation  of  two  proper  finite  functions,  x(t),y(t), 
can  be  represented  as 


^  C  ^  ^  -T  b  d  T')  +  -  «  eO  5 

'T*  " 


*4  C°s(2n* r) *  %  S+tiJL* t) 


(A. lb) 


where  a,b  c.d  are  Fonrier  components  of  two  functions: 


XM=  ra„  C- i-.* (Z2!l  t)  <*.2> 

'Tt  ^  M 

y<*>:  rch  Cos  + Jt,  Sc^CUUlO  (a. 3) 
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The  Fourier  coefficients  are  unique  for  each  set  of  records. 

i 

The  quantities  and  fi  f,  have  a  probability  density  functions 

given  approximately  by  [Donelan  and  Pierson,  1981] 


f {Pm  P,*  -  (  ^Vj 


(A. 4) 


+  zm 


SC/;*// 


(A. 5) 


J2_  _ 


where  5  (/)  is  the  'true',  but  unknown,  cross-spectrum  of  the  random 


procats  which  is  assnaed  to  ba  approximately  a  stationary  Gaussian 
proeass. 

Equation  (A. 4)  is  a  chi-square  distribution  with  two  degrees  of 
freedoa  with  an  unknown  parameter, 

.  2 

Given  just  one  sat  of  tine  profiles  the  values  of  ^  nust  ba 
snoothad  over  frequency  as  in 


rtrl C  ^ 

Pm*  — —  Z"  P* 

2K  +  I  K 


or  more  generally 


(A. 6a) 


(A. 6b) 


where 


(A. 7) 


so  as  to  obtain  a  snoother  function  for  the  cross-spectrua  estiaate  by 
aeans  of  the  assumption,  which  nay  not  always  be  correct,  that  the  true 
spectrua  is  slowly  varying.  If  the  cross-spectrua  is  slowly  varying, 
then  the  values  of  will  be  approxiaately  distributed  according  to 
a  chi  square  distribution  with  2(2R+1)  degrees  of  freedoa.  Successive 
estiaates  will  not  be  independent.  Those  eleaental  frequency  bands 
that  are  2R+1  bands  apart  will  be  independent.  With  five  pooled 
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could  b« 


staple a  tad  a  triangular  moving  window,  the  value 
shown  as 


(A. 8) 


If  one  gets  an  average  of  10  individual  staple  records  of  *  then 

'  of 

equation  (A. 5)  san  be  estiaated  by  100  degrees  of  freedoa. 

For  a  chi-square  distribution  with  2  degrees  of  freedoa,  it  can  be 
shown  that 

(o./e  3  <  2  S  -91  )-O-70  (A.  9 ) 

Sh 


and  since  the  value  of  ^7  is  known  from  the  FFT,  it  follows  that 


(A. 10a) 
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Given  just  one  value  of  ,  the  value  of  the  cross-spectrua 

is  not  known  to  within  a  factor  of  58  at  the  90  percent  confidence 
level.  For  the  S  point,  triangle  moving  average,  10  different  saaple 


—  r 

records,  the  estimate  of  V  has  100  degrees  of  freedom  so 


P  ( 11- 13) 
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so  that 
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The  tree  cross-spectrum  is  known  to  within  about  ±  23 
or  s  range  of  47.9  percent  at  90  percent  confidence  level. 
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APPEND H  B 


Theoretical  Development  of  'Events*1  Effect  (based  on  Moore,  1981) 

For  the  theoretical  development  a  simplified  model  involving  a 
oniform  capillary— wave  distribatioa  is  assumed  (Section  5.4).  It 
further  is  assnaied  that  randomly-  located  specular  points  occur.  The 
waves  are  assumed  to  be  one-dimensional  traveling  in  the  direction  of 
the  wind. 

The  Bragg-scatter  component  is  made  ap  of  a  coherent  snm  of  the 
voltages  associated  with  scatter  from  each  individual  capillary 
wavelet.  The  voltage  received  from  an  individual  capillary  wavelet  is 
Vc.  The  voltage  associated  with  the  scatter  from  an  individual 
specular  point  is  Vp.  The  total  number  of  capillary  wavelets  in  a 
footprint  is  '77^  and  the  total  number  of  specular  points  is  • 

Assume  a  resolution  cell  with  length  L  and  width  V.  The  radar  is 
at  a  height  h  and  has  both  vertical  and  horizontal  beamwidths 
The  transverse  dimension  of  the  footprint  is 


(B.l) 


and  the  longitudinal  dimension  is 


L  -  -  f*  h  (D.2) 

since 

\l  —  CoS  0  (B.3) 
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Consequently  the  footprint  tret  is 


/)  X  ^  !l  (B.4) 

Co*  *fi 

The  total  voltage  received  from  the  capillary  waves  is 
proportional  to  the  nuaber  of  them  because  of  the  coherent  additions, 
that  is. 


Vr<  -  & 

Consequently,  the  square  of  the  received  voltage  is 


Note  that  the 
proportional  to 
given  by 


scattering  coefficient  is  an  ensemble  average  that 

z 


Vyc 


Consequently  the  scattering  coefficient 


is 

is 


<r,l  s/t  <  vrl >// y* y/s  «*•» 


There  i*  ■  constant  including  the  various  radar  paraaeters  and  the 
other  quantities  have  been  designated  above. 

Since  the  'events'  occur  at  random,  the  voltage  due  to  them  do  not 
add  coherently.  Rather,  the  power  must  be  added.  Thus,  the 
mean-square  received  voltage  for  the  events  is  given  by 
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(B.  8) 


where  is  the  ness  number  of  events  per  footprint.  This  nay  be 

less  than  1  since  many  footprints  aust  be  averaged  to  obtain  the 
enseable  average  of  (B.8).  Using  this  expression  we  nay  determine  a 
scattering  coefficient  for  the  events: 

<r  *  _  sM** 

/>  /> 

The  total  scattering  coefficient  is  the  sua  of  those  due  to  the 
capillary  waves  froa  (B.7)  and  those  due  to  the  events  froa  (B.8). 
Thus 


cT 


£L 


(B. 10) 


For  a  given  length  L  of  the  illuainated  area,  one  aay  obtain  nc  and  tip. 
n0  is  siaply  the  length  divided  by  the  Bragg-resonant  wavelength  on  the 
surface. 


/yi£  &  (B.  11) 
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Up  is  the  ratio  of  the  length  L  of  the  footprint  to  the  mean  spacing 
Sp  between  events 


Substituting  these  values  into  (B.8)  gives 
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which  way  be  rewritten  in  teras  of  the  capillary  wave  scattering 
coefficient  and  a  correction  term  as 


<r*=  4(<Vc>^y*»Wu  <£>*_ 

'  4  >  %  L  y 

or 

r  °=  ( '+  —  <r»r| )  <b.  14b) 

The  first  factor  in  the  second  tern  is  an  unknown  quantity  because  we 
do  not  know  the  ratio  of  the  nean  power  returned  for  an  event  to  the 
aean  power  returned  froa  an  individual  coaponent  of  the  Bragg 
scattering  capillary  waves.  The  quantity  is  also  an  unknown  at 

the  present  tiae,  although  it  aay  be  possible  to  find  a  value  for 
^  .  At  present,  however,  these  unknown  quantities  are  siaply 
coabined  into  a  single  unknown  factor  Q.  The  second  factor  in  the 
second  texa  involves  paraaeters  of  the  systea,  and  is  readily 
detexained.  The  third  factor  depends,  of  course,  on  the  angle  of 
incidenoe.  Thus,  we  nay  write  this  equation  as 

<r°-<rc%(/+  Q  (B.is) 

Although  the  value  of  Q  is  unknown,  we  can  coapare  results  for  tower 
and  aircraft.  Ezaaples  are  shown  in  the  tables.  Table  B.l  gives  the 
variation  of  the  ratio  of  the  scattering  coefficient  including  the 


effect  of  the  events  to  that  ignoring  the  events  and  including  only 
capillary-wave  Bragg  scatter  for  the  tower  experiment.  Values  of  Q  of 
100  to  50,000  are  presented.  If  Q  is  only  100  the  events  have  little 
effect  on  the  scattering.  However,  a  value  of  Q  of  10,000  gives  a 
difference  of  over  10  dB  with  and  without  the  events  at  20  degrees, 
falling  off  to  a  negligible  value  at  80o. 

Table  B.2  gives  a  comparable  set  of  calculations  for  an  aircraft 
at  3000  m  altitude  except  that  the  minimum  value  of  Q  is  1000  and 
computations  are  carried  out  up  to  Q  *  5  x  10-6 .  If  one  assumes  that  Q 
is  somewhere  between  10,000  and  20,000,  in  accord  with  the  observation 
that  tower  experiments  give  20  degree  measurements  between  11.7  and 
14.5  dB  higher  than  aircraft  experiswnts,  we  can  check  Tables  B.l  and 
B.2  to  see  if  this  is  really  the  case.  For  Q  “  10,000  at  20  degrees 
with  the  aircraft  the  difference  is  only  .23  dB  between  considering  and 
ignoring  the  events.  Tet,  this  was  the  case  where  there  was  more  tnan 
a  10  dB  difference  at  two  heights.  For  Q  *  20,000  the  effect  of  the 
events  is  still  less  than  0.5  dB  for  the  aircraft,  yet  it  is  more  tnan 
13  dB  for  the  tower. 

Thus,  one  can  use  this  to  explain  the  difference  between  tower  and 
aircraft  measurements  on  the  assumption  that  Q  is  between  10,000  and 
20,000  for  the  condition  of  measurement. 


1 


V 

20  dB 

30  dB 

40  dB 

50  dB 

60  dB 

70  dB 

80  dB 

100 

0.56 

0.23 

0.11 

0.06 

0.03 

0.01 

0.002 

1000 

3.77 

1.90 

1.0) 

0.53 

0.26 

0.10 

0.02 

10,000 

11.71 

8.13 

5.57 

3.60 

2.07 

0.94 

0.24 

20,000 

14.58 

10.79 

7.93 

5.54 

3.47 

1.72 

0.47 

30,000 

16.29 

12.43 

9.45 

6.88 

4.52 

2.38 

0.69  j 

50,000 

18.46 

14.55 

11.47 

8.72 

6.08 

3-45 

1.09  ; 

The  variation  of  the 
including  the  effect 
events  and  including 
for 

TABLE  B.l 

ratio  of  the  scattering  coefficient 
of  the  events  to  that  ignoring  the 
only  capillary-wave  Bragg-scatter 
tower  experiment 

1  1 

1 

i 

1  ' 

20  dB 

30  dB 

40  dB 

50  dB 

60  dB 

70  dB 

l 

80  dB  j 

1000 

0.02 

0.01 

0.004 

0.002 

0.001 

0.0004 

0.0061  j 

5000 

0.12 

0.05 

0.02 

0.01 

0.01 

0.002 

0.0005  I 

10,000 

0.23 

0.09 

0.04 

0.02 

0.01 

0.004 

0.001  j 

20,000 

0.46 

0.19 

0.09 

0.04 

0.02 

0.01 

0.002 

5x1  O'* 

5.76 

3.22 

1.82 

1 .0 

0,50 

0.21 

0.05 

5xl06 

14.58 

10.79 

7.93 

5.54 

3.47 

1.72 

0.47 

TABLE  B.2: 

A  comparable  set  of  calculations  of  the  racio  of  the 
scattering  coefficient  including  the  effect  of 
the  events  to  ignoring  the  events  and  including  only 
capillary-wave  Bragg-scatter  for  an  aircraft 
at  3000  m  a  1 1 i tude 
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APPENDIX  C 


Computer  Programs 

The  mala  goal  of  these  computations  is  to  calculate  the 
theoretical  value  of  0  based  on  instantaneous  slope,  determine  the 
spatial  data  profile  from  the  time— data  profile  and  ultimately 
calculate  the  proper  cross-correlations  both  in  the  time  and  spatial 
domains. 

The  processing  routine  consists  of  a  main  program  (EEN3J31,  two 
main  sub-routines  (CCA3G3  and  HH03L3 )  and  a  large  block  of  calculation 
routines.  Figure  C.l  shows  the  block  diagram  of  the  processing 
routine. 

(a)  EEN353  reads  the  first  set  of  raw  data  and  feeds  it  through 
HH03L3,  which  controls  and  provides  the  inputs  to  the  calculation 
routines.  Calculated  values  will  be  located  in  the  buffer  of  EEN353. 
When  the  first  set  of  calculations  is  complete,  the  second  and  third 
will  be  fed  through  HH03L3.  Calculated  results  will  be  averaged  by 
means  of  CCA3G3  routine.  Note  that  in  cross-correlation  calculations 
the  results  will  be  smoothed  and  averaged  in  the  frequency  domain  and 
then  will  be  transformed  to  the  time  domain  (Figure  C.2) . 

(b)  Low-pass  filtering:  FLTRR:  This  block  low-pass  filters  the 
wave  height  time  profile  by  means  of  an  FFT  routine.  It  transforms  the 
h(0,t)  from  the  time  domain  to  the  frequency  domain,  chops  all 
frequency  components  which  are  greater  than  .35  Hz,  and  transforms  the 
chopped  spectrum  to  the  time  domain  (Figure  C.3). 

(c)  Slope  calculation  SSLLPP:  This  block  takes  the  filtered  wave 
height  and  computes  the  spatial  derivatives  by  means  of  an  FFT  routine 
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(Figure  C.4).  Note  that 
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(>,-*)  _ 


A  n, 


feZj  Ut*  ^fr'  <r 

‘  J 

*  t-  * 

*<:  j  l*  £ii  e 

4‘  w7 


«t  =  &„"7  IhillM 

^  'h'X  / 


(C.l) 


(C.2) 


(C.  3 ) 


(C.4) 


i 
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(d)  Theoretical  (T  0  Calculation:  This  block  consists  of  three 

different  routines,  DIECO,  DDSTPP,.  and  PWRSK.  DIECO  is  a  standard 

routine  which  calculates  the  sea-water  dielectric 

constant,  £  ■  and  R  ,T  ,  Fresnel  reflection  and 

r  r 

transaission  coefficients  from  Q  ,  (antenna  pointing  angle), oC  •  (sea 
slope),  X  ,  (frequency),  ^(polarization  Ha  or  W) ,  sx  (salinity)  and 
T  (water  teaperature)  (Figure  C.3). 

PWRSK  is  a  routine  which  calculates  tne  value  of  S(K) ,  based  on 

(ji  »  (average  windspeed),  ^ ,  (frequency)  and  (9  ,  (local  angle  of 

incidence) . 

DDSTPP  is  a  routine  which  computes  the  theoretical  scattering 
coefficient  based  on  Chan  and  Fung's  model  [1977J 


r°(t)  -  r  i‘iHrp}Z  S( k) 

e 

and  also  ,  ^£*£2/  (measured  (J  0  over  theoretical  ). 

*i***fr  i/(0 


(C.  5 ) 


Final ly. 


it  calculate!  the  RMS  aodulation  'index*  based  upon: 


O'* 

row 


( b/v}  c 


TJfforjr 


( C.6 ) 


/"(i)  —  O* rt>tcx/  _  { 

CLf 


(C.7) 


(C.  8) 


where  N  is  number  of  samples  per  set. 

(e)  Spatial  Transformation  routine:  This  block  consists  of  two 
routines.  FRCFGN  and  XXFFCC.  FRCFGN  generates  the  proper  Fourier 
coefficients  and  makes  a  table  which  can  be  used  by  XXFFCC. 

XXFFCC  transforms  the  time  profile  data  set  to  a  spstial  profile 
(Figure  C.6) . 

“%(*/*)  r/fe  Tcc  e  (C.9) 

3(^V) 

=:  /^C  3  (C.10) 

c 

Note  that  if  the  transformation  is  not  made  this  block  can  be  ignored. 
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(f)  Cross-Correlation  routine:  This  block  consists  of  fonr 


routines:  MMOOTT,  ZZGGNN,  CCS323,  end  CCL3T3 . 

MMOOTT  subtracts  the  scan  from  the  sample  set  (Figure  C.7). 

(C.U) 

“  A/ 

1*  —  j*  (C.  12) 

ZZGGNN  adds  a  string  of  zeroes  at  tne  end  of  each  file.  This 
provides  an  interpolation  process  when  tne  FFT  routine  is  used  for 
oross-correlation  calculation. 

CCS3L3  and  CCL3T3  process  end  calculate  the  cross-spectra  by  means 
of  FFT.  The  result  of  averaged  nj  samples  of  cross-spectra  are  saved 
in  EEN353. 

(g)  The  final  cross-correlation  functions  are  provided  when  CCA3G3 
smooths  the  provided  values  of  cross-spectra  by  means  of  a  five-point 
moving  average  and  transforms  it  to  the  time  or  spatial  domain.  The 
smoothing  process,  five-point  moving  average  improves  the  estimation  of 
the  'true'  value  of  cross-spectra,  (Appendix  B) 

Miscellaneus  Proarams: 

(1)  Soclout : 

This  program  calculates  the  upper  limit  of  the  maximum  value 
of  scattering  coefficient  which  can  occur.  It  calculates  the  mean  of 
scattering  coefficient  and  adds  about  ^«^dB  to  the  mean  and  compares 
this  value  to  all  the  sample  points  of  measured  scattering  coefficient. 
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Figure  C.7 

Block  diagram  of  cross-correlation  calculation  routine  (CCS3L3) 
Note  that  the  final  transform,  from  frequency  to 
time  domain  will  be  completed  by  means  of  CCA3G3. 
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If  the  value  is  (raster  than  this  'threshold',  a  specular  event 
(Section  S.4)  will  be  replaced  by  a  weighted  average  of  its  S  adjacent 


points.  Note  that  for  90  percent  of  the  time  the  scattering 
coefficient  over  a  range  from  5.8  dB  above  the  mean  envelope  level  to 
11.9  dB  below  the  mean. 

(2)  MDTRNC 

This  program  calculates  MTF,  Section  3.0.  It  simply 
calculates  tne  smoothed,  averaged  cross-spectra  of  wave  height  and 
power  return  over  the  smoothed  averaged  power-spectra  of  wave  height 
and  mean  of  power  return.  For  obtaining  the  'stable'  and  smooth 
modulation  transfer  function,  the  result  will  be  smoothed  by  means  of 
seven  point  moving  average. 

4  fV)  H  V)  > 

!l> 

(3)  Coherence: 

This  program  calculates  the  coherence  function  for  power 
return  and  wave  height.  It  simply  finds  the  average  cross-spectrum  of 
wave  height  and  power  over  the  square  root  average  of  wave  height  power 
spectrum  and  square  root  of  average  power  return  power  spectrum. 

c\/)=  f 


(4)  FOPBT 

This  is  an  FFT  standard  Fast  Fourier  Transform  routine.  It 
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tikes  a  data  sat  and  will  transform  it  from  tine  to  frequency  or  frosi 
frequency  to  tiaa  domain.  Tie  direction  of  tie  transformation  depends 
on  iow  tie  variables  have  been  set. 


Y  (»  ) 


-  zr  «  * 


2T  Y  ( on)  e 

/tw 


bfl  h/rr> 

/V 


(C.15) 

(C.16) 


Call.FOUKr  (DATA,  Num.a.b.c.d) 


where: 

Data  is  tie  name  of  the  array 
Mum  is  the  number  of  Data  elements 

a,b  if  a«l,  b»— 1.  Tie  Data  will  be  transformed  from  time  to 
frequency  domain. 

a,b  if  a*l#  b«l.  Tie  Data  will  be  transformed  from  frequency 
to  tine  domain. 

o*0  if  imaginary  part  of  Data  is  zero. 
c*l  if  imaginary  part  of  Data  is  not  zero. 

4*0  if  No.  of  element  is  the  power  of  2. 
d«l  is  no.  of  element  is  not  the  power  of  2. 
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0013  PARAAETER  FC1  »30,FC2«1  J,FC>I3  ,FC4«1 0,  FCJ-11  ,FC4»14 

0014  PARAAETER  FC7«0I,FCS»09,FC?*14 

OOI3C  PARAAETER  FRQ»1 3. ,POL*l • ,  1HET  1*30., 00*21  *44 

001  AC  PARAAETER  FRQ*I3.,POL*I.,TETA-40., 00-20. 44 

0020  PARAAETER  FRO»14.1,POL  «1 .,  1H£IA-40.,U0*I.S1 

002 1C  <IF  POL* t .  POLARIZATION*!,  IF  P0L»2.  P0LARIZAII0I*9V> 

0022  PARAAETER  LF>0.0,MF>.33 

0023  1NTESER  S,U,U,l,F>OI,  S2,U2, 82,02, P0I2 

0024  REAL  TATATS>,XAIA(S),nn(2,O),ZC0SlO,S),ZSH<l,S>,FCOSC9,S>,FSIRCI,S> 
0027  REAL  0ATA(2,U),EATA<2,U>,flLT<2,Sl,CAtA<U>,SATA(S) 

0021  REAL  SL0P(2,S)  ,FREI(S)  ,ALPIA(S),  JUAK(2,S) 

0029  REAL  CRSTSU 1)  ,CISOSU< I)  ,CISTSU< 85 , CRSOSL  (U 1 ,  CRSDSU1U ) , CROSSC  I) 

0030  REAL  CRSOSL(II)  ,CRSLSU(U)  ,CRSIU(U)  ,CRSSS(  I)  ,CtSiH(U) 

0031  REAL  SATA(2,U)  ,ZATA(2, 0)  ,QAT4<I)  ,D8HIU) 

0032  S2«S 
0033  0*0 
0034  U2*U 
00 C  02-0 
0034  PaC*POM 
0037  PRINT, "POH»"  rP OM 2 
0034  REUIND  TCI 
0034  Ell  >0  FC2 
0014  EUIND  FC1 
0034  R EUIND  FC4 
0034  EH  40  FCJ 
0014  EUIND  FC4 
003 1C  SO  TO  997 

0019  CAIL  FILTER CFCT  ,FC2,S2r  LF,HF,  IA  TA  ,F  ILT ,  CA  TA ) 

0040  CALL  SLOPE  <  FC2, F  C3  ,S2,  SLOP  ,FE  I,  ALPHA,  JURK) 

0045  CALL  DISTRN<FC1,FC3,FC4,FC3,FC4,S2,FIQ, POL, TETA.UO, DATA, LATA) 

0047  GO  TO  917 

0049  CALL  XFUNCC  FC2,  S2 ,02, TATA,XA7A,2C0S,ZS IN, F COS, fSlN,FlNl) 

OOSOC 

0032  CALL  XrUMC(Fa,S7,02,TATA,lATAlZCOS,ZS!N,rC8S,FSZI,rilT> 

0033  CALL  XFUNC1FC4, S2, 02, TATA, XATA,ZCOS,ZS1«,FCOS,FSIN, PINT) 

0034  CALL  XFUNC(FC3,S2,02,TA7A,1ATA,2C0S,  ZSIN,FCO$,fSIN,FINT> 

0033  CALL  XFUNC(FC4,S2,02,TATA,2AIA,ZCOS,2S1N,FCOS,FS1I,FIIT) 

0037  997  CONTINUE 

0031  CALL  AN0«T(FC2,S2,RATA)  ' 

0039  CALL  ANOUT<  FC3,  S2,RA7A> 

0040  CALL  ANOOT (FC4,  S2.RATA) 

0070  CALL  ANOUK  PCS,  S2, RATA) 

0900  'AU  ANOUK  FC4,S2, RATA) 

0090  CALL  ZRE»< FC2,  2*$2tCATA) 

0092  CALL  ZREI(FC3,2<S2.CATA) 

0094  CALL  ZIGEN IFC4 ,2*S2, CATA) 

0094  CALL  ZIGEN(FC5,2»S2ICATA) 

0090  CALL  ZREMFC4,  2*$2,CATA) 

O099C  GO  TO  998 
0O99C  PRIIT,"P0«2',P0«2 

0100  CALL  CISSLK  FC2,  FC],FC4,FCS,FC4,FC7,FC8,FC?,U2,  92, OAT  A,  IAIA,CISTSL 
0 1050,  CRSOSU,  CRSTSU,  CRSOSL,  CRSOSU,CRS  SSL,  CROSS,  CRSL  SU  ,CIS8U  .CASSS  ,Cr<SW1,:AIA, 
0I07SSATA, OUR,OA!A,PON2, 

O10llHGT,SLP,SSO,SST,TWI) 

010  9  971  CGITIIUC 
01  IOC  STOP 
0115  RETUtN 
0121  EIO 


1.5.1 


.•’'KSKateti: 


,-lTRS 


.00  j  SUSSQuTIKt  FlLTENtTCl  ,r  C2.S,LF,iir  .DATA.r  !L  *  ,wATA) 

ooor  ifiTtGtR  re: ,rc:,s 

* 0 0 C  SEAL  LF,:if 

ooic:  ? akametcS.  s»i2* 

0000  REAL  DAT A ( », S 1 , F iL H 2 ,S J 
0021  SEA*  CATA.S) 

0030  SEAL  SUH, SUN 

0040C 

002  OC 

coio: 

0070  real  run 

0075  NUN-S 

0077  ?tIHT,"NUH*",NUN 

0081  FRE0*.3S 

0032  U)N*INT  (FSEB-NUN  ) 

0083  FEI»T,LU« 

0083  rtUN-l./FLOATINUN) 

0087  SUN-O. 0 

OOfO  00  10  K»1 ,NUH 

0100  READ<30,IOO)J,»HI,PUT,$IGO 

0110  OATAH  ,K)*iiHT 

01  IS  CAT*(8)*#HT 

0120  DATA(2,K)*0.0 

012S  FtLT(1  ,0*0.0 

0124  FILT  <  2,  It )  *0.0 

0127  SUMSUN-HATAH.iO 

0 130  10  CONTINUE 

0125  AV£«SUN/NUN 

0134  ? SINT, “SUN* , SUN 

0137  PRINT , "OUT* ,AVt 

VI 3a  CALL  F0URT(CAIA,SUN,1 ,*t ,0,0) 

u I  40  DO  11  K>1, (NUN/2) 

0142  if  d'.'jt.iun)  90  t4  444 
0144  FILTH  ,K*1 1*1 . 

0144  FILTC.NH 1*1. 

0144  r.'LTI  I, *««-«♦»  )»FUTI1,K»1) 


01  JO  ?ILTi2,«Uh*iC»l  '»-rILT(2.«+l ) 

0152  11  CONTINUE 

0153  44*  CONTINUE 
0154  FILTH, 11*1. 

0154  FILT(2,n«0.0 
0158  FILTH, NUN/2>»0.0 
0140  FILT(2,NUN/2)*0.0 

0200  10  12  K’l.NUH 

0210  DAT*(1,«)«0ATA(l,K)»FILTn,»)*NUN 
0220  I)ATA(2,K)»I)ATA(2,K)«FILT(2,«)*NUN 

0222  fill(T7X>«M»*mTrT5*T" 

0224  filL(2,IH«4u«*f  ilti2,k> 

0228C  1IRITEH4, 100  )K,0AIA(  1  ,K)  ,IATA(2,K) 

0227C  14, 101)4-1 ,  filtt  1  ,l> )  ,  fill* 2,4) 

0230  12  CONTINUE 

0240C 

0250  C 

0240  CALL  FOUNT) 0ATA.NUN,t,I,1,0) 

0245  CALL  F0URT1FUT, NUN,  1,1,1,01 

0270C 

0280  SUN«0.0 

0200  l>0  14  X-l  .NUN 

0205  5UN«SUN49ATA( 1 ,KI 

0300  UNI  TEH  2, 1000  IK.DATA1 1 ,  K ) 

0302C  -Pit*H  1,100)4-1  .FILTH  ,k),FIlT<2,K) 

0305C  UNITE  H  3. 1 0O!K,CATA(K> , DAT AH ,K),9ATA(2,K),FILT(K) 

0310  14  CONTINUE 

0315  PRINT, *SUN*, SUN 

0320  PRINT .  *UNTf  ILT* ,  SUN/NUf! 

0330  1000  F0NRAT1V) 

0333  100  FOIWATIV) 

0337  NETUNN 
0340C  STOP 
0350  £N» 


L 


y  1-.  ■ . mJfcMUifccim 
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L  5SLLPP 


00)5  SUBROUTINE  SLOPE ( F C2 ,FC3  .S, SLOP ,FR£3, ALPHA, JUNK ) 

00 1 0  INTEGER  S.FC2.FC3 
0)23C  DIMENSION  FILT(S) 

0)30  REAL  SLOP! 2,S) ,Fft£Q( S ) , JUNK (2,3) , huPHAI S ) 

0C50  INTEGER  NUN 
055  REAL  NUN 
0000  NUN*S 

005  NUN«I. /FLOAT (NUN) 

0070  PI*3. 1 415! ? 

0100  SUN»0.0 
0103  FQ‘.3S 
0147  REUINt)  12 
01 10  DO  30  r,-l  ,MJN 

0125  reaoiij.iooij.nm: 

0130  SLdF< I ,  K  J  *'JhT 
0140  SLO?u,id*0.o 
01 50  auN=ii..l*UHT 
OlaO  57  COnTIhUE 
0170  AVHIsSUN/NUn 
0140  PRINT, AVHT 
0135  PRINT, -SUM* .SUN 
0 1 90C 

020)  CALL  FQURTISLQP ,NUN ,  1 ,  - 1 ,)  ,0) 

0205  PRINT, "O.K.!  >!!!!!“ 

0210C 

022)  DO  101  1*1, NON 

024)  FREQII MFLQAT ( I-l ) /FLOAT (NUN ) 

0250  101  CONTINUE 

025 1 C 

0252C 

u2S3C  DO  112  «»1, NUN/2 

0254C  FILTiKl*1,/d ♦  <FREfl<  K J /r 0 j ** la  1 

025SC  FILT iK)*1 . 

0234C  SLOP!  1  , K) *r ILTIK )*SL0P 1 1  ,L ) 

0258C  SLOP! 1 .NUN-K-M  )«SLCP( I .NUN-K-M  )«f ILTCK) 

02aOC  SLOFi 2.NUR-K* 1 1*$L0P<  2,  NUN-KM ) <F1LT ( X ) 

02a2C  S^OP42.K»*FILT<K1 *SL0F(2,K) 

0264C  112  CONTINUE 
0245  JUNKd,  11*0.0 
0240  JUNlfl 2, 1 1*0,0 
0206C 

0270  DO  102  K* I, NUN/2 

0280  Z*(2.*PI*FREQiKM  1*2.  *F  I  *fRZ3li'.M 1  )/(T.S1 ) 

02S5  IFlFRESIKM)  .SE.  5.!  2*0. i 
02?0  jUNK ( 1 . K ♦ I ) *-Z*HUN  *SLO? iZ.KM ) 

■o::o  juiiKit-.rrn jxnras .jVGr nr 

C3I0  jun;:: i ,;,un-Kr i »-o*;:ii ,i;»n 
■2320  JUNK  t -K JUNK 1 2 , 1  *  I ) 

'.'32!  102  CC.'.TMJc 

0332  DO  Li.  ::*1  .H5.1 

'7234  sLCFd.KMfiL'N'SLJPi  1  ,K! 

033a  3LU;  i  2,11 !  *NI1N>SL0F  (2,0 
0338  SC  CONTINUE 
0J40  JUNK(:,iflU»/2)4l)-).9 
0343  CO  107  K*I ,NUH 

0347C  URITE(15,4CS)rRE0U)  ,5lCT  d  ,K)  ,2LJP< 2,ld  ,2 ( K ) , JUNK d 
034?  107  CONTINUE 


vlaOC 

O3’0  CALL  .*GUET(SL0Pt.;a.1.;.1,1,O) 
O250C 

7  2  T  0  CRcl  TCLRIt  wwKl« «  NUN, ! ,;,!,)/ 

?  4  C  C  2 

74  0  DO  1)3  1*1, NUN 

j42f'  alpha ; ; >  =a;an i junk ii , i > i 

0*23  AL*  HAi  I  ’  *  d 39 , ;F I  "■  ‘ALPHA! I ) 

103  CONTINUE 
. * 50  DO  10*  K*1  .NUN 
0*a)C  LA  I  TE'OR,  20)  IK- 1  ,SLJF.  I  ,K,  ,  JUKI! 
74a*  UR ITE 1 1 3, 300  IK , ALFnA IK J 
"470  10*  CONTINUE 
04*5  DO  F CRN*;. Vl 


0* jo  «oo  F::..v.r< 
;  *35  so:  rc pfisi' 
.  :t-j  II".  -  D  k.-.T, 

i .  i 


.  ' .  4  i  ,1  .2  • 


4  L  .*  N  7. .  L  1 ,  r  R  c  I '  4 
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> 1 .1  5 . 3 


Id,  JUNK  (2,  Id 


JopY  available  to  DTIC  dv'F 
permit  fully  legible  juprcdu 


loC 


v  I'USTrtn 


vwC  3u»f'wwii*c  4>»ii7R»'~Cl.r22*fL*..*C2frC*t,3.f*w,rCw,T«c7w,w.,0«i«,,.M»»# 

Wv7  I  )i  ”  2  C£a  r  C .  ».;23.rC,»,rCi,rCa,3 

000^  *  £  ml  rRu,rULt*Hc"A,uO 
00 1 v  COhPccA  vi&C0fA*»nn 

vOZO  ^chl  *2>2/*cfti>3i.) 

jOjO  h*l  vAT*i 2 , 5 j , £*7A *  2 . 5 » 

'.'035  REAL  SS.X,60 
003"  REAL  lAaOh 
.'04  0  T»:S6. 

0043  3«c -35. 

005QC  r  ?\£u-=  1 0 . 

0055C  TH£TA*40. 

wieC  00*5. ’e 

005 "C  jG*8.6 

0035 C  uOB* .  '0 

wS  *C  w.  ,».77 

OOovC  LAhfcA»2. 307o52 3 

Ovol  lANBn*jO./?A» 

wo  2  r  h  I  iNi  tr  5u«i»a0m  « i  n£T«  #f  uL  ,  w 

wo3  ?  i*4,*A7AN4  1  •  * 

0045 

004/  Z£o-*v .  0 

0048  SUH*6.0 

0070  A£UlHu  13 

0075  REWIND  36 

0060  DO  1 0  K*1 1 i>tU(l 

0696  RE  AK  1  3, 1 00) UU,  ALPHA 

0110  Rt80i36.bO6)g,»rti  ,r  ••  ,5150 

0120  6cTn*TntTA'*if 

Gl»2  OAAi4*»cTh 

lp  i  23  »»w»'n*-AL.;  nA»r*l/1  66  . 

0.25  dc"A*.;  **d£Th/» 50. 

0126C  F-ftlNTfi£TH 

012/  a*4  .  »F1  «S I K  \  *E7mj  '  »L«rtL>A) 

Oi  30  AA-AcOui  <J\  \  2  .  /  r  i  /  *  * » Z»r  */  wAAl'A  /  3  # ) 

O'»ov  CALl  criZOCN v T ,$Ac  ,r R5 ,  SANA,  McC0»6ii95Vf  62 , EMI 52 ' 

wl70C 

Ol  BO  FER5*<CAiSvftH>, *4  CAPS  C  ftH) ) 

0 1 35C  P  R InT ,r ERS 
0l?0  BMAlGG 1 6 < f ERS) 

0 1 92  t5s58*HL05  •  v  ( ' C  u5\ 5  c  T  A / )  **4 > 

Oi  *3C  ?  ft i  A • . r> 

0195  CALI  Shlv6|A.5ai 
0262  C wsS5 

0.6 jC  r 4  in T ,CC» uO 

02*6  62**hl05  1 6 \  i ./ 31  n\  vc  •  A  / ) 

0226  Ec5mA*d3*C^4»50 

06223  i'hTh(1 

0225C  URI TE (  19,  f  GO) K*  i ,  AA ,  BE* , CC  •  DIi ,  EE « EERS, DATA  1 1  ,  i» ) 

0235  Ip  i  AlF  HA.lI  .i.w)  cm  'AM  ,il)  *$I5C*Luo,«tr  nn, 

*237  lr  v  mlp  n*.6T  «g  •  •  ,  c*<A.  .  , i  *3  150/  2«i  <  *cr •  .*  i 

**  .v  2E3*IE3»i'*THi  1  ,K) 

0250  5un*5UH*c*TA< 1  tK* 

*255C  6RITE1  1 7, 1 06 >»*-i  , » v.  •  AwiG' Oi  51  ou > 

5*51  w*  i  Tc  i  1  i  f  *06  )K*i  r  1  HE  "A ,  Ac)  ha*  1 66  •.* )  I  »G*flA  f  mA,  pp ,  X ,  CC ,  M* ,  **«  *  a  '  ,*>.  ,  MEC3 

3240  10  CONTINUE 

o:ro  ::5*:eg,ploat(nu«) 

;  2 3 0  Sort* SUH/ FLOAT ■  i»U rt • 

029,)  p  4710*50*1 IE6 
.■£ ?2  r *  I  wT ,  RATIO*  ' r RATIO, 3uH, 2EG 
l 9 5 >1  « AT  10-  * . 

0-;4  "tMil'i1  « •* 

..**0 
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0003  SUBROUTINE  Dl 5TRN  t  r  Cl  ,FC3 , ?C4  ,?C5 ,  f  C6 , 5  , 7  R3  ,f  Gl  f  THt7«,  yv .i»A7A,£* .  » , 

«^w w7  m.tGcS  FCl •FC2«fi* «r CS'F 2**3 

wOur  At/tt.  f  nu»rWta,  umLTk»*v 
wiu  COrt?LcA  DicCw*.%«  #iri,/\ii,6V 

co:o  real  E*i$si2> 

0030  REAL  vfcTAi  1.5 J . £#?*■' 2 «6' 

VVji  At"W  44'A'UU 
Gvj/  ntm.  LMfl8A 

0040  7*23$. 

00*3  Snl*33. 

0C50C  rR£Q*10. 

0055c  THe;„a40. 

005  6C  UG*5.7b 
•>057C  tiO*8.S 
005 8 C  UG*6.2S 

oo59c  uco.rr 

0C66C  lA«06*2* 2076923 

0061  iA«DA*30./rR6 

00«2  P8IK7  ,r  *Q,innb*t  7HE76, FOi. ,  JO 

0063  f  I  ®  4 .  *  AT  n*>»  <  1  #i 

0O6S  MiM*S 

0067  ZcG*0«G 

0068  5urt=0.0 

0v/0  nt m3 NO  13 

0073  RtiitN#  30 

0060  M  iO  »»^i  f  /«uN 
0090  AEAK  13, 100>UU,  ALr-H* 

OHO  R£A9<J0,o00>JfUNT,?UTf5i$0 

0120  BETA*THtTA-»lr«A 

0122  &AiiA*BETA 

0123  nl?  n6-i»tf  nA*PI/l  90 . 

0123  BEiMfi^iETA/ t»0« 
ui2*C  fiirijfiti" 

Cl 2/  A* -ir.  *r  1*52 ft (BE  .  Ai /'u6r.*Ai 

C 1 30  AA*ALGGf  w<  i5*/f  Z/«<  (2*?I/lArt2'6>»*5>i 
01  40C  AA*AA*ALG&  *  0  n  CO j  i  i  c  i  A*  **4» 

vioO  LAwt  EnioCN^  Tr  S6L»r SJf  GAMA,  uItC0,R6 .  RV, R'2  f  CrtISS  ) 

oiroc 

0H6  GV*ftV-M  (CuS(BE7A>  >*-*2  )*  .3-m  t 1  >*0  *4  I .  -i  i  .*  BIECu/  >H ;  3i«(*£7A>  i  *0 

0198  FtRS*  *  CABS 1 5V >  '*tCA65tOV# ) 

Oi  9r  B6-AuG5 i Ot r EK$ > 
w200  CALL  S6\ttvfA,a6> 

02*2  CL-33 

02v5C  rr\i*T  f  CC»  Af  uv 

0 2  * 0  Gv*ni»Git6( i ./alNlvt tAi ) 

0220  EE*»n-*io»Cl  *i*B 
00223  BATAt I , * >*10.**£E 

0223C  URITci 19, 100>K-»  f  AA,BB  ,CC.  DC  ,  SEVERS  ,  DATA!  1  ,K) 

023*  IF t «Lr  AH.t 7 , y. 0/  cAmt » , r* )*5IG£+C0a  i *u/ «* > 

0227  Ir  \  ALf  nn««r  *  ,<J  •  v  »  *.*  i  6 t  •  f  »  «  s-iOu/CG»  '•  *»w? n6  t 

02***  2E3  =  I2G*DA7ii(  t  ,n  i 

0230  $'JH-$UA*t  AT  At  1 #rw 

0233C  wA I7£  v  i  •  vO » r»-  i ,  1  v •  *WIW«  i v \ a Iuu / 

v257  7£i  T  3*  »0v  , 7ii£TA# A^r'iiA*  I  80.  >'*i  « GArt*  f  6 A,  ,  iV ,  t  f  CC « sfi*9  i«ATA»  i  ,  K> » *>  I  ECO 

0260  10  CCriTINUE 

0270  2£u*2t«/FLOATirtU«) 

0230  3iirt*SUtt/FuuATi  NOrti 
».*2*0  S AT  IO*SUH/ 2to 

rp<ii*i  f  ,A»»7iO«5l>*»»2EG 

02 »3C  RAT iu*  1 . 
w-27  A  ni  wiNi*  i  * 

02 pv  i*M«*v.v 
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DOSTVV  (com*d.) 


■*p 


vi?7  DO  14  ri*1 

0340  0A7A(  i  ,* )-i»riO*&w7A«  .  ,rJ 
030 1C  &«fA( 1 ,K>*10. *ALuO IOC  DATA ( 1 ,K) 

G302C  tATA( 1 *K)*fO . *hwuu 1 Oi £hTA< 1 ,K) 

0303  UU*A*Sv  v£*7Ai  I  ,K)/  0*<m  » #r%*  /*i . ) 

304  $Un*50fUuu*uu 

03i7C  ilu<ci?iiUih)ailkrrt(lthj 

0333  JRITE  \  I  ’j,  1 40 it*-  i , 3A7AC  1 ,4) 

030?  liftITK  1 1  , 10O1K-1 , DATA*  i  ,*  > 

0310  tiftlTEi  14,  1 0 0 > K  — .  ,uw 
03 1 Z  iinXTtC  17,  »0v»n*  1 , Uu 
0310  1 3  CU»7it«ui£ 

0323  dMS*SGfiT ( SUA/NUH > 

0324  PRINT, *RA3  OP  4<t'«’*,RrtS 
0330  200  P ORrtATlVi 

0333  400  rOR1AT(tf> 

0337  IOO  FORNATvV) 

0340  400  FORMATi V) 

0342  700  PORHmmV) 

Oo« 0  500  POraAi  <  13 *3ii*F7«4t3X,r 7  •  4 ,3X*f v  » 4#3X  #?7 . 4, 3X #PT, 4 1 3*  »f  7.  , 3* » ?7.  A ,3a ,P7. 4 f  3X  ,r  *  4W f  3a»? 7 ,  i  < ,r .  • « 
v330  RtTuurt 
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l  i I ECO 


7430  SUBROUTINE  EHIGCHCUf.?  .SAL.PREv.THEIA,  OIECC.EH.RV , 

94?  0  1  R2,EK;5S) 

9500* 

9510*  PURPOSE 

9520* 

”530'  i.  IKE  E(;Z33KV  RCUTIuE  CALCULATES  T rtE  3ISLEC7A.C 

9540*  CC.'.STAXT  Or  SEA  WATER  AT  ANT  rKitSL£i«CY  lr, 

9550*  THE  KICXCUAVE  BAND  USING  TNc  ULc'c.  EXPRESSION 

9340 1  ( DEE Y£,F. ,  "POLAR  SCLECOlEs',  „.Cii£K.  PATS., 

9570*  VCL.  1,  1?4i).  THE  C OC“F I C TcTATS  OF  THIS 

V580*  EXPRESSION  ARE  FROM  THE  OOST  RECENT  EMPIRICAL. 

V 590*  MODELING  (MY  BiiDlIES/  L.  i.’LEIN  AND  C.  SUIFT 

7405*  (KLEIN, L. A,  C.T.  SUIFT,  AN  IMPROVED 

9410*  MODEL  FOR  THE  DIELECTRIC  CUSSIALT  OF  SEA 

9420*  WATER  AT  MICROWAVE  FREQUENCIES' ,  IEEE 

9430*  TRANS,  A«It*niAS  ARS  P RuF-SA".,  V0u.  AP-25( 1 J , 1 971 ) , 

9440* 

9450*  The  FSES.1EL  REFLECTION  COEFFICIENTS  ARE 

9440*  CALCULATED  USING  nETMD  BY  STRAITON 

9470*  < STRAITON, J . A ,  ELECTROMAGNETIC  THEORY, 

9480*  AC  SSAJ-HILL  DOOR  CO.,  I94W 

9490* 

97S0*  i.  THE  POWER  REFLECTION  COEFFICIENTS  AN 9 

9710*  SXISSSVITY  ARE  UETSRiilKE&t  MANUAL  OF  REMOTE 

9720*  SENSING,  AMES.  SOC.  OF  PHOT OGRANMETRY ,  CHAPTER  9). 

9730* 

9740*  INPUT  ARuUMEXTS 

9750* 

9740*  TEMP  •  SEA  SURFACE  TENPERATIIRE,  KELVIN. 

9770*  SAL  «  SALINITY,  0/00. 

9780*  FREQ  *  FREQUENCY  ,fiNZ. 

9790*  THETA  •  NORMAL  INCIJENCE  ANSLE ,  DEGREES. 

9100* 

9110*  OUTPUT  ARSUNENTS 

9I20* 

9830*  9IEC0  •  DIELECTRIC  CONSTANT  (COMPLEX). 

9040*  RH.RV  *  FRESNEL  REFLECTION  COEr FICENTS, 

9850*  •  H-  ANI  V-  POLARIZATIONS. 

90*0*  R2  •  POWER  REFLECTION  COEFFICENT, 

9070*  •  1  *  H-POL,  2  *  V-POL. 

9080*  EMISS  *  EMISSIVITY.  »  l.N-POl. 

9090*  *  *  2,  V-POL. 

9900* 

9910* 

9920*  CALCULATION  OF  DIELECTRIC  CONSTANT 

9930* 

994*  DIMENSION  K2C2r.ERI3SI7r 

9950* 

9940*  COMPLEX  VARIABLES 

9970* 

9980  COMPLEX  DIECO.RV,  RH 

9990*  COMPLEX  LUMMY  VARIABLES 

10000* 

10010  COMPLEX  0UMMY1 ,DUHMY2,IEP1 , TEK2 , COST , SO 

10020* 

10030*  DEFINE  DIELCIRIC  COnSTANY  *T  INFINITE  FREQUENCY ,EF In. 

10040*  AND  PERHITIVITY  OF  FREE  SPACE, E». 

1 0050' 

•:e i'j  DATA  EFIN, £0/4.9,  S.354E-I2/ 

10070* 

.5080  data  Fi/3.i4!59/ 

9  ft  *.  0  A  • 
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01  ECO  (cont'd.) 


1 9070* 
101 00* 
tone* 
10110 
10150 
l  0  I  40 
10150* 
10100* 
10170* 
10180 
10190 
10200 
10210 
10220 
10230* 
10210* 
10250* 
10260 
10270 
10280 
10290 
10300 
1 0310* 
1 C32  5  * 

ic::c* 

10340 
10320 
T  *3AC 
10370 
10383 
.:3*c* 

■  w«*v  W  * 

1  :•*;  s* 

*5423 
13430 
15443* 
1 0452 
1  3463 
1 041*3* 
TC433 
10490* 
10500* 
10510* 
10520 
10530 
10540 
10550 
10560* 
10370* 
10580* 
13590* 
10653 
!  o  1  3 
'5*23 


vb4 :  13 

3*23* 


1  ua.*v 


InlTIALlZE  CERTAIN  9A8IAHES 


T  *  leOf  -  273.16 
DELIn  *  25.5  ■  T 
U  *  2.  •  ?I  *  FSE3  *  1.0  £59 

EXPRESSION  708  IOMIC  CONDUCTIVITY  OF  SEA  UATER.CMD 

S  »  SAL  •  (0.182521  -  t .44I92£-3*SAL  *  2.09324E-3»SAL**2 
3  -l.2S20S£-7«SAl**3) 

I  •  2.0336-2  ♦  ).26aE-4«0ELTA  ♦  2.464E-6*0EtT»**2 
t  -SAL  *  (  1.844E-5  -  2 . S3 1 E— 7  *  BELT  A  .  2.53IE-3»D2L1a**2j 

CONO  *  S  *£XPl-OZLTA*l> 

EXPRESSION  F08  THE  STATIC  DIELECTRIC  C CHS T ANT  , -3. 

F  *  87.134  -  1.949E-!*!  -  1 .2766-2*7* >2 
8  ♦:.4916-4»T»-J 

A  *  1.0  *t.4!3E-5*5.UC  -  :.i5aE-2*SAL 
J  ♦*.2!£-5»SAL«‘2  -  4.232£-7*SAL**3 

ES  *  F  »  A 

RELAXATION  TME,  TAU 

TT  *  1.7486-11  -  6.034E-I3«T  ♦  1 .104E-I4*T»*2 
J  -3. ! 1 I£-17*T*»3 

IT  *  I.  ♦  2.2S2E-3*SAL*T  -  7.43ic-4*SAL 
3  -7.740£-»»3AL**2  ♦  1 . 16S£-8*S*L**3 

7AU  *  TT  .  8T 

SOL '.'2  iZDTE  EXPRESSION  FOR  DIELECTRIC  CONSTANT, "lECO. 

OuNHTl  *  CNFLX (1.4.  il*TAUl 
S3WIY2  *  C«PLX(O.3,C0NC/(U*£0)) 

TEFl  «  CNPLXIEFIK.O.O) 

TEr 2  *  CHFLX (ES-EFIN ,  0.3) 

DIECO  •  TEFI  ♦  TEFC/DUKNYl  -  ;c,1fl72 

CALCULATION  Of  fRESNt..  REFLECT  10*  COEFFICIENT, EiiJRV 

COST  «  CHPLXI  C0S( THE1A/S7.294) ,  0.0) 

SO  •  CSOSKDIECO  -  CSPLX ! i IN i  THETA,  5 7. 294)** 2 ,  0.0)  ) 

RH  •  (C5ST  -  SOI  t  (COST  ♦  it) 

RV  »  CISCO  •  COST  -  SO)  ( 3IEC0  •  COST  *  SO) 

CAcCotATC  l:  OUER  REFLECTION  CuErr  ICIEHTS.R21 '  *H,2*V> 
AND  EtllSalVITT,  EHlSSi ) *K,2*V1 . 


R2  ;  i  !  »  f.n  •  CCt.uGtAiit 

*2' 21  *  p. 1  CO.twutEVt 

ptf  1  t*  i ,. 

C  l'.  S3 . 1 ,  *  *  5, 2 , .  i 

*  w*9  •  *  >  •  w  w 


RZTUh 

too  J 


Copy  available  t~  DTIC  tie":?  •■'I 
permit  {ally  legible  iry-^lurt^ 
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L  XXFFCC 


01 00  SU E>K0UTI NE  XFUNCI FC , S, i , iAIA, XATA , ICOS.ZSlN f FCOS.FSIN.FIHT ) 

0290  INTEGER  S,0,FC 
yZIOC  FARAXETER  5*118,0*94 

0390  REAL  TATA<S>,XATA<SI,ZC0S<a,S),ISlN<a,S!,FINTI2,S) 

0310  REAL  FCOS(B,5),F5IR(0,S) 

0326  REAL  NUN 
0530  P 1*4 ,*ATAN< 1 • ) 

0333  G*9.S1 

0337  »0*9.0 

0340  NUI1»S 

05S0  NUN*I ./FLOAT (NUN 1 

0333  ii*2.*PI*NUN 

0339  PRINT, U 

0337  PRINT, NUN, FC, 5,0 

033?  REUIND  FC 

0390  80  13  l-I.NUN 

0393C  AA*C0S(U*1  .  *K) 

03A9C  AA«SIN(2.«PI*.17*IO*COSU.2*Ni»SlA(2.A*lt> 

0398  READ (FC,1 900 1VV, AA 
0370  TATAU)*AA 
0371c  PRINT, J 
0372  AQ*AO*TATA(K)*NUN 
0373  XATA<l)>«.0 
0373  13  CONTINUE 
0300  DO  13  K«1, NUN/2 
0383  FINTU  ,K)*4.( 

0390  FIIiT(2, 0*0.0 
0400  13  CONTINUE 
0403  REUIND  1? 

0407C  pri4t,*0000000?0?0" 

0410  to  It  1*1, NUN/2 
0413C  v*u*N*N/S 
0417C  V*FLOAT(K) 

0420  IG  30  J*1,NUH 

0427C  U*FL0AiCJ> 

0428C  U*1.3»J 

0439C  rCOS<K,J)*COSill*?*U> 

0440C  FSlNiK,J1*SIN<U*V*U) 

0444C  ICOSCK, J) *CuS \ 9*(*J ) 

0449C  ZSIN(K,J>*SIN(U*K*J) 

0447  READ) 1 ?, lOftOkAA.SI ,CC ,00 

0448  FC0S1K, J)*AA 

044?  F5IN(K,21*9« 

0430  ZCOSU,Ji*CC 

0431  ZSINvil.  J)*iiD 

•y 433C  liRfTE  1  Js,  1  COOT  C  ,'X 
v’490  20  CONTI/ilis 

0430  10  CONTINUE 
04»OC  prj«t,''????0?(?l?l?" 

0390  DO  39  K*  1 ,  Hurt/ 2 

03  IOC 

0520  >0  40  J*1 ,NUN 

0330  FINTd  ,K)*FiNT  ( I  ,R)  *NUN*TATA  <  J)*ZCOS<K,J> 

0340  FINT(2,K)«FINT<3,0*itUN*TATA<J>*ZSINU,J) 

0330  40  CONTINUE 

0590C  «*ITE(1*,:00)*,FiNT(1,K),FINTi:,i() 

0379  30  CONTINUE 

03S0C 

0590  DOSOl*l,NUR 
090SC 

OoiO  DO  90  J*1  ,iNUf1-2i/: 

0420  XATA4L  )*XATA(L  >  *2. *r INTI1 , J) »FC05( J, t )*2.*FIXT(2,jI*F$in(J,lI 
y«39C  PRINT, CATAiL) 

0433C  pri»t,/co*< j,l ;,f »l»< /.  i  i 

y.4 „  iV  CONTINUE 

044'C  f Klni , *  1 ! 1 ' 1 • 1 1 M 1 ■ ’ ! H ' ■ ■ 0* ' 

„t39C  rxlNl ,CAIA(Ll 
.ctO  IA!A(L;*4ATA(L)*A0 

0»33  IATA(L»*IATA(L1*riNTi 1 , ( NUrt/2i 1 .FC054 (NUN/2 , ,1 > 
ci.  •a;A(U»4AI4>.,*F INt,2,(NUn/2)  ).P5IN«(NUA/2),Li 

y»39  30  CONTINUE 

O940C 
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XXFFCC  (cont'd.) 


WO’VU 

6?!  PIIIiT,“»0*",flO 
0697  DEUIiiD  FC 
0700  DO  70  l* I, NUN 
0701C  U*4.»14.*L/F10»T(NUM) 

07 0J  U«30.*3.H 

O703C  K£*J(l2f!00)J,X,Z,iG,HH,yu,K* 

0704C  XcAv(99, 1000 >U,T 

0705C  1*3II<4.*M*?I-*.  17*.17*U/8)*C0Si4.34*U/5>*$ifci4.7i*U/'6> 

0706C  Y*C0Si 1.44*1) 

0707c  r*cos<v*u*u 
0708C  T*04T4(L) 

0709C  #ir»T-XATA(L) 

071 0C  ORITE ( I 5»i 0O)U , ThTAIL ) ,T .XATAiLl , OIF , r III ( I ,L< ,7 In 1(2,1 ) 

0< iii  niITc;l7,300)l,iA70(i.J , l,CAIA(L) >95,01^, HM ,ABS(CATA(L)~6G) 

0717  URITElFC.IOOOIO.Xii.Aa) 

0718  IFtrC.EB.i4l  ilRITCC  18, 1 OOO )U, XAT  A( l ) 

0710  70  CONTINUE 

0725  100  F0fm»T(F7. 3, »X, f ». 5, 4X, Ft. 5, 4X, f 9. 5, 4X.F9. 5,91,79.7,41,79. 7) 

0728  300  F0m7U3,4X,F9.S,4X,F9.3,4X,r?.5,<X,f4.5,4X,F?.5,4X,F?.J,4X,F».J> 
0727  200  F0»HA7<:3,4X,F9.3,4X,F9.3> 

0729  1000  rQRfl47(V> 

0730  RETURN 
07330  3T0P 
0740  EN8 
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L  r'4r.$/. 


OGtO  ivEACGTIi.c  3«'(wG,a,5S; 

0020  REAL  53,130,1 

0030  iri . 6. i.i.u’G KANO.  (UO.LE.  10. 1  >  SS*-4.*I 1  5**1061  OU)  *4106 I  0( .  I  782 » 

0040  1F< ( 1 3.LE.U0) .AND. (UO.IE.22. ) )  SS— 2.0f04»4LOOIO(I)»4LOOIO( .04*01 

oosoc 

0040C 

0470C 

OOIOC 

0090C 

OtOO  RETURR 
0110  ERD 


L  FRCFGN 


0290  PARAMETER  5*121,9*44 

0300  REAl  ZCOS (0, S) , :SI4i 9 ,S ) ,FCOS ( 9 , 5 ) ,FS i« i 1 , S ) 

0320  REAL  NUN 

Ojjfc  PI»4.  *4T«N(  I .  I 

viii  O'T.OI 

0040  KUfl*5 

uj50  ,/FlOh.  iRUR> 

0313  ••*2.*f3*n»X 
0324  FAINT ,U 

0410  DO  19  R*I,NUN/J 
0415  0*U*<iK/0 

>.'4I7C  V-ri.OATif.1 
0420  DO  <0  J*!,nON 

042/C  U*FLOAT(J) 

0420  J*30.»3.»J 

0430  FC5SU,  ,J)*C03(U»O*OI 

0440  FSIR1I.,  J)*SIn(U»V*U) 

0444  2COS(K,J1*COS(U»K«J) 

044*  ISIN.K,J!*$I*<y*».Jt 

04  4/  URITE(l»,lOtOlFCOSii9.J>,rSIR<*tJi.iC«iif.,j..:i;,4.i.,Jl 

0450  20  COUTIHUE 

0460  '9  COMTTNUe 

\j 499  i 090  F0nRA1,/1 

390  $70? 

OHO  END 


162 


L  NMOOTT 

0005  SUBROUTIRE  NN0UT(FC7,5.RATA) 
0010  INTESER  FC7,S 
0011  REAL  J'IATA(S) 

0012  n«n>t 
0015  SUR'O.O 
001*  REUIND  FC? 

0020  do  10  k<1.nu» 

0030  READ1FC7, I001J.SB 
0040  SUN«SUM«89 
0075  10  CONTINUE 
0077  vMu«/float(aua> 

0078  ?RINT,-M£AN«',U 
007?  REUIND  FC7 
0080  100  FORMAT (V) 

0012  DO  12  K«I,NUN 

0004  REA8(FC7.100)J,BI 

008A  >ATA(K)«IB-V 

0018  12  CONTINUE 

008?  REWIND  FC7 

00?0  DO  1 A  K«1 f NUN 

0092  URITE (FC7, 1 00 ) K  —  I , NATA(K) 

0094  U  CONTINUE 
0090  RETURN 
0100  ENI 


L  ZZGONN 

0010  SUBROUTINE  2RGEN(FC8,S.CATA) 
0020  INTESER  FCS.S 
0025  REAL  CATi(S) 

0027  NUI1«S 
0030  Jo  10  k«1,NUN 
0040  CATA<K)»0.0 
0050  10  CONTINUE 
0055  REWIND  FC8 
0040  DO  12  K«1 , NUN/2 
OO’O  READ(FC8 , 1 0 0 ) J , B I 
0010  CATAOI'BI 
0090  12  CONTINUE 
0092C 

0091  REWIND  FC8 

00?S  DO  14  K*  1 , NUN 

0100  URITE1FCI. IOO 1N-I .CATAIK) 

0105  U  CONTINUE 

0110  100  fornaltvl 

0120  RETURN 

0130  tad 
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*t  CCS3L3 


10  SUOtOU T 1NE  CISSU<FC2.FC3.FC0,FCS.FC4,FC7,FCfl,FC9,3,U,DATA,EATA,CtSISl 
203.CRSOSU , CRSTSU^RSGSL.CRSDSUtCRSDSL,  CROSS,  CRSLSU, CRSUli,  CRSSS,  CRSMI,  2A1A, 
3ttSATA,0UH,QATA,L0», 

3UH8T.SL7,  SOO.SGT,  (HOD) 

32C  PRINT,  •L0N\lOI 

44  1NTESEI  FC  2 ,  FC  3  ,FC  I ,  FC  3,FC* , FC7 ,F  C8  ,F C? ,S ,U 
St  REAL  9ATA(2,SI,£A1A12,S> 

40  integei  aui,ioa,ui 

70  REAL  CRSTSL(UI,CRSCI$U<U),CRSTSU(U),CtSQSldl>,CI5ISI(()  .CISISUI) 

80  REAL  CROSS) U ) ,CRSLSU(U) ,  CRSUU 1 S I , CRSSSv  9)  ,C tS/lli (U ) 

90  REAL  Z AT A ( 2, SI  ,  SAT  A  (2,0)  ,QATA  (U I , DUftlU ) 

1 00  LAI1*L0N 
no  seiiid  FC2 
120  REUINt  FC3 
130  REUINt  FC4 
140  REUINt  FC3 
130  RESINI  FC* 

143  OUN>S 
170  QH1«U 

180  PRm,n.O»-*,lOI,,NIR1*",NUN1"OUI«“,QI1l1 

194  10  34  X«1,N0h 

200  REA01 12,100)  J,  OHT 

210  IATA11  ,K)»UHT 

220  EATA<1,K)»«HT 

230  DATA<  2,1)  >0.0 

240  EATA12,K)»0.0 

242  HtT-«07*«HT*UHT 

203  3  0  CO  IT  IN  OE 

23 OC  30  UI1TE0 1 1,  IOO)  IATA11  ,K  1,  CATA11  ,K) 

2*0  REUINt  12 

271  CALL  COROL1 IAIA  ,EAT I,  CROSS,  ZAIA  ,SATA,  1AIA,  tuft ,NUI  ,11UI) 

200  10  35  K«1,0UN 

290  CRSUUOK)-CRSUU(K)  •CROSS)  I) 

300  33  CONIINOE 
0312  DO  330  «»1  ,N9H 
0300  SEADd  0,100)J.IIFF 
030*  DAT l<  1,1)  *D IFF 
0307  EAT A<  l,«)*DIFF 
0308  OAT A(  2, 1)  *0. 0 
0319  EATA(2,K  1*0.0 
0310  TMI-TMOD+IIFF.DIFF 
0311  330  CONTINUE 
0312  REUINt  14 

0313  CALL  COROLIDATA, EAIA, CROSS, ZATA, SATA, DATA,  JUN, NON, DUN) 

0313  00  333  K*l ,  OUlt 

031*  CRSIW<IO*CKSMII<«>  •CtOSSIId 

0317  33S  CONTI  HIE 

0319  REUINt  II 

32  0  80  1  0  Ot.NJN 

330  SEADd  I,  IOOU,SICT 

300  REAK  13,1  OOIJ,  SLJP 

330  IATA11  ,K  )»SIGT 

3*0  EATA<  1  ,K  )»SLOP 

370  SATAO  2,K)*0 .0 

380  EAT*(2,I)  ■0*0 

0380  S  IT  *S  OT  *S  IG  T*  SI  OT 

0388  SL  9-  SLR*  SLOP  *51.07 

390  10  CtNTIAUE 

400  CALL  COtOLIDATA.EAIA.CROSS.ZATA.  SATA,  DATA,  DUrl ,  NUN.UUMI 
410  DO  13  K»f ,QU« 

400  13  CRST SL  (I!  )»CRST  SL  (11 )» CROSS!  I) 

0403  IE  HI  AO  11 
430  (CU1 ND  13 
440  DO  29  K«1,NUH 
470  REAK  13, 1 0O IJ ,3107 
4 tO  DAT 4<  l.l)«SLOP 
4  90  EAIA<  1 ,11  )•  SL  OF’ 

500  IA 1A  (2  ,K  >*0.0 
514  IA!A(2,i:i«4.0 
524  77  CONTI  HUE 

534  CALL  COROLIDATA. CATA.CtOSU.CATA. SATA, DATA. DUn.NOA.aum 
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CCS3L3  (cont'd.) 


340  SO  47  Id.aurt 

331  49  CIS SS  IK  )»CISS5(K) ♦CROSS (K) 

340  REUIND  13 
37  OC 

380  0«  81 

37#  REAK  13,101)  J,  SLOP 
400  IEA»(12,1 OOIJ.UNT 

A 10  MTAI1  ,K)»$LOP 
A  20  EAIAI1 
630  DA1A<2,K)-0.0 
640  EATA12.K1-0.0 
A50  18  COin  IDE 

A At  CALI  COIOUDATA  ,EAT  A.  CROSS, 2ATA.  SATA ,QAT  A,  5IJN,  AUA,  OUI) 

A 70  00  99  K*1  ,QU41 

A  SO  99  CRSISU4  I)  •CISLSIK  ()*CROSS(() 

A 10  REUIND  10 
700  REUINI  12 
710  10  14  1=1, HUM 
720  ME  AO  (1  0,  lOOIJ.SIGO 
730  REAK  12,10O)J,UlfT 
740  0  4TAH  ,KJ»SISO 
730  £ATA(  1,1)  *UIT 
7A0  OAT  4(2,11 *0.0 
770  EOT A(  2, 0  4.0 
0773  SGO=SGO+SIGO*SIGO 
7 tO  14  CONTINUE 
770C  ' 

800  CALI  C0R0KD4TA, EAIA, CIOSS.ZATA, SA.IA, GATA,  KJO. NUN, OUM) 
810C 

820  00  13  X«1  ,QIR1 

830  1  3  CRSOSU  IK  )=CRSOSU (K  !♦  CROSS  <K) 

840C 

830  REUINI  11 
8A0  IE  IINI  12 
87  OC 

880  DO  1 A  K-1.NIM 
890  READ<I1,100)J,SIGT 
700  CAO<12,1IO>J,U«T 
910  IAT441  ,K )  *SIGT 
720  EATAI1  ,ld»UHI 
730  EATA<2,K)=0.0 
740  DA  TA  (2  ,K  1=0 . 0 
930  1A  COniNUE 
7A0C 

770  CAU  CGIOUIATA.EATA, CROSS. 2AM, SATA. QATA, OUI. AUI.DUI1 
7  IOC 

970  K3  17  K=I,QUM 

1000  1  7  CISTSU<«)«CI5T31KI)*CI(JSS<K) 

1  01  OC 

1020IESINI  11 
1030  REUINI  12 
1040  REUINI  10 
1030  RE  OX  ID  1 3 
1060  DO  II  K>I,*UK 
1O70REACK  10, 1  OO  U  ,S  ISO 
1  OGO  IEAIH1  J,lOO)J,SlOP 
1  070  OATAU,IU«STGO 
1  100  EATA4 1,1)  •SION 
1110  IATAI2. I>*0.0 
1120  E4T4<  2.1)  *0.0 
1  130  10  CONTINUE 
I  I40C 

1  1 50  C«IL  COROLIDAlA.COTi,  CROSS.  2A M.GAr 4.  OA M.OUa.min.  0U«) 
I  R.OC 
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CCS3L3  (cont'd.) 


540  DO  49  K»i,0UM 

SSI  49  CRSSS  IK  )*CRSSS(K  )+CRQSS(K7 

SAO  REMIND  13 

57  OC 

580  01  81  K-t,«JI 
591  RCAK  13,101)  J.SLOP 

6oo  ka»<i2,i  oou.urr 

6  tO  DATAd  ,U)»5LSP 
620  EATAd  (KIbIMT 

630  OATA(2,KI»4.l 
640  EATA(2,K)-I.I 

631  e  coitiiue 

661  CALL  COROL<DATA  ,EATA,  CR0SS.2ATA,  SA  1A.QAT A.XJI.IUA.OUA) 

670  DO  99  K*1 .QUH 

610  99  CRSLS04  It)  <SSLSU<  ()+CRUSS(R) 

690  REUINO  10 
700  REUINI  12 
711  90  14  R»l,NUh 
721  lEADdl,  I001J.SIGD 
730  ftEAK  i2,i«ou,uirr 
740  OiTAd  ,K)»SI50 
7S0  E4T4(  1.R)  >UI4T 
76 I  DATA(2,I)*0.0 
770  EAT 4<  2,  R)  *0.0 
0773  SG  *■  SG  0+  SIG4»S  1G0 
710  14  CGITIMJE 
79 4C"  * 

SOI  CALL  CO  ROL<  DATA,  LATA, CROSS  ,2  AT  A,  SATA  ,□  AT  A,  DUI.NUII,  OUH) 
8I0C  ... 

820  00  13  X«1,QW 

810  IS  CRSOSUOd-CRSUSLHKHCRDSSIK) 

840C 

830  REU04I  11 
861  IE  UNI  12 
87IC 

880  OO  16  K-1 ,NIH 
890  KEA04  11 , 1 0O  )J  ,S  IG  T 

910  REACH  2,1 10) J,URT 

911  |ATA(  1  tK )  »S  IGI 

920  EAtAd  ,ld»l*4T 
931  EAT  A<  2,  R)*l.  I 
941  M1A(2,K)»0.0 

931  16  COITINUE 

96 1C 

970  CALL  COIOL<SATA,£ATA, CROSS. IATA.SAU, GITA, DUI. «UI,OUI) 
9D0C 

991  IQ  17  K-I.OUh 

1040  17  CRSTSUt  I)  «CRSISU(R>*CR0S5(K) 
tltOC 

1020IEUNI  11 
1010  RCUD4I  12 
1040  REUINI  10 
1051  REUIIO  1] 

1061  00  11  K»I.«UR 
1 070RCAD1 10»  1 40  IJ,S  IGI 
1  0S4  IE AIK1  1.140)  J.SLOP 
1  J90  3A!A<  t,K)«SICO 

1  101  EATAC  1,1)  »SL09 
1114  1ATA(2.1)  *0.0 

1120  EAT A<  2. * >  «0 .0 
I  <31  1G  CONTI  RUE 
1  MIC 

mo  CALL  C3R0K0A1A, CATA. CROSS.  IAIA, SATA. IA1A, ClIJA. riw.ouiu 
'  I<.0C 
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L  CCL3T3 


10  SLUROUTINE  COOOl<  tAIA  ,c»T4,  CROS  S,  IAT 4 ,S4TA, 3AT4,  SUIf  >UAr  OUI) 

30  INTEGER  BUN, NUN 

33  LUN>NUH 

33  JIHCUN 

40  REM.  RUN 

43  REAL.  SATtC.SUl) 

30  REAL  ZATA(2,N1»1),D»TAi2,NUf1>,E4TA<2,»U8)  .CSOSSIOUM) 

40  REAL  QATA(QUH) , SUN (SUM) 

70  NUN-1. /rLOAT(NUN) 

80  «  1  2  H-l,  BUI 

82  SATA11 , KI-O.O 
84  SATA  (2  ,K  1-0.0 
90  12  CONTINUE 
IOC 

102  10  880  K*1 ,  NUN 

IOC  UlITEd  0,1  000)  DAIAd  ,ICI , MIA (2 ,l< ), EATA d  ,1) , EAT* (2, 1 ) 

104  980  CONTINUE 

110  CALL  F0LHT(O*T»,lU«,  1,-1, 0,0) 

120  C.'.LL  FOURI<  EATA ,L IW,  I , -1 ,0,0 1 
130C 

140  »  10  1-I.NUH 

130  Z»T*<  1,«).(D»TA<1,»)«€»T»(  l,l)«OAT*C2,l)«£*TN<2,K)l/LUN 
1  40  ZATA(2,K)*<0ATA(2,K)»EAIAd  ,K)-JATAd  ,K  )»EATA(  2,K)  1/LUM 
170  10  CONTINUE 
110  DO  112  l-I.NUI 
184  CIOSSIId-IATOd  ,K) 

188  CUSS (NUI*td>Z4T4<  2,0 
200  112  CONTINUE 
210  RETURN 
220  ENI 
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L  CCA3G ' 


10  SUBROUTINE  CRSLINI DATA.EATA, CROSS,ZAIA,SATA ,QAIA,DUH,NUN ,OUH> 

39  INTEGER  QUN.HUN 
33  LUN*NUN 

33  JUH«OU» 

40  SEAL  NUN 

43  REAL  SATA ( 2 , SUN » 

39  REAL  ZA1A(2,NUH),DATA(2,NUN>,EAT4<2,NUII),CIOSS(QUH> 

99  REAL  DAT A (DUN) f SUN (OUR) 

A3C  PRINT,“NUN»-,NUN 

70  NUN»1 ./FLOAT (NUN) 

71  P*. 33333333 

80  DO  12  K*1,QUN 
82  SATA(1,K)*0.0 

84  SATA(2,K>*0.0 

70  12  CONTINUE 
92  DO  10  K*1 ,NUN 
94  ZATAO  ,K)*CROSS(K) 

94  ZATA<2,I(>*CR0SS(NUH«K) 

90  10  CONTINUE 

100  DO  23  K»3,  (NUN/2)*2 

102  ZATA<1,K)*P«(f*(ZATA<I,E-2)*ZArA<l,K«2)l*ZArA<l,*)*2.«R»(:ATA(l,E-l) 

1 93t»ZATA( 1 ,K» I ) ) > 

104  ZATA(2,K)*P«(r«(ZATA(2,R-2)»ZATA(2,K«2mZATA<2,«x2.«P*(ZATA(2,T-l) 

1 0S1»ZATA(2,KM ) ) ) 

107  23  CONTINUE 

109  ZATA11 , 1  )«.3«<ZATA(  1 ,  1 >*P*ZATA( 1 ,3)  *2.  «P*ZATA(  1,21) 

110  ZATA<2,t>*.3«(ZATA(2,1XMZAtA<2,3)*2.»P»ZATA(2,2>> 

111  ZATA1 t,2)»(3./8.)»(P*ZATA1l,4>*ZATA(l ,2)  *2.  *P  MZATA1 1 , 1 )  ♦  Z  AT  A  <  1 , 3) ) ) 

112  ZATA<2,2)*(  J./8.  >*<P«ZATA(2, 4HZATA12, 2>*2.*P<(ZATA(  2, 1  )*ZATA(2,J)I) 

1 14  ZATA< I , (NUN/2) >*.5»(ZATA( I ,  (NUN/2)  XP«ZATAt I , INUN/2>-2>*2.  »P*ZATA1 
11311, (HUN/2)-!)) 

117  ZATA(2,<NUN/2)>».34<Z9TA<2,(NUH/2))*P*ZATA(2,(NUN/2>-2>*2.*P»ZATA( 

11832, (NUN/2J-I)) 

120  ZATA(l,<NUN/2)-l)«<3./8.>»(ZATA<1,(NUN/2>-IXP*ZATA(1,(NUN/2)-3)*2.»P*UATA 
12211 1,(NUH/2)-2)*ZATA(1,(NUH/2)))> 

124  ZATA(2,  (NUN/2 )-1  1M3./8.  )•( ZATA( 2,  (  NUN/2 ) - 1 )  *t *ZATA( 2,  ( NUN/2 >-3 M2. *9* < ZATA 
1261(2, (NUN/2) -2 >♦ ZATA (2, (NUN/2) )> ) 

130  DO  ISA  K*l, NUN/2 

132  ZATA ( I ,KM  NUH/2 ))*ZATA(1,K> 

134  ZATA(2,K*(NUN/2) >»-ZATA(2,K) 

1 36  ISA  CONTINUE 

172  DO  43  K=1 , NUN/2 

174  SATA(I,K)«ZATA(1,K) 

17A  SATA( 1 , QUA* 1-K)* ZATA ( 1 , NUN+1 -X ) 

178  SATA(2,0UN«1-I(MZATA(2,NUN»I-K) 

ISO  SATA(2,KMZATA(2,K) 

192  <3  CONTINUE 

240  CALL  rOURT(SATA,JUN,l,-17i;fr- 

2S0  DO  IS  K-1.QUN 

300  13  DUH(K)*$ATA( 1 ,K) 

310  100  FORAAT (V ) 

320  DO  16  K*1 , GUN/2 

330  CROSS (K)*DUN(K+(QUN/2)) 

340  CROSS (0UN/2+K)  sDUN(X  ) 

330  16  CONTINUE 
3A0  RETURN 
379  END 
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I  SFCLQUT 


0010  PARAMETER  NUN»I28,J0*. 14*6*6*7 

00 1  SC  PARANETE*  IR5H0LJ*"  10.0 

0020  DEAL  C AT A <  NUN) , EAT A <  NUN )  r  OAT  A  (NUN) , I AT  A (NUN) 

0423  REAL  OATA(NUH) 

0030  INTEJER  SAT  A  <  HUB  > 

0032  INTEGER  AA 

0033  DO  78  J»l,5 

0034  SUNS»0.0 

0034  SUNP'0.0 

0037  $UNST*0.0 

0031  SUf1PT»0.0 

0030  SUNPT *0.0 

0041  REUINI  30 

0042  10  3*  1*1 ,NUN 

0044  «EAI<40,100)AA,I8,CC,3D 

0046  URITE(3O,10O)AA,JB,CC,8D 

0048  34  CONTINUE 

0040  REUINI  30 

0030  JO  10  I* I, NUN 

0060  HEAD (30. 10O)SATA(K),CATA(K),£ATA(K),DATA(K) 

0074  ZATA(K)*DATA(K) 

0073  OATA(K)*CATA(K) 

0077  SUNS*SUHS* 1 0 . • AL08 1 1 ( IATA<  K) > 

0084  10  CONTINUE 

0000  USR0LD*<;UXS/F10AT(N8H))«6. 

0003  SUH$*0.0 
0100  to  20  Cl, NUN 

0110  IF< lO.'ALQGI OiDATA(X>  > .LT.TRSHOLD)  SO  TO  20 

0120  U1*4.0 

0123  U2-0.0 

0130  U4*0.0 

0133  U3>0.0 

0134  IF (8.EQ. 1 )  GO  TO  344 

0140  IF ( 10. *ALOG10(DATA(K-I ) ) -LT.TRSHOLD)  U2*2.»08 
0143  IF4K.EQ.2)  GO  TO  344 

0130  IF(IO.*ALOGIO(DATA(«-2)).LT.TJSHOED)  «I*QQ 

0133  IF(X.EO.NUN)  GO  TO  334 
0134  344  CONTINUE 

0133  IF(10.*AL0610(DA(A<K«I)).L1.IRSHOLDI  U4*2. >00 

0163  IF (K.EQ.NUH-1 )  GO  TO  334 

0140  IF(10.|ALOG1O<DATA(K*2)).lT.TISHOlD)  US*JQ 

0163  334  CONTINUE 

0170  UU*U1 »U2+U4»U5 

0173  IF(UU.EQ.O.O)  FRINT ,K 

0180  DATA<K)«(U1*DATA<K-2>*U2*JATA(K-1 )  »U4*  JATA1K*  I )  *U5*DATA<  C2  >  >/UU 

01  TO  EAIA(N)*(U1*EATA(K-2)*U2*EATA(«-I)*U4*EAIA(«M )»U5»EAIA(K»2) )/UU 

0200IF<DATA(K).N£.ZATA(X))URITE(3I,100)K,1 9. *41061 0CZATA4K) 1 f  tO.*ALOCIO(DATA(K>3 

02033, I0.*ALQGI0(8ATA(K) ),10.*AL0G10(EATA(K)),TR5H0LD 

0214  20  CONTINUE 

0212  DO  43  K>I,NUN 

0214  SUNS*SUNS+ZATA(K) 

0213  SUNP»SUNP*QATA(I() 

0214  SUNST*SUN$T*DATA(K) 

0217  SUNP1*SUNPT+EATA(K) 

0218  43  CONTINUE 

021?  SUNS*IO, *AL0G1 0  <  SUN S / F L 0 AT (  NUM  J ) 

0221  SUNP* 1 0. *ALOS 1 0 ( SUNP/FIOAT (NUN ) ) 

0222  SUNST*1 0. •  AlOGIOt  SUNST/F10AT1NUN1 ) 

0224  SUNPT*IO.*ACOG10($UNPT/FLOAT(NUN1) 

0223  PRINT , "TRSHOID*  .TRSHQLI , *SUH$" , SUNS, “3UNP’ ,SUNP, "SUNS T“ ,SUN$T , ’SUNPT" , SUNPT 
0230  DO  30  K*1,NUH 

0240  URITEC33, 1 OO)SATA(K) ,CATA(K) ,EATA(K ) ,DATA<K3 
0243  UIITE134 , 1 OOIK ,OATA (I) 

0230  30  CONTINUE 
0233  UKITEC21 , 1 OO) J.TRSHOLD 
0264  78  CONTINUE 
02*0  100  FORNAT(V) 

0270  STOP 
0280  EN8 
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•L  MBTRFNC 


08 

0010  INTEGER  B,T,RAU 
0020  INTEGER  NUA.AVE 
0030  PARAMETER  V»I28,U*128 
0040  PI*3. 141393 
OOSO  Q*. 33333 
0040  NUM*U 

0070  PARAMETER  5*1280 
0030  REAL  CRH<V),MORM(V) 

0090  DIMENSION  ERH< V) , AVRFUN ( V ) 

OtOO  DIMENSION  DATA(2,V),£ATA(2,V) 

OHO  DIMENSION  8ULK  C 2  ,S > 

0120  DIMENSION  AVCRSC.V) 

130  DIMENSION  PURSf EC(V) , PUS < W ) , GAAS(U) , CROSPURIV) , FI  IV) 

0140  DIMENSION  fUSPECm,RfSm,CNRENCE<V>,fRfQ(y> 

0150  DIMENSION  GMASPECUD ,CR0SPEC<2,VI ,RFUNCC2,VI .RFSPECtV) 

0140  DIMENSION  AORFUNC(U) .AUFROtU) .AVFIIU) ,AVPUR(U> .AUGAAIU) 
0170  DIMENSION  AVCSPURUJ) ,  AVC0H1 U ) 

OI80C 

0190  PRINT,  NUMBER  OF  PRCESSED  TERMS  TO  IE  A0ERAGED/1AUE*?) ' 
0200  REAO,AVPTIN 
0210  AVE*AVPTIN 
0220C 

0230  SUMPT*O.0 
0240C 

0210  DO  111  K s I , S 

0200  READ! 1 4,900) J,UHT, PUT, S IS 

0270  PRTN*PUT 

0290  SUMPT*SUNFT*PRTN 

02*0  PRTN.PUT 

0300  BULK <1,0 *FR7N 

0310  8ULK(2,K)*UHT 

03:0  111  CONTINUE 

033CC 

0340  DO  124  J*1,N1)M 
0350  AURFUH(J)*0.0 
0340  AUFI(J)*O.0 
0370  AUPURJ J>*0.0 
0380  AVGiiAl J ) «0.0 
0390  AVCSPBRiJ 1*0.0 
0400  AVCRSH  ,J>«0.0 
0410  AVCRS (2, J) *0. 0 
0420  AVCOHIJmO.O 
0430  124  CONTIHUE 
0440C 

04JC-  AU£PT«SUHPr .'S 
04)0  PRINT, AVERT 
0170  DO  115  r-I.AUE 
130  RAU* if- 1  mNUH 

01T0C 

0200  DO  113  CH.ilUM 

•3510  DATA!  1  ,K  1  -  DuLl.U  ,  (RAU  +  K ) ) 

0530  DA7A(2,i!/*0.0 

0530  EATAi 1 ,K * *BULK ( 2 , 1RAU*K ) > 

:-3 40  eata ( :,i; > »o.o 
0530  113  CONTINUE 

;3.;c 

03*0  CAL2  FOuR7(DATA,NUI4,  I ,  *1 ,0,0) 

33 30  CALL  FOUST ( EA TA ,NUH , I , - 1 ,0,0) 

0"  DC 

0.-  0  DO  i::  1*1  ,n’Urt 

0  GMA3rECC'siEATACl,I)*EATA(l,I)*EATA(2,I)#EATAI2,I)  )/NUN 
0)20  ?  'JA5PEC  ( I  )*<  DATA!  I ,  I )  »DATA( ! ,  1 1  f  DATA!  2  , 1 )  •DATA  I  2  , 1 ) )  /  RUM 
0:30  CR03FEC{1,I)*(»AIA(l,I)*EATA(t,IHDATA<2,I)*EATA(2,I)l/!1UN 
3)1 j  C533f ECi2,1 1*(DATA(2,I ) >EATA( I , I)-DATA( I , i >  *E A T A 4 2 , 1 ) ) /NOM 
3)30*  ? U3r EC  1 1 > •PURSf EC II ) 
umO  122  CONTINUE 

*oc 
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MDTRFNC  (cont'd.) 


OoSO  DO  211  J*1,flUH 
0  K*J-I 

O'OO  FSEO(J) ‘FLOAT  <K)/FI.O»T(NUn) 

'tO  211  CONTINUE 

cr:o  do  ii*  d*i .nun 

730  AVCSSil ,D1*AVCRS(I , B) ♦CROSFECt I ,S> 

'40  AVCRS 1 2,B) «mVCSS(2 ,D) ♦CROSF'EC 1 2,0) 

0730  OV?UR( D> *HVr'JR(DI ♦PUS SPEC! D) 

"iO  AVGrlAt  9>  *hVGHAC  S>«GMASPEC  CB  > 

0770  117  CONTINUE 
0780  113  CONTINUE 
•5790  PRINT, '•0C0r»O' 

0 S 0 0  AVSC»0.0 

OCtC  DO  1 1 6  1*1, NUN 

0020  hVFRG  i  I  *  *F?.EO  ( I ) 

Oi*5  A'.'P'OSd  .*A»P!IR(I>/4VE 
OS  - '''  AVGNA(  1 1  *  AV  GilA \  I )  / AVE 

Avcssd.ii'AvcRsu.ii/AiiE 
cs ii  avcrs *avcrs ( : , i ' /ave 

RFUNC !  i ,  I  >*A7C!S< I ,  I  '/(«V’EFT*AVGN»<  I )  1 
:::o  'r ukc c r , : ) «Avcss»r , i > »‘«:avept*awgi»a<  i > > 

OS  '  DC  AVr I ( 1 1 *hTAN i RFUNC i 2 , 1 ) /RFUNC (1,1)1 

O’OO  AVRFUNi  I ) *S5RT i SF'JNC ( 1 , 1  )«RFUNC (1,1) «RFUNCI 2, 1 )  ‘KFUNCI 2 , 1 ) ) 

OVD  !H5R;i;i)*iAVEPT»»tfGrt*U)i 

0-20  C R H <  I > *SSRT ( AVCRS i 1,1) •AVCRS ( 1 , 1 ) «AVCtS (2, 1 )  * AVCRS ( 2, I ) ) 

0730C  ESH;i)sAVRr  'uN(D*?.3/(  (2‘PI*AVrR0(  I)  14*2) 

v?40  114  CONTINUE 

C.--OC 

0'"v0  DO  203  J-l  ,NUi1 
O’"0  K*J 

OT50  USITE (O', 350)!:, AWFRQ(J), CRH ( J) ,NORN( J) ,AVG«A( J) , AVERT, AVRFUNU) 

OT’O  203  CONTINUE 
1000  DO  210  1*4,  (NUN-41 

.  o ivtusi i > *r<* < o *«vfus(i-2 > ♦: *o  »avfu* <  1- n >avfur( i i 

1004  S*2*Q*AVr  Ur.  i  I  *1  i  *0 1 AVFUA ( I  »2) ) 

i :  0  NORM ( I ) *£)■» (0 iNORIli  1-2)  *2*O*N0RR(  1-1 ) •NORIM  I ) ♦2»0»N0««(  I » 1 1  *0‘NOI(H( I  *2) ) 

I OCOCRNI I )  »3 >  ;  0-C?.N( I -2)  * 2 '0 >CRN(  I- 1 )  *CRH(I ) *0  *2  *CRH (1  *  I  )»0*CRH(  I  *2 ) ) 

1030  210  CONTINUE 

i::ic 

1 032  DO  214  1*4, (NUN- 4) 

1 034  AVRFUNi 1 1 sCRKI I ) /MORiK  I i 

1034  ERH( I i *AUSn)iii I)  »9.6/{  *2'RI  >AVFR0( I ) )  >*2) 

1*33  214  CONTINUE 
01  040  DO  443  >0, NUN-3 

01041  AI2RFUI1  (1(1  *0  •  ( 0  <  1 AVRFUIM  N-2)  ♦AVRFU(4(K*2)  MAVRFUIM  K )  *2.  ‘OMAVRFUN 
0;043i;K4l)*AVSrUNiK-1))> 

01043  443  CONTINUE 
i  :■  4S  50  213  1-4 ,  ( NUM-4 ) 

10CCC  NT.  ITti03,32?)  I  ,hVFRQ(  I),AVRFUNd),ERHd),CRHd),fOOO*NORNd),t  OOOOOOOO* AVFURf  1 1 
0 ' 132  Nf.;TE.:'i,iiO)AVFSOiI),IOOO.*NO«H(I) 

0  734  NR ITc I J ; , 1 00 ( AUF RO ( I ) , 1 OOOOOOOO . » AVPUR ( I ) 
i w 213  CONTINUE 
•  -*3  D’j  0*7  i. *3,  •  NUfl-8  * 

1 :  »3  .RITE .  1 ) ,  I  DO  >»'.'? RO  I K ) ,  AVRr  UN  ( K ) 

o;.*  •>**  cont;., u: 

1 ,  u  »■*"  r  vi.iIm  >  1  ,*,  1*  ,f.'  iJ,4v,r?.«,4X,r?,2) 

. .  * ?  S3-  FCRiiaT .  V> 

1 2 ~CC  333  r0f..iitTvI3, 4* .rO.S,  4A.F7.2 , 4X,F9 .2, 4X,F»  .2, 4X,F»  .2) 

I’.COC  *00  r0F..i,T.I4,4*,F9.3,4A,F?.7) 
too  for.iat..'' 

■  ..0  ico  *  or.  ,'r*  7 1  v » 

I  0  ror.;...:.:2,4»,r3.3,4X,r7.3,4X,F9.J,4»,F9.1,4X,F».J,4X,F7.3) 


Copt  available  to  DTIC  does  not 
pennit  fully  legible  iepioductioa 
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L  COHERENC 


OOO!  P4RAIE7EI  1*  121,  U*  121,  S*  1210 
001 0  INTEGEI  l,F,NAU 
0020  INIE5E*  «U«,jK/{ 

0040  P!*3.  14  I5P3 

0050  Q«. 33333 

0061  NUH*0 

0045  REAL  UN 

0030  REAL  CRH(V),NO«N<U> 

0099  IIMENSION  ERH(W),4V«FUN(V) 

OtOO  DIMENSION  DOT 4(2, V), EOT 4(2, U) 

0110  DIMENSION  BULK  2,  SO 
0120  M  NEKS  ION  4VCRS(2,V) 

110  DIKE  NS  ION  PU«SPEC(Ut,PUS<l/>  ,G4ASiV>  .CIQSPOR  (U )  ,FI  (U  l 
0140  DINENS10N  FUSPEC(V)  ,RFS (V 1 ,CHRENCE ( V) , FREQ <Y 1 
0130  DIMENSION  GMASPEC  ( 0 )  ,C*DSPEC  12, 11)  ,I<FUNC  (2,0  > .  IF  SPEC  (V 1 
0140  DIMENSION  AVRFUNC (U ) ,AVFRO<  U>  ,AVF  I<  11)  ,4VPUI(  II)  ,AOGIA(U) 
0170  DIMENSION  AVCSPUIt(U)  ,AUCQN< II) 

0100  HUN*I  ./FLOAT  INIIM) 

0190  PRINT,  'NIN1IEI  OF  PICESSEI  TER  IS  TO  IE  AVEIAfiEI/ (A«E>? )-' 

0200  READ,AVPT!N 

0210  AVE*«VPTIN 

0220  EM>1  ,/FlOtT  (AVE ) 

0  230  SUSP  T»  0.0 
024 OC 

0230  00  I1 1  K-1  ,S 

0240  REAK  14,200)  J,  UH1,  Ml,  SI  CO 

0270  PRIM- PUT 

0210  SUhPT  •SUMPI+PRIN 

0270  PRTN*PNT 

0300  B»L>II1,K)«PR!N 

0310  IUUI<2,K)-IHT 

0320  111  CONTINUE 

033 OC 

0340  DO  120  J*1  ,NUM 
0350  A0RF0NUI*4.0 
0340  A4FI<J>«0.0 
0320  AVPLIR(J)-0.0 
0380  AUH1A<J)*0.0 
0390  AVC3PUR<J)*0.0 
0400  AUCRS4  I,  JX.O 
0410  «UCRS(2,J)«0.0 
0420  AUCDN(J)*O.Q 
0410  124  CONTINIE 
0440C 

0450  AVEPUSUOPT/S 
0440C  PRINT, AUEPT 
0470  DO  115  F*1  ,AVE 
480  R4U<(F-1)«NUM 
04  IOC 

0500  DO  111  N*1,NUM 

0510  OATAO,K)*IUIK11  ,<!*»►«)  • 

0320  DAT *(  2,i;i»0.0 

C530  E4TA1 1  ,K)*MJU(  (2 ,( IAII+N)  1 

0S40  EAT4(2,K  1*0. 0 

0550  1  13  CONTINUE 

0540C 

0570  C4LL  FOURT (DOT N, NUN,  I ,  "1 ,0 ,0 i 
0580  CAU  F0Uir(E4T4,Ufl,1, -1,0,0) 

0390C 

0400  DO  1 22  t«1  ,NUM/2 

0410  SM4SPCC<I)«(EArA(l,IHE4T4<  I ,  I»  •ENTA<  2. 1)  *€4T4<  2, 1)  I 'NON 
0420  PUNSPECl 1 1*1  DATA! 1 , I XDAT 4<  1 , 1 )  *04T 4(2.1'  *1)41 At  2, 1 1  l/NUM 
0430  CROSPECd  ,1)*(DAIA(1 ,1)«E4T4v  1 . IJ  *0*T«\  2, 1>*E»T4*  2, 1)  1  .'NUN 
0440  CR0SPECC2.  II  >  ( OA  T  A  C  2 , 1 7  «C  AT«<  I,  I  >  -OAT  «<  1 , 1  )•£*  IA<2, 1 )  >  /«»« 
0430  PUSPEC II ) •PUR'JFEC ( 1 1 
0440  122  CONTINUE 
0670C 

OotO  DO  211  3*1, NUM 
0490  «*J-I 

O/OU  F  4C0  ( K 1  »CL  OA I  (  « )  /F  10 4  F  <«#M  I 
710  211  CONTI  HUE 
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COHERENC  (cont'd.) 

0721  DO  117  0*1  *1*1/2 

0734  CC«S4RT(  (CUOSFECd  ,1)  1*(CU01P£C(1  ,D>>  ♦  <  CROSTECI  2,01  )>MCR0SP£C<2,J) )) 
0733  AVCRS<1fB)«AVCRSd,#)»E»A«XI 
0734  AVPUR  ( I)  *AVP4R  (D)  +  EV***  (PURSPECI  >)  ) 

0737  AVGHA<8>«4*B14<0)+£V4*MGIIASP£C<D)> 

0774  117  CONTI  IDE 
0773  113  CONTINUE 
0777  SUM *0.0 
0734  DO  344  K«1  ,NUM/2 

0783  4VCR3C  2,1)  -SORT  (AUGfW(K) )  tSUR  1<  «VFUI<  ( )  1 
0794  CRH  (K)  MUCUS  (1  ,K )/ WCR$<2,K) 

0793  CRH(K)  «CRU(K)*CUH(KI 

0793C  FRIM7,fREQ(K),CRU<lt)  ,AUCUS(1  ,K),A4CRS(2,K> 

0804  URITEI08, I40)FREQ IK)  ,AVPUR (K)  ,4U6HA<  I)  ,AVCRS(2  ,K),  CRIH I) 

0803  SUN«SUH»CRH(K) 

0814  344  CONTINUE 
0812  AUEME>SU»*2.  >nilN 
0114  SUM).! 

0814  DO  989  K*1  ,NUH/2 
0118  SU*-SUHCR8<()*CRU<()-AUEIAEE>AUeUAtE 
0824  989  CONTINUE 
0822  S7DEf«SU*»2.*41UN 

0124  PRINT,  'COHERENCE  IEAN*  ‘ .A4EIASE , "  COIEIECE  ST4NI.  DEV.«/,STDEV 
0123  140  FORMAT  ( 3d  X,Ft  2.9) 1 
0134  240  FORIAT(I) 

0134  STOP 
0844  Eld  * 
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«L  FFT 


1000 
1330c 
1300c 
1370c 
1380c 
1370c 
1400c 
1410c 
1420c 
1430c 
1440c 
1430c 
1  440c 
1470c 
1480c 
1470c 
1300c 
1310c 
1320c 
1330c 
1340c 
1330c 
1340c 
1370c 
1380c 
1370c 
1400c 
1410c 
1420c 
1430c 
1440c 
1430c 
1440c 
1470c 
1480c 
1490c 
1700c 
1710c 
1720c 
1730c 
1740c 
1730c 
1740c 
1770c 
1780c 
1790 
1800 
1810 
1820 
1830  I 
1840 
1830 
1840  2 
1870c 
1880c 
1870c 
1900 
1910 
1720 
1730 
1740 
1730c 
1740c 
1770c 
1780c 
1770c 


SUIR0UT1NC  FOURTdATA.Nh  ,40111,  ISI5N,  IF09I1,  UO»K  > 

bounded  by  5»2»»<-b) *aun( fact«r( j)*«t .3) ,  where  6  ti  the  nunber 
of  bitt  in  the  floating  point  friction  and  factor(j)  am  tho 
prino  factor!  df  ntot. 

progrin  by  aornan  breaner  fron  tho  baaic  progran  by  chhrlea 
radtr.  ralph  altor  suggested  tM  idtb  far  tho  digit  rovoraal. 
nit  lincoln  laboratory,  auouat  1747.  thia  ia  tho  faatoat  aad  neat 
voraatilo  version  of  tbo  fft  k«o*n  to  tho  author,  ahortor  pro- 
grana  four l  and  four2  reatrict  dintnsian  lengths  to  oouora  of  two. 
at*--  1000  audio  tranaactiona  (.nine  1747),  special  laiuo  on  fft. 
aoo  nath  explanation  in  pdp-13  docunentation 


tho  diacroto  fourior  Irani  for*  places  throe  restriction!  upon  tha 
data. 

1.  tho  nunbor  of  laput  data  and  tho  ounbor  of  traoaforn  values 
nuat  bo  tho  aano. 

2.  both  tho  input  data  aad  tbo  traaaforn  valuea  oust  roproaont 
oquiapacod  point!  in  their  respective  donaioa  of  tine  aad 
frequency,  calling  theae  spacing!  dottat  and  doitaf,  it  nuat  bo 
true  that  doltaf>2opi/(an( il*deltat>.  of  courac,  del  tat  mod  aot 
bo  tho  aano  for  ovory  dinonsioa. 

3.  conceptually  at  loaat,  tho  input  data  and  tho  traaaforn  output 
roproaont  ainglo  cydoa  of  periodic  functiooa. 

exanplo  1.  thrre-dinens ional  foruard  fourior  transforn  of  a 

conplox  array  dinonaionod  32  by  23  by  13  in  fortran  iy. 

dinonaion  data (32, 23, 13) ,uork (30) ,on(3) 

conplox  data 

data  nit/32,23,  13/ 

do  I  i*l,32 

do  I  j-1,23 

do  I  1»I,I3 

1  datali  ,  j,h  >*conple*  value 

call  fourt(data,on,3,-l ,1 ,»ork) 

oxanplo  2.  ooo-dinonaiooal  forward  transforn  of  a  real  array  of 

length  44  in  fortran  it. 

dmooaion  data(2,44) 

do  2  1*1,44 

data! 1 , ilnroal  part 

2  data(2,i>*9. 

call  foort (data, 44, I , -l ,0,0) 

dinoaaioa  data< 1 1 ,on( 1 1 ,ifact (32) , »ort (I) 
connon  /nan/  non 
t oop i*4.283 1 83307 
if  (isdin-l  )920, 1,1 
otot*2 

do  2  idin«l , odin 
if ( on (idin) 1720.720,2 
otdt'niot *on( idle) 

nain  loop  for  each  dinonaioo 

np1*2 

do  719  idin>l  ,odin  , 

»*no( idin) 

np2*npi *n 

if (n* 1 1730,790,3 

factor  o 

doconpoao  input  into  n  discrete  tine 

aanploa,  n/2  odd  and  n/2  oven,  keep  deconpoaing  until 

n/2  contain!  too  aonplta. 


FFT  (cont'd.) 


2010  s 
2020 
2030 
2040 
2030  10 
2040 
2070 
2000  II 
2000  12 
2100 
2110 
2120  20 
2130  30 
2140 
2130 
2140  31 
2170  32 
2180 
2100 
2200 
2210  40 
2220 
2234  30 
2240  31 
2230 
2240  40 
2270c 
2280c 
2200c 
2300c 
2310c 
2320c 
2330c 
2340c 
2330c 
2340c 
2370c 
2380c 
2300c 
2400c 
2410c 
2420c 
2430  70 
2440 
2450 
2440  71 
2470  72 
2480 
2400  73 
2300 
2310  74 
2320 
2330 
2340 
2330 
2340 
2370 
2380 

ISM  « 

2410 

2420  03 

3430c 

2440c 

2434c 

2440c 

2470C 

2480c 


otuo*np 1 

lf*l 

idiv*2 

>4UOt>4/101V 

ironen-idivoiouot 
if (lquot-idiv)30,1 1,11 
if(iren>20, 12,20 
nlwo»otwo«ntuo 
Mi4U4t 
jo  to  10 
idiv*3 

ltuol»n/idiv 
iro««#-idi v.ipuot 
iff iouot-idiv)40,3l ,31 
if(iren)40,32,40 
if»cttif>«idiv 
if *if *1 
•■iquot 
*o  to  30 
idiv»idiv+2 
jo  to  30 

if <iro*)40,3t ,40 
n tuo»»t»o»ntuo 
10  to  70 
ifaet(if)«e 

separate  four  cioo* — 

1.  complex  traaoforn  or  real  traosforn  for  the  4th,  Slh.ftC. 
dtaensiom;' 

2.  real  traasforn  for  the  2«d  or  3rd  dtaensioa.  aethod— 
transfor*  half  tho  data,  supplying  tho  othor  hair  by  con- 
jujato  syanotry. 

3.  roal  transfora  for  tho  lot  diaeosioa,  n  odd.  aothod-- 
tramforn  half  tho  data  at  tach  state,  supplvmj  tho  other 
half  by  conjugate  svmetry. 

4.  roal  transform  for  tho  1 s t  dinension,  n  oven.  aethod-* 
tramfora  a  conplea  array  of  length  *72  uhose  real  parte 
are  tho  even  nunbered  real  valvei  and  ahote  taaoinary  parts 
are  the  odd  nuaborod  roal  valuta,  separate  and  supply 

the  second  half  by  eon.iujate  svwietry. 

*on2*np I • ( np2/ntuo ) 
icoseet 

a  f  4 idin-4  >71 , 90 ,90 
i f  < i forn  >72,72 , 10 
icas»*l 

if ( idin-l 173,73 , 10 
icate*3 

if (ntuo-npl >90,40,74 
icate«4 
»tuo«nlue/2 
•••/2 
no2*np2/2 
ntot*atot/2 
i«3 

do  80  j*2,ntot 
data! j  ><data(i) 
i1Mj*npl 

tf( icase-2> 100, 93,100 
1 1  rn’j.npO*!  1  *nprev/2) 

shuffle  on  tho  factors  of  tvo  in  n.  as  tho  shufflinj 

can  bo  done  by  smple  interchange,  no  «orkin-;i  array  is  oeodod 

doconposod  input  nust  bo  shufflod  to  take 

cart  of  bit  reversal  to  output. 

oxanole  f0.f4,f2tf4.fl,f3.f3,f7  .ill 


175 


2470c  product  coefficient!  -pO , « 1  ,j2,g3,j4, $5  ....jn 

2709c  fourier  cm  thea  bo  bona  it  place  eithaut 
2710c  additional  cart. 

2720c 

2730  100  if<atvo-npl 1400,600, 110 

2740  110  np2hf*np2/2 

2730  j»l 

2740  do  ISO  i2»l,ap2,aon2 

2770  if < j-t2J 1 20, 130,130 

2780  120  i1«aa»i2*nan2-2 

2700  do  123  il»i2,ilnax,2 

2(00  da  123  i3*i I ,ntot ,»p2 

2810  j3*j*i3-i2 

2820  tenpmdatat 13) 

2830  tenpi*data( i 3* 1 ) 

2840  datal i3)*data( j31 

2830  data(i3*1  >*data< jl*1 ) 

2840  data< j3)*te*pr 

2870  123  datat j3*1 >»teapi 

2880  130  **ap2hf 

2800  140  if < j~n) 130 , ISO, 1 43 

2000  143  j* j-n 

2910  »*n/2 

2920  if(«-nan2)l30,l40,l40 

2030  130 

2940c 

2930c  aain  loop  far  factor!  of  tut.  ptrfara  fountr  transform  of 

2040c  length  four,  uith  ant  of  leagth  t«o  if  atadtd.  the  teiidle  factor 

2070c  u*txp<  iiigne2«pi*sqrt  (-1  )*a/(  4anaai)  I .  chect  for  u*itign*!grt(-1 ) 

2080c  and  repeat  for  mu lgnaaqrtl-l  l*coa.io-j«te(»>.  aote  j»s«irt(-l> 

2000c 

3000  non2t>non2«non2 

3010  ipar>atao/np1 

3020  310  if ( ipar-21330,330,320 

3030  320  ipar*ipar/4 

3040  go  to  310 

3030  330  do  340  il«l,itrng,2 

3040  da  340  j3»il ,aoa2,npl 

3070  do  340  kl3j3,ntot,nan2t 

3080  k2*kl*non2 

3090  tanpr*data(L2) 

3100  teapi»data(k2*t 1 

3110  data(t2)*dala(kl >-lr«pr 

3120  data(k2*l >*data(l I  *  1 >-le*pi 

3130  datalkl  )»data<k  I  lUanpr 

3140  340  data! k 1*1  )>data(klO)»tenpi 

3 1 30  330  nnai*aoa2 

3140  340  if  1imi>:-np2hf >370, 400,400 

3170  370  lnax*na>:0(aoa2t.anaK/2) 

3180  i  f (nnax-nan2) 403, 403,380 

31*i>  380  thata**tuopi  >flaat  <  aoa21 /float  <  4*««<a>:> 

3200  if  1 1 1 ign ) 400,310,310 

3210  390  thetan-theta 

3220  400  ur*co!( theta) 

3230  xiMinl  thela) 

3240  »itpr>-2.*ui >*» 

3230  uitpi*2.eur*«i 

3260  403  do  370  l*nan2rl*a>, »oa2l 

3270  n»l 

3280  if ( enai'ton2) 420,420,410 

3270  410  u2r<ur •vr-ui»»i 

3300  u2i<2.>«r>ui 

3310  u3r*u2raur-*2i «*i 

3320  »3i«a2r<«i«u2i«ur 

3330  420  do  *30  Ifafjllrng;? - 

3340  do  330  ,i3*i  1 , nan2,apl 

3330  tnin«,i3*ipar*n 

3340  if(nnax-aon2)430, 430,440 

3370  430  k«ia«j3 

3380  440  kdif *ipar«nnai 


FFT  (cont'd.) 


3390  430  kst*p«4akdif 

3400  4o  320  kl«k«U,tt*t,k»t*p 

3410  k2»k1*kdif 

3420  k3*k2»kdif 

3430  k4»k3«kdif 

3440  i f ( «*a«-no*2) 440,140, 480 

3430  440  ulr*data(k1 >»data(k2) 

3440  ulladata(k1«1)*data(k2*l> 

3470  u2r*data(k3)rdata<k4> 

3480  u2i*d*ta(k3*l I*data(k4*l1 

3490  u3r*data(k1 )-data(k2) 

3300  u3i*data(kl«t)>data<k2*1) 

3310  if<isijnl470,475,473 

3320  470  u4r>data(k3*1 >-data(k4»l I 

3330  u4i»dala<k4)-data<k3> 

3340  9k  ta  310 

3330  473  u4r«data(k4«l>-data(k3*U 

3340  u4i*data(k3)*data(k4) 

3370  jo  ta  310 

3380  480  t2r>w2radatatk2)-«2i«data<k2»U 

3390  t2i*u2r»dataU2*1  )**2i*dat»(k2> 

3400  t3r*wradata<t3)-ai«data(k3»1 ) 

3410  tSikuradatalkl-M  )«ui*data(kl) 

3420  t4r*«3r*data(k4)-u3i*data(k4»1 ) 

3430  t4i*w3r*4ata(k4«1 )*«3i*data<k4> 

3440  ulr*data(k1 )*l2r 

3430  ■  uli*data<k1H)+t2i 
3440  u2r*t3r*t4r 

3470  «2i»t3irt4i 

3480  ujr«data(kn-l2r 

3400  «3i»dala<kt»1l-t21 

3700  if( iii^n) 490,300,300 

3710  400  «4r»t3i-l4i 

3720  *4i»t4r-l3r 

3730  )0  to  310 

3740  300  «4r»t4i-t3i 

3730  *4i»t3r-l4r 

3740  310  datatkl )>utr«u2r 

3770  datalklrl )«ul i*u2i 

3780  data(k2)*v3mi4r 

3700  4ata(k2r1 )*o3i*u4i 

3800  data(k3)>ulr-u2r 

3810  data(k3tU«u1i-u2i 

3820  datalkl )<u3r-n4r 

3830  320  4ata(k4*1 >»u3i-u4i 

3840  k*i*»4*<k«t»-j31* j3 

3830  kdif’ksttp 

3840  if <kdif *np21 430,530,330 

3870  330  coatinu* 

3880  «»»nax-« 

3800  If(isi9nl340,330,330 

3000  340  to*pr«wr 

3010  xr»*Mi 

3020  *i«-t*apr 

3930  jo  ta  340 

3940  330  t*«pr««r 

3030  ar»ui 

3940  wi*t*npr 

3970  540  if (a-lnax>343,343,4IO 

3080  343  t*npr*wr 

3090 

4000  370  ui*»iawstpr«t*npr*astpir«i 

4010  ipar«3-ipar 

4020  Mlitiuomm 

4030  fo  to  340 

4040c 

4030c  naif  loop  far  factors  not  **ual  to  tso.  apply  tko  taiddl*  factor 

4040c  »Mxp(  is  ljoa2.pi  nqrt  ( -  I )  •<  J2-I  >M.i1-.i2)/(ap2*ifpl  1 ),  than 

4070c  porfora  a  fourtar  traasfara  af  Iaaqth  ifacttiO,  aatiaq  u<*  of 
4080c  coaja^at*  synaatriot. 

4090c  10*  Olffarancot  ta  noth  nata  coapatrr  dacunantatipn 

4100c 


Ml 


FFT  (cont'd.) 


4110  400  if <»lo0'0p2) 405, 700,700 

4120  <45  ifpl»no«2 

4120  if »1 

4140  np1hf>np1/2 

4150  <10  ifp2«ifpl/ifacttif> 

4 1 <0  j1rnj«np2 

4170  if(ica»«-3XI2,41 1  ,<12 

4180  <11  j1rnj*(np2+i7pl 1/2 

4100  j2*tp*np2/if»{t(if 1 

4200  jt  rg2-< j2*tp*ifp2>/2 

4210  <12  j2nin*1*ifp2 

4220  if ( if  p  I -np2) 41 3, 640, 440 

4230  <15  <o  <35  j2*j»nin, ifpl ,  ifp2 

4240  thola«-t»opi*f laitl j  2 - 1 )/fIo»t(n<2> 

4250  if(itifn)<25,<20,420 

4240  <20  thotja-tkoia 

4270  <25  *inth»im<ttiota/2.) 


4280 

•ftpra-2.*iiKth*<inUt 

4270 

»«tpi*> i» ( thoti ) 

4300 

ar»»«tpr+1 . 

4310 

ui'ustpi 

4320 

jl«in*j2*ifpl 

4330 

do  433  jlojlmn, jlr»9,ifp* 

4340 

itnaxaj l*i 1r*o-2 

4350  ' 

— 

-<o-  <30- 1  r=-jT;msrrJ  • 

4300 

do  430  i3*i 1  ,n»ot ,op2 

4370 

j3oa**l3*ifp2-npl 

4380 

do  430  j3«i3,j3«»x,»p1 

4370 

tonpr*datat j3) 

4400 

data! j3)*d»ta< j3)»or- data!  j3*i >»ui 

4410  <30 

datal jJ+l >»tonpr*»i»data(j3*l l»»r 

4420 

tonpravr 

4430 

upaor«wotpr*ui  •a»tpi*i<r 

4440  <33 

oi*to«pr»w*tpi»ol7«»tpr*wi 

4430  <40 

thotaa-tuopi/floall  if  act!  if) ) 

4460 

if!  isijn) 430, 443,443 

4470  <45 

tkota*-th*ta 

4480  <50 

th*ta/2.) 

4470 

u«tpr*-2. •«inth*«ioth 

4300 

mtpi**m(thota) 

4510 

k*top*2*o/ifacl<  if) 

4520 

kran)*k< top*(  if  act!  if  I/2HI 

4330 

da  478  il»l  ,i  ln>1,2 

4340 

df  478  i3»i  1  ,ntot,iip2 

4550 

do  470  kmii»1  ,kranj,kstop 

4360 

j1«ai»i3»jlr»-3-ifp1 

4370 

do  400  jl  “i3,,l  laax,  ifpl 

4380 

j3oa>aj 1 f ifp2*opf 

4370 

dP  480  j3«j1 , J3na«,n»1 

4<00 

j2na»»,i3*ifp1*ifp2 

4410 

k>kaia*<  j3-j1»<  jl*i3)/ifact(  if)  1/npltif 

4620 

if (kaia-1 1453,453,445 

<<30 

<33  s««r*0. 

<440 

tuoi«4. 

4630 

do  440  j2*.i3,  j2»ax,  ifpJ 

444  0 

suarMunr  *<al  a(  ,t2 ) 

4474 

660  +d»ta(  j  2  ♦  1  ) 

AitO 

workUMoimr 

4470 

»ork(k*1 )*»uoi 

4700 

jo  to  480 

4710 

44 

3  kcon.l*)  *2  Mn-kmn+l ) 

4720 

i2».i2»a« 

4730 

svnr«data( j2) 

4740 

fum<data<  j2*1 1 

4730 

oldir«0. 

4740 

0 ld% 1 >0 • 

4770 

j2«j2-ifp2 

4780 

470  ttnprMunr 
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FFT  (cont'd.) 


4700  tonpi«su<ii 

<800  aunr*tuoar*aunf-oldsr«data<  i2) 

<010  *oni*twowr*auai-al4ai>data(  j2+1 ) 

4820  oldir>tt«pr 

4030  oldti*to*pi 

4840  j2*j2-ifp2 

4050  if< j2-j3)i?3, 473,470 

4040  473  tanpf«yraauiir-ol4sr+dat»( .12) ' 

4870  ttApi«m<4uni 

4880  york (l )»tonpr-t»npi 

4800  »ork(kcon.i)*ilApr»t«npi 

4000  tt«pr«wra*u»i-oldsi*data<  j 2*  t ) 

4010  tt*pi*aia«unr 

4020  york <k«1 >«t*«pr«ttnpi 

4030  uorMkeonj+l  >*t»npr-t»«pi 

4040  480  continue 

4030  if (knin-l 1405,405,484 

4040  483  yr*ustpr»1. 

4070  ui*u*tpi 

4000  jo  to  400 

4000  484  t*«pr*wr 

5000  yrawravitpr*wi*a«tpi*yr 

3010  ui*toapr*u«tpi*ai<uitpr«wi 

3020  400  tuawi-ur+ur 

3030  if (ic>f«-3)402, 401,402 

3040  401  if(ifpl-np21403, 402,402 

5050  402  k*1 

5040  i2n»x*i]*np2'npl 

3070  Oa  403  i2*i3,i2nai,ipl 

SOSO  dita(i2)*york(k) 

3000  data(i2+l >*uork (k*l ) 

3100  403  k«k«2 

3110  $a  ta  400 

5120c 

3130c  coaplota  a  rail  traaaforn  la  th«  1st  dmoaaioa,  n  odd,  by  coa- 
3140c  judato  aynaotrita  at  »ach  ataio. 

3130c 

3140  403  .i3nax«i3*ifp2-apl 

3170  do  407  j3*i3, .ilnax.apl 

3180  j2aai»j3*np2-j2atp 

3100  do  407  j2*j3,.i2«ax, ,i2atp 

3200  jlaa«*j2*.|lr<|2-ifp2 

3210  j1cnj«,i3*,j2n»:i*.j2stp' j2 

3220  do  407  j1*j2,.tl«»x,ifp2 

3230  k*l ♦ jl • i3 

3240  data! j I >*uork (k ) 

3230  data!  j  1  •  U*»ork(k  *1 ) 

3240  If(j1-j2)407,407,404 

3270  404  datal j lc nj) *work Ik) 

3280  datal.ilcnjal  )«-work(k*1 ) 

3200  407  j lcnj*j 1cnj-ifp2 

3300  408  conti nua 

3310  if»if*l 

3320  ifp1<ifp2 

3330  i  f ( ifpl -apl >700,700 ,410 

3340c 

3330c  coaplatt  a  raal  traaaforn  m  th»  lat  dinonnon,  »  rvrn,  by  con- 
3340c  jufatt  aynaotrita. 

3370c 

3380  700  ]0  to  (000,800,000,701), icaat 

3300  701  nkalfM 

3400  n«n*n 

3410  tfi«ta*~taopi/ float)  n) 

3420  if (iaidn)703,702,702 

5430  702  th«ta*-Uota 

3440  703  ainU»ain(tb«ta/2.) 


3430 

aalpf»-2.*ainth*ainth 

3440 

aatpi *» ia( thpta) 

3470 

«r«ystpral . 

3480 

ui'ratpi 

3400 

i«ia«3 

3300 

jain*2an(ialf‘l 

3510 

fo  to  723 

n?0  710 
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FFT  (cont'd.) 


3330 

da  720  i*iflin,atot,np2 

3340 

4unr«( datai  i  )«data< j  >1/2. 

3350 

«u«i*(daU(i*t  )«diU(  j*1  ))/2. 

3340 

difr«< datai  i> -datai  j )  >/2. 

5370 

dif i*< datai i*l >-Jata<  ,i*l 11/2. 

53B0 

taxpr»wr»»u«i*ui • d i fr 

3390 

ttnpi*ui*iu*i-ur*difr 

5400 

datai i  >*»unr*ta»pr 

3410 

data(i»1  >»dif  i*tmpi 

5420 

data! j >  »iu»r- tanpr 

3430 

■  J  a 1 a  <  j ♦ 1 >»-difi*ta«pi 

5440 

720 

j*j*"p2 

3430 

5440 

jnin»j»i«-2 

3470 

ttapr«ur 

5430 

•r*wr*nttpr-«i  >mtpi  »ur 

5490 

»J*ttxpr*M»tDl*ul*»stpr*ui 

3700 

723 

if(i«ic-jnin)710,730,740 

5710 

730 

if  fi*i9J>l?31 ,740,740 

3720 

731 

da  733  I*lnin,ntat,np2 

3730 

733 

datai i*1 l»-datal  1*1 1 

3740 

740 

>p2*ap2*ap2 

3730 

ntot*nlot *atot 

5740 

j»ntol*1 

3770 

iaa»*ntat/2+l 

3780 

745 

mit*inax-2«nf.alf 

3790 

i*inin 

3300 

$a  to  735 

3810 

730 

datatj)<data(i) 

3020 

•data <  j ■>  1  )«-dala< i«t ) 

5830 

733 

i«  i*2 

3840 

j»j-2 

3830 

if  ii-iiian)  730,740, 740 

3840 

740 

datai j>*datai ini n) -data! i»ii+l) 

3870 

data! j»l 1*0, 

5880 

if( i-j)770,780,710 

5890 

743 

datai  .i)*dataii> 

5900 

data) j+l >»data( i*l ) 

5910 

770 

i*i-2 

3920 

j»j-2 

3930 

if ( i-ia la  1773, 773, 743 

5940 

775 

data! j)*data<  mint*datal ) 

3930 

data< j*1 )>0. 

3940 

i»ax*inia 

3970 

90  to  745 

3980 

780 

datail >*dataill*datai2> 

3990 

data(2)*0. 

4000 

4010c 

9a  to  900 

i020c 

conplfto  a  raal  traaaforn  for  tn« 

oOJOc 

conjuaait  tynaatriat. 

4040c 

4050 

800 

if  ( i  1  rn<|-tipl  1103, 900  ,f00 

4040 

803 

da  840  i 3* 1 ,ntot,np2 

4070 

i2«»i»i3*np2-np( 

4080 

do  840  i2*i3,i2nax,npl 

4090 

inin*i2*i 1rn9 

4100 

i»a»*i2»npl-2 

4110 

jnax»2«»3*npl-ixjit 

4120 

if (i2-i3 >820, 120,810 

4130 

810 

,inax*jnai«ap2 

4140 

120 

if ( 1  din* 2) 850,830,330 

4150 

830 

j>.ma<*ap0 

4140 

do  840  i>imn,iaa>,2 

4170 

datai i)«data( j) 

4180 

dataiirl )«-data<  j  ♦  1 ) 

4190 

840 

j*j-2 

4200 

830 

j*.i»a« 

4210 

da  140  i»inia,ina::,np0 

4220 

datai 1 >*data( j ) 

4230 

dataii»ll»-dala(.i«1 ) 

4240 

4230c 

340 

,l*,i-np0 
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FFT  (cont'd.) 


4260c  tnd  of  loop  oi  each  dinontion 
4270c 

6290  900  Bp  0 » n p I 

6290  npt*np2 

4300  910  nprov*n 

4310  roturn 

4320c 

6330c  :wt  coot  htr#  only  i«  orror 
4340c 

4330  920  writ#  (6,  9101)  nOin 

4360  rttufn 

6370c 

4310  930  uritf  (6,  9602)  ldin,  no  i idin) 

6390  -  ■  rtturn  • 

4400c 

4410c 

4420  9401  fornat  (3</> ,  10('  C»>' )  ,5*.  't  o  <1  t  i.  n  0  l", 
44301  '  o  x  o  c  u  t  »  d' , I Ox , 'only' , i2, '  dinoniion.  ) 

4440  9602  fornat  <3</) , I0< '<*>')  ,Jx, 'f  t  u  r  t  to  l',3t, 
64300  '»  *  t  c  u  t  o  d' , lOx, '  len-iUi  of  dinticiot  I',i2, 

64600  '  it  only'.M) 

6470  ttd 
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